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1.0 INTRODUCTION 

The objective of this research was to develop analytical tools 
capable of economically evaluating the cyclic time-dependent plasticity 
which occurs in hot section engine components in areas of strain 
concentration resulting from the combination of both mechanical and 
thermal stresses. The techniques developed are capable of accomodating 
large excursions in temperatures with the associated variations in 
material properties including plasticity and creep. 

The overall objective of this research program was to develop 
advanced 3-D inelastic structural/stress analysis methods and solution 
strategies for more accurate and yet more cost-effective analysis of 
combustors, turbine blades, and vanes. The approach was to develop a 
matrix of formulation elements and constitutive models, three 
increasingly more complex formulation models and three increasingly 
more complex constitutive models. 

The three constitutive models were developed in conjunction with 
optimized iteration techniques, accelerators, and convergence criteria 
within a framework of dynamic time incrementing. These consist of a 
simple model, a classical model, and a unified model. The simple model 
performs time-independent inelastic analyses using a bilinear 
stress-strain curve and time-dependent inelastic analyses using a 
power-law creep equation. The second model is the classical model of 
Professors Walter Haisler and David Allen (Reference 1) of Texas A&M 
University. The third model is the unified model of 



Bodner, Partom, et al. (Reference 2) . All models were customized 
for a linear variation of loads and temperatures with all material 
properties and constitutive models being temperature dependent. 

The three formulation models developed are an eight-noded 
midsurface shell element, a nine-noded midsurface shell element and a 
twenty-noded isoparametric solid element. Both of the shell elements 
are obtained by "degenerating" 3D isoparametric solid elements and 
then imposing the necessary kinematic assumptions in connection with 
the small dimension of the shell thickness (References 3 and 4) . The 
eight-noded element uses Serendipity shape functions and the 
nine-noded element uses Lagrange shape functions. The eight-noded 
element uses Gaussian quadrature for numerical integration, with 
nodal and surface stresses being obtained by extrapolation/mapping 
techniques. Lobatto quadrature is being used with the nine-noded 
element to effectively provide for direct recovery of the stresses 
and strains at the surfaces and node points. The eight-noded element 
has an excellent combination of accuracy and economy in the normal 
element aspect range encountered when modeling most hot section 
components. The nine-noded Lagrangian formulation overcomes the 
shear locking problem experienced when the element 
size-versus-thickness-aspect ratio becomes very large. The 
twenty-noded isoparametric element uses Gaussian quadrature. 

A separate computer program has been developed for each 
combination of constitutive model-formulation model. Each program 
provides a functional, stand alone capability for performing cyclic 
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nonlinear structural analysis. In addition, the analysis 
capabilities incorporated into each program can be abstracted in 
subroutine form for incorporation into other codes or to form new 
combinations. These programs will provide the structural analyst 
with a matrix of capabilities involving the constitutive 
models-formulation models from which he will be able to select the 
combination that satisfies his particualr needs. 

The program architecture employs state-of-the-art techniques to 
maximize efficiency, utility, and portability. Among these features 
are the following: 

(i) User Friendly I/O 

e Free format data input 

e Global, local coordinate system, (Cartesian, Cylindrical, 
Spherical) 

e Automatic generation of nodal and elemental attributes 
e User-controlled optional print out 

J 

Nodal Displacements 
Nodal Forces 
Element Forces 

Element Stresses and Strains 

(ii) Programming Efficiency 

r 

e Dynamic core allocation 
e Optimization of file/core utilization 
e Blocked column skyline out-of-core equation solver 
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(iii) 


Accurate and Economical Solution Techniques 
e Right-hand side pseudoforce technique 
e Accelerators for the iteration scheme 

e Convergence criteria based on both the local inelastic 
strain and the global displacements. 

The ability to model piecewise linear load histories was also 
included in the finite element codes. Since the inelastic strain 
rate could be expected to change dramatically during a linear load 
history, it is important to include a dynamic time- incrementing 
procedure . 

Three separate time step control criteria are used. These are 
the maximum stress increment, maximum inelastic strain increment, and 
maximum rate of change of the inelastic strain rate. The minimum 
time step calculated from the three criteria is the value actually 
used. Since the calculations are based on values readily available 
from the previous time step, little computational effort is required. 

These formulation models and constitutive models have been 
checked out extensively against both theory and experiment. Figure 1 
shows the correlation between Bodner's model in the eight-noded and 
mid-surface shell element (MSS8) and both experiment and other 
predictions (Reference 5) . Figures 2 through 6 illustrate the 
predictability of the classical Hiasler-Allen model. Figure 7 shows 
a comparison of both Bodner's model and the simple model to both 
experiment and independent predictions (Reference 6) . 
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These nine programs, both source (Fortran 77) and compiled, have 
been installed and checked out on the NASA- Lewis CRAY-1 machine. The 
interactive deck generator has been installed on the NASA-Lewis 
AMDAHL machine. 

Table 1 shows the lines of source code for each of the nine 
computer programs. These numbers do not include the interactive deck 
generators. 



Table 1. Lines 

of Source 

Code 


[ — 
I 



Elements 


1 


20-Noded 

8-Noded 

9-Noded 

i 

Simple 

8300 

13.800 

17 . 900 

1 

Constitutive 

Haisler-Allen 

9200 

16.300 

19 . 000 

Models 

Bodner 

7300 

1.3,800 

17 . 600 


Since these programs use dynamic core allocation, they can be 
[— recompiled to size for any specific machine. They are presently loaded 

for 10 7 bytes of core. At this size, the maximum problem would be 

approximately 4000 nodes and 1000 elements, and 24000 degrees of 

I 

freedom. 

I 

r 
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Figure 1. Displacement Controlled Cycling Results. 
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Figure 2. Load Histories for Plasticity Example. 
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Figure 7. Comparison of Bodner Mode to Simple Creep. 


2.0 TECHNICAL PROGRESS 

The first activity in this program was the performance of a literature 
The pertinent results of this are given in Appendix A. Based on 
the results of this survey, three constitutive models and three 
formulation models have been developed. 

The final versions of the computer programs for the 3D Inelastic 
Analysis Methods contract have been installed on the NASA-Lewis CRAY. 

There are nine separate programs, each with a different combination of the 
3 element types and 3 constitutive models. The element types are a 
20-noded insoparametric brick element, an 8-noded midsurface shell 
element, and a 9-noded midsurface shell element. The constitutive models 
are a simple isotropic hardening classical plasticity model, a 
sophisticated classical plasticity model with combined kinematic and 
isotropic hardening, and a unified model. Both classical plasticity 
models also have a secondary creep model combined with them. 

The major features of the programs are described below. These 
features are generally available in each of the programs except as 
specifically stated. 

The physical geometry of the finite element models is described by 
nodal locations and element connectivities. The shell programs provide 
for automatic generation of nodes and elements if requested. Since 
quadratic shape functions are used in all the element types straight or 
curved boundaries that can be described by a quadratic are easily modeled 
simply by appropriately placing the nodes on the element edge around the 
boundary. No special information is necessary. 
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Loadings can be specified in the form of nodal concentrated loads, 
displacements, and/or temperatures or on an element level as 
pressures for the 20-noded brick element and as pressures and line loads 
in the shell program. Global loadings such as gravity and rotational 
loads can also be specified. Loads with the same value can be easily 
applied to large numbers of nodes and elements since series of evenly 
incremented node and element numbers can be specified. 

All of the programs provide for linear ramping of loadings. Initial 
and final conditions for the load case are input and the intermediate 
conditions are obtained by linear interpolation. The intermediate points 
can be specified in a couple of ways. A number of even load increments 
can be specified with or without a reference to time. Also available is 
a dynamic time incrementing scheme which allows the program to calculate 
time increments and that user input maximums of stress change, inelastic 
strain change, or inelastic strain integration error are not exceeded. 

In this way the program will use large time steps for load intervals 
causing small inelastic action and small time steps for loadings causing a 
great deal of inelastic action to occur. 

The constitutive models are capable of predicting inelastic responses 
for constant or variable temperature conditions This is done by including 
terms that arise due to thermal variation in the derivation of these 
models and by allowing the user to specify the required inelastic data at 
up to 10 different temperatures per material type. There can be up to 3 
different material types per model. Each of the constitutive models is 
capable of modeling response due to simple loading as well as reverse or 
cyclic loadings. 
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The elastic material properties (elastic modulus, Poisson's ratio, and 
thermal expansion coefficient) can be specified at up to 10 different 
temperatures per material type for orthotropic as well as isotropic 
materials . The program will then linearly interpolate between the values 
at the given temperatures. 

Several numerical techniques have been included in the programs in 
order to speed up execution time, convergence, and make efficient use of 
available memory. All of the programs provide for out of core solution of 
the system of equations so that larger systems can be solved. The memory 
available for solution is determined by the program and the system of 
equations is broken up into appropriately sized blocks, stored on file and 
solved a block at a time using the available core. The shell programs 
also use a dynamic core allocation scheme which makes maximum use of the 
available core. This is done by the program surveying the input for the 
problem to be solved and reserving just the amount of core needed to store 
the required information and still have a large sized block of core 
available for solution of the equations. Information will be stored on 
file as necessary in order to leave a solution block of core available. 

The shell programs also have a provision for the user to specify nodes 
which are fixed (with zero displacement) throughout the problem. The 
equations associated with these nodes are then eliminated from the system 
of equations, thus reducing the problem size. 

An Aitken acceleration scheme which causes a more rapid convergence to 
the inelastic solution is built into the 20-noded brick and 8-noded shell 
programs . This scheme uses the two previous iteration solutions to modify 
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the present iteraction toward the apparent converged solution. This 
scheme is used with the plasticity and unified constitutive models, but 
net with the seconday creep model. 

These programs cover a large spectrum in sophistication in their 
ability to predict inelastic response. The problem to be solved, the 
desired degree of accuracy and available amount of effort will dictate the 
element type and constitutive model combination to be used. A great deal 
of overlap occurs, especially for simple problems with simple loadings, 
but as more complexity is added one combination may be better suited than 
another. Many capabilities are available in these programs and can be 
used to great advantage in a wide variety of problems. 


2.1 CONSTITUTIVE MODELS 

Three separate constitutive models have been implemented in each of 
the three formulation models. These constitutive models represent a broad 
perspective in the present state of the art in constitutive modeling, 
running from a simplified model based on the classical plasticity theory, 
to a much more sophisticated and advanced treatment of the classical 
plasticity theory, and finally a unified model where plasticity and creep 
are viewed together as inelastic action. A separate seconday creep model 
is combined with the two classical plasticity models in order to expand 
their capability to the time-dependent domain. 



2.1.1 SIMPLE PLASTICITY MODEL 

Classical plasticity theories propose the existence of a function 
describing the onset of yield and leading to an uncoupled rate-independent 
of deformation. There are four common ingredients in the majority 
of the incremental classical plasticity theories: (1) a yield function 

distinguishing inelastic and elastic deformation, (2) a relationship 
between the stress increment tensor and the elastic strain increment 
tensor, (3) a description, of the rate- independent plastic strain increment 


such as the normality condition, and (4) a work-hardening rule describing - 
the evolution of the yield function under mechanical loading. 

The simplified plasticity model is an incremental classical plasticity 
model which works with effective stresses and strains in the constitutive 
law and expands these effective quantities to component form through 
application of the normality condition. 


The yield function used in the simplified model is the Von Mises yield 
function and is given by: 



where : 


K) = 3/2 S. . S., 



0 


( 1 ) 


s ij = °ij ' 1/3 <° u ♦ °22 * ' 33 ) 

( 2 ) 

is the deviatoric stress tensor, and: 

cem the ®aximum effective stress; that is, the yield surface 
radius . 
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The above shows that the yield function accounts only for an expansion 
of the yield surface about the origin (isotropic hardening only). If this 
function is less than or equal to zero, there is no plastic action. If it 
is greater than zero, plastic action occurs. 

The relationship between the stress state and elastic strain state in 
the simplified model is characterized by an equation relating effective 


stress and elastic strain: 


-E - 


s r 


where : 

F = 


2(ln0 


( 3 ) 


( 4 ) 


is the effective modulus, which differs from the elastic modulus due to 
the difference in the definitions of effective stress and effective 
strain. The effective strain is defined as: 


J 2 


- f [( 


e n - c „) 2 + ( - e 
2 


2 ( C 12 * C 23 


22 ) * ( C 22 " C 33 ) + ( C 33 

)] 1/J 


’11> 


( 5 ) 


2 * ‘ 31 2 


where engineering shear strains are used. 

The effective stress is defined as: 


». * }= [(°11 - 'llf * (°22 ‘ a 32) 2 'll) 

(" 12 2 * ° 23 2 * S 31 2 )]’ /2 


( 6 ) 
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as 


The plastic strain increment is calculated using effective quantities 


Ac 


P 

e 


* AC 

e 



(7) 


where: 

Ae e is the total effective strain increment, and (8) 



is the effective strain increment. 

To expand this effective plastic strain into component form and also 
ensure that the plastic strain increment is normal to the yield surface, 
the Prandtl-Reuss Normality Flow Rule is used. This is given by: 

dcf . = dX S. . ( 9 ) 

ij 


where the scale dx can be calculated from the effective plastic strain 
increment. To see this, rewrite the Flow Rule as: 



( 10 ) 


( 11 ) 



and: 


2 S ij S ij 


Co,)' 


( 12 ) 


(13) 
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So: 


dc ? = | dX a 
e 3 e 


(14) 


dX a 



(15) 


Finally, the work-hardening rule is fully isotropic as mentioned 
previously. Thus, the maximum effective stress value is retained as the 
yield surface size for use in the yield function. 

For a typical load case, the finite element equations are used to 
obtain displacements due to the structure loadings and previous inelastic 
history. These displacements are used to calculate total strain 
increments. From this, a stress increment is calculated with the initial 
assumption of no plastic action. The stress state is then used in the 
yield function to check whether any plastic action occurs in the load 
step. If not, calculations proceed to the next load case. If plastic 
action does occur, the effective elastic strain ACg is calculated and 
subtracted from the total effective strain increment to obtain the 
effective plastic strain increment. This is expanded to component form 
using the flow rule, and an incremental pseudoforce is then calculated. 
After going through this procedure over the entire structure, the 
pseudoforce increments are applied to the structure with the other loads, 
and the process is repeated until convergence is obtained. 
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2,1 ‘ 2 Ha.isler-Allen Cla ssical Plasticity Model 
Professor Walter Haisler and David Allen of Texas A&M developed their 
classical plasticity theory along the sane lines as other classical 
plasticity theories. However, they have added the capability to 
accurately handle thermodynamic loadings in conjunction with multiaxial 
mechanical loadings. Previous classical plasticity models have not been 
used to model rate-dependent or temperature-dependent media, until 
recently. 

The yield function used in the Haisler-Allen model is assumed to be of 
the form: 



■ o. .) * k 4 
ij 



( 16 ) 


where : 

Sij *= the second Piola-Kirchoff stress tensor 

a ij = coordinates of the yield surface center in stress space 

k = a characteristic radius dimension describing the size of the 

yield surface in stress space 
-P 

de = the history of equivalent uniaxial plastic strain, and 
T = temperature 

The yield surface described by this equation defines an area stress 
space enclosed by the surface where all deformation is elastic and thus 
recoverable, whereas on this surface, inelastic deformation is allowed. 
The form of the yield function suggest that the yield surface may both 
translate kinematically by means of the translation tensor «ij and 
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expand insotropically due to changes in the radius dimension k, thus 
producing a combined isotropic kineatic work-hardening rule. Note that 
any thermally induced changes to the yield surface will be totally 
isotropic in nature. 

Differentiation of this general form of the yield function gives a 
statement of consistency during plastic loading: 







dl ? - 2k |* dT = 0 


(17) 


where the term 
3F 



represents 

3F(S . .-o. .) 

U U 
as. . 
ij 

evaluated at Sj_j- c:^j, which can be seen to be equivalent to 
8F 

d(S . ,-cr. .) ’ 
ij ij 

Since during neutral loading (one in which the yield surface remains 
unaltered) , the plastic strain increment and da^j are zero, it is 
apparent that a statement governing loading can be defined to be: 


ir -»#«>• 


whereas unloading is described by: 


3F . _ 3k . n 

3S“ dS ij " 2k 3T dT " 0 


(18) 

(19) 
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Two commonly accepted yield functions of the above form are the Tresca 
and Von Mises yield functions. The Haisler-Allen model uses the Von Mises 
yield function, which can be described as: 

fOV = \ [(», - c 2 ) 2 . «r 2 - o/ ♦ C°3 * <V 2 1 * k 2 (20) 

where o 2 , and o 3 are principal stresses and k is the current 
uniaxial yield surface size. This yield condition supposes that yielding 
is a function of deviatoric stress only. Utilization of this yield 
function limits the model to analysis of initially isotropic media. 

The stress tensor may be related to the elastic strain tensor by: 

S . . = D?. (z - \ 

i J i J an \ am am am am ) 

( 21 ) 

where E__ is the Green-Lagrange strain tensor, D^. is the 

i jmn 

temperature-dependent elastic constitutive tensor at time t, and the 
collected terms in parentheses represent the elastic recoverable strain. 

In addition, superscripts P, C, and T denote plastic, creep, and thermal 
strains, respectively. 

Since the model is of incremental form, the above equation must be 
incremental i zed. Thus, 


dS = s t+At 
ij S ij 


sj. • 


'e - E P - E C 
^ am am am 


-E T V 

an/ 


E - 
an an 


24 


which can be rewritten in the form used in the model as: 


dS 


ij 


a D t+At (dl 
ljnn ^ an 

* dD. . fz z 

ijmn \ am 


d E P - d£ C - dZ T ) 
am am mn j 

E Pt - E Ct - E Tt ) 

mn mn mn / 


( 23 ) 


where the superscripts t denote parameters measured at the start of a lead 
step and the superscript t+At denotes measurement at the end of a load 
step. The terms preceded by d represent the increment in these parameters 
during the load step. 

The flow rule used in the model is the flow rule associated with the 
Von Mises criterion, and is given by: 


dE P . =d*rf- 
13 3S ij 


( 24 ) 


where d X is an unknown scalar to be determined by the consistency 
condition. This equation is commonly called the normality condition, it 
can be seen that the plastic strain increment is assumed to be a 
projection on the outer normal to the yield surface, which is scaled by 
the parameter dX . 


The final expression required to complete the const_tutive law is the 
hardening rule. According to Ziegler's modification of the Prager 
hardening rule, a tensorially correct statement describing translation of 
the yield surface during a load increment is: 


da. . s d (S. . - a. .) 
ij m ij iy 


( 25 ) 
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where the yield surface translation scalar, dy, contains the temperature 
and plastic strain history-dependence of the yield surface translation 
tensor, o^j. The translation scalar is described as: 

(26) 


aT-dS.. 
8S. . ij 

dp = -ii- 


3k 3k . -P 

2k r= dT - 2k — r dc 

37 3 £ P 


(S 


am 


„ ^ -SI- 

OUX'' 3S 
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These components described above are used to derive the constitutive 

law used in the Haisler-Allen model. The resulting stress-strain relation 
is: 


dS. . = C . IdE 
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Also, H' is the slope of the uniaxial stress-strain diagram for the 
temperature at the end of the load step, and Scr/s T is the change in 
uniaxial stress due to a temperature change for the total uniaxial strain 
at the start of the load step. 

It can be seen that the form of the constitutive law is an effective 
modulus C* multiplied by a strain increment added to the change in 
the effective modulus multiplied by the elastic strain at the start of the 
step. The term dP^j is a correction term which arises from estimating 
H' using the uniaxial strain at the start of the step. 

The general sequence followed in the implementation of this model in 
the finite element codes is to solve for displacements from the given 
loading and previous history. These displacements are then used to 
calculate an increment in total strain which is fed to the constitutive 
model along with the temperature increment. A stress increment is then 
calculated initially assuming zero creep and plastic strain increments. 
This stress increment is used with the yield function to check for plastic 
loading. If there is none, the results are correct and the procedure for 
the load step is complete. 

If loading is predicted, the effective modulus ( c i j mn) > the 
effective modulus increment (dlC^j mn ). , and the correction factor 
(dPij), are calculated and used to determine the stress increment. The 
yield surface radius (k) is then updated and the translation of the yield 
surface center (da^j) is calculated to complete the procedure. 

This process is carried out over the entire structure resulting in a 
pseudoforce increment caused by the plastic strain increments. The 
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process is then repeated with the just-calculated pseudoforce increment 
(along with the other loads applied to the structure) until convergence is 
achieved. 


2.1.3 SIMPLE CREEP MODEL 

To make possible the calculation of time-dependent response with the 
two classical plasticity models, a simplified creep model has been 
included with these plasticity implementations. This has been done in 
such a way that the mode'.s can be used separately in creep-only or 
plasticity-only applications, or they may be used simultaneously for 
certain applications (generally in material testing) where both 
time-dependent and independent responses may occur. The implementation 
allows for equal time increments or dynamic time- incrementing as developed 
for use with the Bodner model. 

The simplified creep model is a linear model based on secondary creep 
only and is represented as: 



where Q and r are constants. The integration scheme used in the 
implementation of this model is a trapezoid rule (using the stress value 
at the beginning and end of the time step) to compute creep strain rates 
at those times, so that 


Ac 



(32) 


where : 


•c 

c i 


QC V 


(33) 
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2 . 1. 3 . 1 EXAMPLES WITH PLASTICITY/CREEP MODELS 
For the classical plasticity models, uniaxial stress cases have been run 
to check out the constitutive routines. One set of four plasticity-only 
cases are very useful since they involve both thermal and mechanical 
loadings in various combinations. Note that the total thermomechanical 
load is the same for these four cases, but different residual strains 
result due to the various loading sequences. The following Figures (8 
through 13) show the material data used, the four thermomechanical load 
histories, and the results. 

An example of the creep-model-only is shown in which the total strain 
is imposed and creep strains vary accordingly. This results in a stress 
relaxation as time increases. The creep coefficients used are chosen to 
simplify the calculations as Q ■ 0.44E-8 and r ■ 1. 

Finally, two examples of plasticity and creep combined are shown. The 
plasticity data used is the same as for the plasticity-only examples, and 
the creep coefficients are the same as for the creep-only example. The 
results are encouraging in that the plasticity results are unaffected by 
the creep, which is the assumption of the classical plasticity theory. 

The creep results are improved on the total strain, given that the loading 
is done slowly enough to introduce time-dependent action but quickly 
enough for time- independent action to occur simultaneously. The use of 
this option is determined by the user, although some initial examinations 
of time-step sizing indicate that for slow loadings when creep-only is 
expected, plasticity does not occur (even with the plasticity calculations 
included) given small enough time steps. 
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Figure 11. Strain Controlled Creep. 



<Q>-— — <Q> Elastic 

□ □ Plastic 

A A Creep 


in/ in 










2.1.4 BODNERjS UNIFIED MOD EL 

Bodner f s constitutive model is considered a unified model since 
plastic and creep strains are included in a single inelastic strain 
measure. Thus, the total strain is the sum of elastic, inelastic, and 
thermal strains. The inelastic strain rate tensor is defined in a form 
similar to the Prandtl-Reuss flow rule. It is: 



Do exp 




(34) 


where Do is the limiting strain rate in shear. The material constant n 
controls the strain rate sensitivity and J 2 is the second invariant of 
the deviatoric stress tensor, S i j. The state variable Z is a measure of 
resistance to inelastic flow and is given an initial value of Zo. 

The state variable evolution equation is: 



where the first term defines the hardening and the second term 
characterizes the thermal recovery. The material parameter Z^ is the 
maximum value of Z and Z 2 is the minimum value of Z obtained in thermal 
recovery. The inelastic work rate is W 1 and m, A, and R are material 
constants. 

The Bodner model presented here predicts only isotropic hardening, 
since Z is a scalar quantity. The material parameters are easily 
calculated from standard stress strain curves and creep tests, and the 
numerical implementation is relatively simple. 
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In order to test the Bodner model, the solution, and dynamic time 
incrementing schemes discussed in Section 2.3.3 within the context of a 
finite element code, a two-dimensional model of the benchmark notch 
specimen (References 10 and 16) containing more than 1,000 constant strain 
triangles has been run (Figure 14 ). This code uses Bodner 's constitutive 
model, and the first cycles of Tests 8, 9, and 10 of the benchmark notch 
program were simulated (Figures 15 through 18 ). The results were quite 
sa ^i s fying, both in the performance of the constitutive model and in the 
economy of the solution. On General Electric's Honeywell 6000 computer 
systems, each of these analyses required approximately 3 hours of CPU time 
at a cost of about $500 each. With similar codes on the CRAY-1 and the 
Honeywell 6000, we have experienced speed ratios of about 50 to 1. 
Therefore, on a machine such as the CRAY-1, we would expect such an 
analysis to use about five minutes of CPU time. The three analyses 
required between 205 and 219 time sub increments. 

It can be concluded that this code and constitutive model show promise 
as a design tool. 



Figure 14* Two-Dimensional Model of Benchmark Notch Specimen. 
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Figure 16. Cycle Test 8. 
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Load, lbs 



Figure 17. Cycle Test 9. 
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Load, lbs 



Figure 18 


Cycle Test 10. 


2.2 FORMULATIO N MODELS 

Three formulation models have been developed, a nine-noded degenerated 
shell element, an eight-noded degenerated shell element, and a 20-noded 
isoparametric solid element. The implementation of the theoretical 
formulation of these elements in computer software has been accomplished 
in such a way as to optimize their utility for nonlinear analysis. All 
have been implemented with an out-of-core modified, blocked skyline 
equation solver so that very large, real world problems can be solved. 
■Emphasis has been placed on numerical accuracy, user friendliness, and 
economy of the situation. 

2.2.1 EIGHT-MODE DEGENERAT ED SHELL 
The analysis of shell structures having complex shapes presents an 
intractable analytical problem in the classical theory of shells. In 
addition, if the shell structure is thick and shear deformation is 
significant, the application of classical thin-shell theory becomes 
questionable. The finite element analysis of such structures is a 
feasible alternative for obtaining numerical solutions. In modeling shell 
structures, curved surfaces and faces are often approximated by flat 
elements with straight sides. For accurate representation of the 
structure geometry, many of these elements are necessary. The need for 
elements with curved surfaces is rather obvious. 

Isoparametric mapping provides a logical means for developing such 
elements. The conventional flat-shell elements are often a combination of 
plane stress and bending elements, but the double-curved shell element has 



been derived by modifying the three-dimensional, isoparametric, 
solid-element formulation in a manner consistent with the shell-theory 
assumptions without restricting behavior to purely thin shells. 

The constraint of a straight normal to the midsurface is introduced, 
and the strain energy corresponding to the direct stress normal to the 
middle surface is neglected to conform to the shell theory. However, the 
normal to the middle surface does not remain normal after deformation. 

This permits the element to experience transverse shear deformations 
necessary for thick-shell applications. 

A temperature-dependent property can be accurately determined using 
element nodal temperatures. The element loads consist of thermal loading, 
distributed or uniform pressure loading on element faces, and body forces 
consisting of centrifugal and gravitational loads. A rotated local 
Cartesian coordinate system may also be defined at each node of the 
element. This feature is useful in constraining the rotational degree of 
freedom about the shell normal. 

Consider a 3-D, isoparametric, solid element where linear edges 
connect higher order faces containing an equal number of nodes. Figure 19 
shows a 16-noded, 3-D, isoparametric, solid element where 
eight-noded parabolic faces are connected by linear edges. The external 
faces of these elements are curved, but the sections across the thickness 
are generated by straight lines. For such elements, the pair of Nodes 
it an< * given Cartesian coordinates completely describe the 

element geometry. In 3-D curved shells, this geometry is represented by a 
middle surface containing the element nodes and the normal to these nodes, 
such that V 3i is a unit normal at Node i, and t^ represents the 
element thickness at Node i (Figure 20 and 21) . 
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The unit vectors V^, V 2 j_, and V 3 j_ define an orthogonal local triad 
at Node i. 



1 


“1 
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STRESS AND STRAIN CONSTITUTIVE EQUATION 

Definition of proper stresses and strains consistent with the shell 
theory is necessary for deriving element properties. Consider the shell 
surface defined by 5 ■ constant. At a point on this surface, construct 
orthogonal Cartesian axes x', y', and z' such that x'y' lie in the tangent 
plane and z' is normal to the x'y' plane. The strain components of 
interest can be defined in the local system x'y'z'. 

{e 1 } * 6y'i e x'y' ' e y V » 1 

The strain e 2 ' has been neglected to satisfy shell assumption. The 
stresses {a'} are related to the strains through the elasticity 
matrix: {a'} » [c'l (e'- e o '} where the elasticity matrix [ c ' ] is a 

5x5 symmetric matrix. 

NUMERICAL INTEGRATION 

The calculation of element matrices and the load vectors require 
evaluation of the integrals over the element volume or the faces of the 
element volume. This is accomplished numerically using quadrature. 
Reduced integration is used in 5 and n directions to improve the bending 
behavior. A (2x2x2) integration rule has been chosen for the eight-noded, 
parabolic-shell element. 
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V 

i 

ROTATED LOCAL COORDINATE SYSTEM 

With the basic five degrees of freedom (Uj., V i# W i , Ti , 3 L ) , 
the element does not permit sharp junctions in the structural models. To 
avoid such difficulty, the rotations y.j > ere transformed into 
0 X1 -, 0 ^ © 21 - . However, the element still does not have rotational 
! . stiffness about the normal to the middle surface. Therefore, this 

rotational degree of freedomn should be constrained at the nodes where the 
surface is geometrically smooth before solving for the unknown nodal 
r- displacements. A local Cartesian coordinate system is established at the 

! nodes requiring such constraining in such a way that one of the axes is 

p normal to the middle surface of the shell. The element stiffness, mass, 

end the load vectors are transformed into this local coordinate system, 
j The rotational degree of freedom about the axis normal to the middle 

surface is then deleted before assembly and solution. 

r 

EXAMPLE CASES 

This element has been implemented in a finite element computer 
program, and the program debugged, verified, and validated. Three of the 

r 

test cases are shown in Figures 23, 24, and 25. These cases are very 
severe tests because the geometries investigated are on the lower end of 
plate theory and thus involve both beam and plate characteristics for 
which exact, closed-form solutions do not exist. 

j 

To demonstrate its nonlinear capability, the compact tension specimen 

shown in Figure 26 was modeled (utilizing the line of symmetry) with 30 
eight-noded shell elements and 114 nodes. The results of the analysis 
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using the simplified constitutive model is shown in Figure 27. It compared 
favorably with other nonlinear finite element analyses (CYANIDE, an 
in-house program) and published data (Reference 7) . 



E - 30 x 10 6 psi 
v - 0.3 
a • b * 2 inch 
t ■ 0.2 inch 


A cantilever plate is subjected to area pressure-load/line-distributed 
loads as shown above. The square plate was modeled by 8 elements and 37 
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nodes. The comparison of displacements between beamplate theory and 


finite-element results are listed below: 


Loading Conditions 
. M, P, Q, 

Iin-lbf/jn) nbf/in ) (psii 


Vertical Displacem ents finch) 
Beam Theory 

Plate Finite-Slement Result 

fA. B. and & £ £ 


1000 

0 

0 

(0.1 to 

0.091) 

0.0979 

0.0963 

0.0915 

0 

0 

1000 

(0.1 to 

0.091) 

0.0953 

0.0946 

0 . 0938 

0 

1000 

0 

(0.1333 

to 

0.1284 

0.1274 

0 . 1232 


0.1213) 


t 


1 [ M ♦ ♦ M ♦ f 



ss 


Figure 24. Example Case 2. 


E - 30 X 10 6 psi 
V 0.3 

a * b ■ 2 inch 
t ■ 0.2 inch 
Q - 1000 psi 


A square plate with simply supported (SS) boundary conditions along all 
edges is subjected to the uniform pressure over the entire plate. The 
plate was modeled by 8 elements and 37 nodes. The comparison of the 
maximum vertical displacement at the center of the plate (Location A) was 
made between the finite-element result and Roark's published solution. 

Finite-Element Analysis: <5^ ■ 0.003165 inch 

Roark's Formula: ^ ■ 0.00296 inch 




Figure 25. Example Case 3. 


E • 30 x 10 6 psi 
v » 0.2 
a ■ 2 inch 
b » 1 inch 
c * 0.2 inch 
Q « lO^ psi 


A parallelogram plate (skew slab) with all edges simply supported (SS) is 
under uniform pressure load all over the entire area. The plate was 
modeled by 8 elements and 37 nodes. At the center of the plate (Location 
A) , the maximum displacement was obtained by finite-element analysis and 
from Roark's published formula as shown below: 

Finite Element Analysis: $ A » 0.00365 inch 

Roark's Formula: $ A « 0.00393 inch 


J 
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DUpUce»«nt. U (10“ 3 Inch) 








2-2.2 MINE-NODE DEGENERATED SHELL 
The basis for this element, as it is for most of the others, is the 
principle of virtual work for displacement increments, written as: 


f * /* - - r ~ 4 

J c 6e dv ■ / f «U dv ♦ J. T 6u 


where u,e ,a J and £ are the displacement-, strain-, stress-, 
traction- and body force-rate fields, respectively. The displacement-rate 
field is related to the nodal displacement rates {u} through shape 
functions appropriate for each individual element by: 

u = [K] {u} 

(37) 

The shape functions [N] are functions of the spatial coordinates, and 
are chosen to satisfy admissibility conditions for u as well as particular 
characteristics for each element. The strain field is found by taking- the 
appropriate combination of partial derivatives of u with respect to the 
spatial coordinates, and may be written as: 

* ■ IN’ ] {u} ,, ox 


The stress-strain law is in the form (see Section 2.1) 


o a [E] e - T 
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Substituting (37) -(39) into (36) and taking variations with respec: 
to {u} yields the element equations: 

(K] {u} * {F IE } ♦ {F 8 } ♦ {F 1 } ♦ {F} 


where : 


m -X 


(N'J (E] IN' ]dV 


(41. a) 


{i IE } = J v (N'] T t dV 


(41. b) 


{F®} = f w [H] T £ dV 


(41. c) 


{f} « f sa [N] T T dS 


(41. d) 


and where {F} are the nodal force rates exerted on the element by its 
neighbors. Generally, some of the body force field is due to structural 
vibration. This ultimately yields a mass matrix given by 

IM] * J y [N] T p IK] dV (41. e) 
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where P is the mass density of the element. The objective of any of the 
element packages is to determine five expressions (41) appropriate for a 
given element, and from geometric, material, and load data, determine 
numerical values for the entries in the particular vectors and matrices. 

The nine-node degenerated shell element formulation is a modification 
of the element developed by Chang, et al. (References 3 and 4). The 
Figures illustrate the geometry of the element. The middle surface of the 
shell is generated by mapping a 2x2 square into the surface using: 

• 

r 0 * ( V V V * tN ( $» H)1 U*)» (42) 

where (x Q , y Q , z Q ) represents a point on the midsurface of the shell 
and {{jc}, {y}, { 2 ].} is the collection of node point coordinates, one of the 
given bits of information in the element development. The shape 
functions IN] are the bequadratic shape functions for the 2x2 square. 

Also given at each node point are the shell thicknesses {h} . 

The thickness is assumed to be measured locally normal to the shell 
midsurface. We determine the normal direction by first obtaining e° 
and e°, vectors tangent to the midsurface at a particular point, and 

T\ 

taking their cross product. These vectors are found by: 


H 

I 

i 


n 

i 


'""*T 


! — ) 


H 


H 


•( * if s [||Ju*}. {y}, {*}} 

•Vtf <*>• <*» 


n 
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The normalized cross-product is labeled v 3 , defined in symbols as: 


v - • X $/!.{ x «„l * (» 3> ,, * 3i j. v 3iJ ) ( W , 

The nondimensionalized coordinate in this direction is labeled 4. 
Thus, any point in the shell is represented as: 

r * Cx.y.z) = (N(5,f|)]{{x ♦ -*C*v 3 J, {y ♦ -*?hv 3 J, {z ♦ -*$hv 3>3 }} (45) 


An individual element of the vector {x « 1/2 3hv 3 for 
instance, is x L + 1/2 Shi V 3'l' where x^hi, and v i 3 are 
the x, coordinate, thickness, and x-component of the v 3 vector at node 
i, respectively. The z dependence, then, is assumed to be interpolated 
in the same way as the midsurface coordinates, without 

using the (£,n) explicit function implied by (44). The three 
directional vectors for the mapping (45) are: 


** a ff = * 2 {hV 3,2 }> {hV 3,3 ] 

en = an SS *^* , 2^ {|n] ^ hy 3,l*’ ^ hv 3,2^’ * hv 3,3^ 
= ic = 2 ^ hv 3,l*’ * bv 3,2*’ * bv 3,3^ 
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The Jacobian matrix (J) is formed by: 


[J] 



( 


( 47 ) 


In general, (J) is a function of ( 5 , n , ; ): however, interest in (J) 
will be at the nodal points i (i»l,...,9), where z » 0, and above and 
below the nodal points, where C 8 1, for a total of 27 points. These 
27 points will be the numerical integration points for the element. 

The ( 5 , n ,C) coordinates are convenient to describe the geometry . 

In addition, e z is parallel to v, a t the nodal points. Note that 
the e z and e n are not necessarily perpendicular to each other, (and 
away from the midsurface, probably not perpendicular to eC) , and thus do 
not form a local Cartesian system at any point. However, from a shell 
analysis point of view, it is desirable to distinguish between in-plane 

behavior. In order to accomplish this, the (1, 2, 3) coordinate system 

• 

is introduced where v 3 is the unit vector orthogonal to the 

* * 

midsurface, and and v 2 are defined orthogonal to each other in the 
plane of the midsurface. The v^ vector is defined by normalizing the 
cross product of v 3 and the Cartesian unit vector most perpendicular 
to it; then v 2 ■ v 3 x 

At each node point of the element, there is defined the six degrees 
of freedom (DOF) for that node point. They are the three translations 
U x ,U y , and U z and the three rotations $ x , Q y , and e z . Note 
that these are defined relative to global DOF to local DOF in the (1, 2, 
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3) system. Actually, the local translations will be identical to the 
global ones, and only the rotations will be transformed. As in most 
shell analysis, the DOF 03 will be suppressed. The transformation at 


node point is 

accomplished through the 

matrix 

given as : 


1 

0 

0 

0 

0 

0 



0 

l 

0 

0 

0 

0 


(IJ = 

0 

0 

1 

0 

0 

0 

(48) 
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i 

V 3,3__ 
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The collection of (T i ) 's along a diagonal is denoted as (T) ; thus, 


for the element 
■local 1 


{i loc,1 }= !T]{u« Ub * 1 } 


(49) 
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The (TJ matrix is 54x54 and there are 54 DOF's for this element. 
Another useful matrix is the tensor transformation matrix [Ti] for 
node point i, given as: 


(Ti] = 



<’m> 2 

< v i!/ 

y 3 > 2 


v iN 

V 1 *2 








V 2^2 



" 3 !i> 2 


<* 3 y 


'-all 

V 3^2 



»*!. 

2v * v * 
1.2 2,2 

2v ^ v 
*1.3 2,3 


<2 

" 2*1 


*at, 

’ 3 |i 

2v 1 v 1 
2,2 V 3,2 

2v 1 v 1 
2,3 *3,3 

’aS 


*^1 



’.S 

2v 1 v i 
3,2 V l,2 

2v ^ v * 
3,3 1,3 

t 3 S 

<z 

+ v i!i 

'iiz 


(50) 
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The degrees of freedom are defined at the midsurface of the shel 
To find the displacement field at any point in the shell, use: 

“ a IK (4. n)J [Q]{u local J (51) 

where here, (N) is the biquadratic repeated over three rows, thus 
making [N] a 3x27 matrix, and [Q] is 27x54 made up of nine 3x6 
submatrices (QJ arranged along its diagonal, each [Qj being: 


(52) 



1 

0 

0 

• 2 


0 

IQJ = 

0 

1 

0 

- 1 °*1 ^2 

i '^2 

0 


0 

0 

1 

• i ’at, 

2 ’M 

0 


Note that u is defined in terms of the global directions and the 
zero sixth column of [qj indicates the independence of "u on . 
Furthermore, [Qj clearly represents the usual shell assumption of 
linear variation of the displacement field through the thickness. 


Strains are calculated in Cartesian coordinates first. The relation 
between the differential operators in (x, y, z) and (r, n , - ) are: 



~3/dx 


’a/84 

(£’] = 

a/3y 

* Uf 1 

a/an 


a / K 


a/a? 

m « 


Denoting as [F‘] a 9x3 matrix with 3 [ f ’ ] arranged on the "main 

diagonal," and [D] , a 6x9 Boolean matrix denoting a linear combination 
of partial derivatives of u, the strain displacement relation is written 
as: 

* * (dj an m iqd m {i * loba1 } (54) 

where use has been made of (49) . The combination of matrices on the 
right-hand side of (54) is denoted as (s 8- - 0041 ] . 

The strain relative to (x, y, z) is generally not interesting from a 
physical point of view. It is desirable to obtain the strain in the 
(1, 2, 3) system where it physically has greater meaning. Thus, using 
(50) , 

[3 = ITJ ta (55) 

where the i subscript refers to the i th nodal point, and the j 
subscript refers to the j** 1 surface (middle, upper, or lower) . 
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The constitutive law (39) must be modified to take into account the 
mechanics of materials assumption that a 33 ■ o. With (E) and supplied 
in the (1, 2, 3) system for the 27 integration points, a partitition of 
(39) is made as follows: 


d 1 


E „ E 


t 1 



c 

1 

c3 cc 


c 

* 

c 

0 




*33 ' 

J 

1 

|t 33 

1 


(56) 


where a c are the five non-zero stress-value fields, e the 
corresponding five strain fields, and the dimensions of E cc , E c3 , 

e 3c' e 33' t c' and *33' are 5x5, 5xl ' lx5 ' 5xl ' ixl, 5x1, and 

lxl, respectively. Solving the lower partition of e 33 and using it in 

the upper partition yields the modified relation: 



(57) 


where : 


H-W-MWH 

vvMH'i, 


(58) 


(58b) 
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Furthermore, it is convenient to define [3 lccal i as 
..-■i ijc 


li with the row corresponding to deleted; then: 


0. 

ijc 


, r* local, 

] = U cc l (C ijc 1 


( 59 ) 


The matrix [3. local ] can be defined as [ 3 local ] 

ij ‘ t ijc 1 

augmented by a row of zeros corresponding to t 33 = 0. Finally, the 

stress in global coordinates is expressed as: 


• global a * T *-1 local j* global j, _ ♦ local 
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Element matrix and vector calculations are now ready to begin. For 
instance for the stiffness matrix, the integrand [k^j] at the 
(i,j) th integration point is: 

[k. .] = [t. local ] fa local i (61) 

iJ l *ijc 1 1 ijc J 

The volume integration will be taken numerically by approximating 
the integrand in ( 41 . a) by evaluation of the integrand at the 27 
integration points and then interpolating these using the triquadratic 
shape functions in ( 5 , n , c ) . The shape functions are integrated 
exactly over the nondimensionalized cube: 

ij “ •f'-l J -1 O) d£ dq d£ (62) 
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where [N ±j 1 is the appropriate shape function for the (i,j) th 
integration point, then: 

( 63 ) 

[K) ■ Z«y [Ky] d«t (Jy) 
orr 

The vector {? } is calculated in a similar manner; defining 

{7..} * [ 3 . l0C4l ] T t local (64) 

1 1JC 1 c 

then: 

< fIE > * h.j °ij d “ (J « 1 <S5) 

The mass matrix calculated this way yields what appears to be a 
lumped- typed mass matrix. 

In general, the value of the body force field is either given or can 
be calculated at every integration point. Let {£^}'be the body force 
load at the (i,j) th integration point. If {jj} is defined as 6x1 
subvector corresponding to the i th node ' s contribution {F fi } , 



The particular body force loads' of interest are due to the weight, 
linear acceleration, and centrifugal acceleration. For the weight load, 
the mass density must be given as well as the direction consines for the 
downward direction, then: 
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where g is the acceleration of gravity. For linear acceleration, the 
three components of acceleration are needed; thus; 


n 




-p 


( 68 ) 


For centrifugal loading, the three components of the angular 
velocity vector, plus a point A on the axis of revolution must be 
given. Defining: 


0 


[ 0 ] 




(69) 


then: 


{f ij l 
lI B 'cig 


ij 
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The traction loading for this element is generally divided into tvo 
types: pressure loading on the upper or lower surfaces and line loading 

along the -element edges. Pressure loading performs normal to the 
surface (in the v 3 direction) . The (x, y, z) components of the 
traction vector at integration point i(i - 1,9? top or bottom surface) 
is: 


(, 


i 

3,1 





(71) 


where pi is the magnitude of the pressure at node i. The 6x1 subvector 
of: 


{ F T } Pres, 


pres 


is thus: 

{*I} pres s pi {t i jT lQ ij lT p* 1 * «Jj j l 

where 



Ni (4, n) <*4dn 


(72) 


( 73 ) 
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For an edge load, it is first assumed that the load is applied along 
the line C ■ 0. For any edge, let the s direction be directed along an 
edge and let the n direction be normal to it. If the edge load at node 
Ej (where the three nodal points on the edge are labeled E^, E 2 , 

E 3 ), is given in terms of Cartesian components as W e j , then the 6x1 
equivalent subload vector {il.} is given as: 

Ej 


{$} 


line = y, 


Ej 




lQ Ej,o r 


Ej 


|«^| 
so 1 


(74) 


Note that s corresponds to c or n depending on the edge. Sometimes, 
an edge load will be given in terms of normal and shear components. For 
shear loading in the s direction: 



as so 1 so 1 


(75) 


For shear loading in the normal (?) direction: 



= W 2 / V 
aC 3 


Ej 


(76) 


For normal loading: 


hi * ^ * ? 3 i 


(77) 



Other loads can be modeled as necessary in similar fashion. It 
hould be noted that "thermal loading" is incorporated into the 
constitutive law by the x term. In fact, ^ way be portioned into sum of 
terms. Each of these are then manipulated in an identical fashion. 

2.2.3 2 O-NODED. ISOPARAMETRIC SOLID ELEMENTS 
The isoparametric solid elements permit the modeling of any general 
three-dimensional (3-D) object, since the elements represent a 
discretization of the object into finite elements which are 3-D 
continuous representations. The basic term "isoparametric" means that 
the elements utilize the same interpolating functions (also called 
"shape functions") to interpolate geometry, displacements, strains, and 
temperatures. It is, therefore, important that the user be aware that 
not any displacement, geometry, and temperature field to be analyzed is 
necessarily compatible with a given element mesh. This is particularly 
true where high temperature or strain gradients occur. The following 
sections discuss the basic element formulation assumptions and define 
the node numbering order and face definitions for pressure load input. 

Shape functions are used to describe the variation of some function 
G within an element in terms of the nodal point values. 

a 

G(x ( y f *) ■ £ ^78) 

i«l 
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where 


G(x,y,z) 


G i ( x i' 2 ; 


H; 


n 


the 'value of the function (such as 
displacement, temperature) at any 
point with coordinates (x,y,z) within 
an element 

the value of the function at node 
point i 

the element "shape function" 
the number of nodes describing 
intraelement variation. 


In order to ensure nonotonic convergence to the correct results, 
shape functions must satisfy several requirements. Satisfaction of 
these requirements results in convergence from an upper bound. These 
displacement function requirements are: 

e They must include all possible rigid body displacements 
e They must be able to represent constant strain states 
e They must be differentiable within elements and compatible 
between adjacent elements. 
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While the above conditions prove valuable for establishing upper 
bounds for solutions, they are not essential. Incompatible displacement 
modes are widely and successfully used. Their principal disadvantage is 
that stiffness may no longer be bounded from above. 

Curvilinear coordinates are introduced into the isoparametric 
concept to overcome the difficulty of formulating shape functions in 
global Cartesian coordinates. Also, generality in element geometry 
definition is obtained by this process. 

A local curvilinear coordinate system ( 5 ,n ) which ranges from -1 
to 1 within each element is introduced in which shape functions are 
formulated. Also, a mapping from curvilinear to global coordinates is 
defined. A typical two dimensional element is shown in Figure 28. 

The same polynomial terms used in the Cartesian coordinates are used 
but with the curvilinear coordinates C.n,; replacing x, y, and z to 
generate shape functions. The C,n and 3 coordinates are the same for 
all global element configurations. 

Structural finite element equilibrium equations: 

[Ml {«£> ♦ tele U£ } ♦ [K] ■ (Fg) ♦ { F$) ♦ (Fj) ♦ (Fc> ♦ iFfn.} (79) 

(Hi • / v p(Hl T [Hjdv "Caaaiacaac" uu aacrix 
[C] ■ Damp i of aacrix 
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(b) Global Coordinates 


Figure 28. Typical Two-Dimensional 
Element. 


n 
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n 

■ i 







(&] • /v(B]*(Dl [Bjdv Sciffnasa matrix 

(Fj) ■ / v l8j^{fj)dv Body forcea 
{Fg} ■ /g[U^]^ ( fgJdS Surface tract iooa 

{Fi> • / v (3l*(Dj (cfMv Iaitial aeraias 

( ? HL) “ /yt#l*tDl (* kl) <1v Honlinaar acraiaa 
whara 

{ui)T • [uj vj uj v 2 w 2 ....] 

{u}T • [u »«] 

U> • [H] {«£> 

{«> ■ (3){ui} 

{«) • - [OV'{c TOT >-{c THERM )-(cni.> > 

{fg}T • [f Bx f By f B ,J 

{fsJ T * tfs* f Sy f S«l 
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This element has been formulated with variable temperature general 
orthotropic material properties. During numerical integration for 
stiffness and equivalent nodal forces due to thermals, plasticity, and 
creep, the material properties at each integration point are evaluated 
at the temperature of that integration point. A Gauss integration 
scheme is used, and the user may choose an integration order of 2 , 3, or 
4 points in each direction ( 5 , n , 5 ) . 

The 20-noded solid has three displacement degrees of freedom on each 
of the 20 nodes for a total of 60 element degrees of freedom. Since 
this element utilizes a higher order displacement (and thus higher order 
strain) function, it can be used -to model larger regions with fewer 
elements. The displacement functions yield a displacement distribution 
which is parabolic on an edge (three nodal points per edge) . 

The node numbering sequence for this element is shown in Figure 
29. The user may define the location of Nodes and Ng 
as desired, thus establishing the element pressure coordinate system and 
the resulting face numbering. 




r 

i 




The interpolating functions can be defined starting with the basic 
corner noded shape functions: 

C 1 • (1 ♦ 5 ) (1 ♦«)<!♦ ;)/3 C 9 - (1 ♦ 5 ) (1 ♦„)(!- c )/8 

c - (1 - 0 (1 ♦ n) <1 ♦ «)/* G U • (l - 0 (1 ♦ „) (1 - ;)/8 

65 ■ <1 - 0 (1 - n) (1 ♦ «)/a c 13 ■ (1 - 5 ) (1 - „) (1 - ;)/8 

c 7 ■ (1 ♦ c ) (1 - n ) 0 ♦ ;)/8 g 15 • (1 ♦ 0 (l - „) (l - ;)/8 

The midside node shape functions are: 

H 2 - (1 - c 2 ) (1 ♦ n) (1 ♦ ;>/* 

H4 ■ (1 - {) (1 - n 2 ) (1 * 

• (1 - c 2 ) (1 - n) <1 ♦ c)/* 

Hg ■ (l ♦ 5) (1 - n 2 ) (1 ♦ ;)/^ 

H 10 - (l - c 2 > (1 ♦ n ) a - c >/* 

®12 • (1 “ {) (1 • n 2 5 (1 ~ ;)/4 

h 14 • (1 - C 2 ) (1 - n) (1 - ;)/ 4 

>16 •<!♦«) (1 * n 2 ) (1 - c)/4 

>17 ■(!♦?) (1 ♦ n> (1 - ; 2 )/4 

>18 ■ (1 “ ?) (I ♦ n) (1 “ c 2 ) /4 

>19 - <1 - O (1 - n) (1 - ; 2 )/4 

>20 * (1 * C) (1 “ n) (1 “ c 2 )/ 4 


! 
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The modified corner node shape functions are: 

H 1 "Ci - (H 2 * Hg ♦ H 17 )/2 

h 3 " C3 “ (H 2 ♦ H 4 ♦ H la )/2 

h 5 *63 - (84 * H 6 ♦ H l9 )/2 

h 7 -67 - (Hj ♦ H 8 ♦ H 20 )/2 

«9 * g 9 “ (Hio ♦ H I6 ♦ H i7 )/2 

H U * C 11 * (H^o ♦ Hj 2 ♦ Hij)/2 

h 13 " C 13 - (H 12 ♦ H U ♦ H I9 )/2 

H 15 * C i7 - (H U ♦ H 16 ♦ H 20 )/2 


This element requires no incompatible modes to simulate bending 
properly, since the displacement functions are complete quadratics. 
Thus the displacements (and strains and temperatures) are computed as: 
20 

u - r 
i • 1 


20 

v - t E i V i 
i ■ 1 

20 

w ■ 1 

i ■ 1 

Given the coordinate system ($,niC ) as previously established, we 
can also now define the face numbering conventions and order of nodes on 
a face. These definitions are needed to established conventions for 
inputting pressure levels on the element and numbering of faces when 
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displaying surface stresses on the faces. These conventions are 
summarized below. 


:e Ho. 

Location 


Nodes end 

Node 

Order on 

Face 


i 

C 

- *i 

HI 

H8 

H7 

H20 

HI 5 

HI 6 

H9 

HI 7 

2 

n 

- .I 

HI 

HI 7 

H9 

H10 

HU 

HI 8 

H3 

H2 

3 

c 

- .i 

HI 

H2 

H3 

84 

84 

H6 

H7 

H8 

4 

c 

- -i 

HI 3 

H19 

K5 

84 

H3 

HI 8 

HU 

H12 

5 

n 

- -i 

HI 3 

HI 4 

HI 5 

820 

H7 

H6 

H5 

HI 9 

6 

c 

■ -i 

HI 3 

H12 

HI 1 

H10 

H9 

H16 

H15 

H14 


This element has been formulated with variable temperature and 
general orthotropic material properties. During numerical integration 
for stiffness and equivalent nodal forces due to thermals, plasticity, 
and creep, the .material properties, at each integration point are 
evaluated at the temperature of that integration point. A Gauss 
integration scheme is used, and the user may choose an integration order 
of 2, 3, or 4 points in each direction (C,n,; ). 

2.3 COMPUTER PROGRAMS 

A separate computer program has been developed for each combination 
of constitutive model-formulation model. Each program provides a 
functional, standalone capability for performing cyclic nonlinear 
structural analysis. In addition, the analysis capabilities 
incorporated into each program can be abstracted in subroutine form for 
incorporation into other codes or to form new combinations. These 
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programs will provide the structural analyst with a matrix of 
capabilities involving the constitutive model-formulation models from 
which he will be able to select the combination that satisfies his 
particular needs. 

The total amount of software developed for this contract is large. 

It will be presented in separate manuals as required. To illustrate the 
program architecture, Appendix B is a Fortran source listing for a main 
calling program. Appendix C contains a listing of the input for the 
various programs. Appendix D contains the listings of the DATA DECK 
GENERATORS programs. These programs can be used to interactively 
generate the input for the structural codes. 

2.3.1 PROGRAM ARCHITECTURE 

The program architecture employs state-of-the-art techniques to 
maximize efficiency, utility, and portability. Among these features are 
the following: 

(i) User Friendly I/O 

e Free format data input 

e Global, local coordinate system, (Cartesian, Cylindrical, 
Spherical) 

• Automatic generation of nodal and elemental attributes 

• User-controlled optional print out 

Nodal Displacements 

Nodal Forces 

Element Forces 

Element Stresses and Strains 
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(ii) Programming Efficiency 

• Dynamic core allocation 
e Optimization of file/core utilization 
e Blocked column skyline equation solver 
(iii) Accurate and Economical Solution Techniques 
e Right-hand side pseudoforce technique 
e Accelerators for the iteration scheme 

e Convergence criteria based on both the local plastic strain 
and the global displacements. 

Figure 30 is a generic flow chart for these programs. 


2.3.2 LINEAR VARIATION OF LOADS 
The ability to model piecewise linear load histories was also 
included in the finite element code. This capability is particularly 
useful when modeling stress strain tests or fatigue loops, and also for 
certain analysis applications. Since the inelastic strain rate could be 
expected to change dramatically during a linear load history, it is 
important to include a dynamic time-incrementing procedure. The term 
"load case" is used to denote a time period for which the initial and 
boundary conditions are defined and vary linearly between the end 
points. 

In order to incorporate linear load histories into this scheme, the 
total displacement vector is decomposed into elastic and inelastic 
components as: 


{ d T ] = {d E i ♦ a 1 ). 


(.80) 
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The vector {d E } is the displacement due to applied thermomechanical 
forces, and fd*} is the displacement due to the inelastic pseudoforces. 
These displacementt vectors can be calculated using: 

{d E } = [Kf l {F}, ( 81 ) 


and: 

{d 1 } = [Kf 1 {F I} - 


( 82 ) 


The elastic displacement are obtained for the initial and final 
thermomechanical loads in the load case. The elastic displacements at 
any time in the load case are given by: 

W ■ W. *(^) ({<*}< - W.) (83) 

r p 

The vectors {d } q and {d are the elastic displacements due to 
the initial and final applied thermomechanical forces. The current time 
in the load case is t, and t Q and t £ are the initial and final times 
in the load case. 

The displacements due to the inelastic strains at any time during 
the load case are given by: 

(W = {d 1 } 0 ♦ 2 {dd 1 } (84) 

Time 

Increments 
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The vector U r } o is the vector of displacements due to inelastic 
strains at the beginning of the load case, and {ad 1 } is a displacement 
increment due to the inelastic strains during a time step. 

The increment in displacements {dd*} due to the change in inelastic 
strains {dc 1 } during a time step is computed using: 

M - ki-'W 1 } (85) 

The inelastic pseudoforce increment {dF 1 } is calculated from: 

{ar 1 } = i /[ b] t [e] (86) 

Elements 


where {dc 1 } is the change in inelastic strain during the time 
increment . 

At the beginning of a load case, the initial and final elastic 
displacements are computed using equation (81) , and the displacements 

due to prior inelastic strains are computed using equation (82) . The 
total strains at the beginning of the load case are recovered for each 
integration point and the elastic strains are computed from: 

{<*} . {e T } - {«“*} - {e 1 } (87) 


where U } are the total strains and {c 01 ^} are the thermal 
strains. 
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Using the current values of the state variables and the stress state 
variable evolution rates. Before entering the time loop, an initial 
time increment is computed and the inelastic strain increments are 
estimated using a forward Euler integration formula. From the estimated 
inelastic strain increments, an initial estimate is made for the 
inelastic pseudoforce increment using Equation (86). 

The usual technique employed with the initial strain method is to 
assume that the incremental inelastic force {AF 1 }, the corresponding 
displacements arid the inelastic strain increments {Ac 1 }, are 

all zero on the first iteration of a time step. The stability of the 
method can be improved considerably when a forward Euler integration of 
the inelastic strain rates is used to make an estimate of {Ac 1 } , 

{AT 1 } and {Ad 1 } on the first iteration. This method results in an 
initial estimate which is much closer to the solution. In sample cases, 
the overall number of iterations was reduced by more than one half. 

The procedure during a time increment is to estimate the solution on 
the first iteration using a forward Euler scheme as outlined above. 

Then displacements, strains, stresses, inelastic strain rates, and state 
variable evolution rates are computed at the end of the time 
increment. The inelastic strains and state variables are integrated 
over the time increment and an improved inelastic force is computed. 

The procedure is repeated until convergence is achieved at the end of 
the time increment. Figure 31 summarizes the logic. 



2.3.3 ~ DYNAMIC TIME INCREMENTING 
In a computer code that allows a linear variation of loads, a 
dynamic time incrementing scheme is very desirable since large 
excursions in stress and inelastic strain rate are to be expected. The 
procedure used to compute the time increments requires a certain amount 
of initial experimentation to determine appropriate time step control 
parameters. However, once this has been done, the procedure works quite 
well and is a tremendous improvement in economy over a constant time 
increment . 

Three separate time step control criteria are used. These are the 
maximum stress increment, maximum inelastic strain increment, and 
maximum rate of change of the inelastic strain rate. The minimum time 
step calculated from the three criteria is the value actually used. 

Since the calculations are based on values readily available from the 
previous time step, little computational effort is required. 

2.3.4 STRESS INCREMENT CRITERION 
A maximum stress increment criterion is used to control the time 
increment during primarily elastic excursions. This criterion is 
necessary to prevent overshoot of the point where significant inelastic 
strain rates begin. The calculation for the time increment is given by: 


bo 

INC 
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Figure 31* Flow Chart of Finite Element Solution Procedure. 
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where ht R-1 is the previous time increment, ( ao x ) max the 
maximum change in effect stress for all integration points during the 
previous time increment, and is the maximum desired stress 

increment. The value for hff INC is program input and will vary 
somewhat, depending on material constants. Typical values are about 15 
MPa. 


2.3.5 INELASTIC STRAIN INCR EMENT CRITERIO N 
The inelastic strain increment criterion controls the time step when 
the stress and inelastic strain rates are not changing significantly. 
This is given by: 


At k = at 


k-l 



I 

INC 


K-l 


oaxj 


(89) 


The maximum change in effective inelastic strain for all integration 
points during the previous time increment is (hr max ) . 

^INC is the maximum desired inelastic strain increment. The value 
for is program input, and typical valus are about 0.000100. 

2.3.6 RATE OF CHANGE OF THE INELASTIC STRAIN RATE CRITERION 
The E 1 criterion controls the time increment when the inelastic strain 
rate is changing rapidly, such as in the "knee" of a stress strain 



curve. The quantity c 1 is a measure of how close the initial 
forward Euler estimation is to the final converged solution. The 
backward- difference formula: 


£ 


I 


i 




/At, 


k-1 


(90) 


n 

: \ 


1 

i 


n 

I 


n 


is used to estimate e 1 ^. The maximum value of c 1 for all 
integration points, g 1 max is used to estimate the next time step 
using: 


At X = 


2At 


k-1 




e 


(91) 



\ 



\ 


The parameter e is the maximum desired percent error by which the 
initial forward Euler estimation is in error. The value for e is 
program input and typical values are about 0.01. Equation (91) is 
derived simply from taking the difference between a Euler integration 
scheme and the more accurate second order Adams-Moulton method. 


H 


-1 


n 

, f 
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TITLE* 


n A T P =n , 5 /H' 3 /Dc: 


TTMC 


i 3 03 


PAGE = 2 




200090C 

200100 SUBROUTINE MAINPG(LTOT) 

2001 IOC 

200120C EQUATION SOLVER BY PACKED SKYLINE , COLUMN SOLUTION SCHEME 

2001 30C STORAGE DYNAMIC ALLOCATION IS USED 

200140C 

200150CALC0M 

200160C 

200170 DIMENSION I A ( 1 ) 

200180 COMMON A (3000) 

200190 EQUIVALENCE (A,IA) 

200200C 

200210 EQUIVALENCE (W1, RPM 123(1)), (W2.R PM 123(2)), (W3,RPM 123(3)) 

200220C 

200230 CHARACTER J0BID*4(18) 

200240 CHARACTER JUNK*4 
2002 50C 

200260 DIMENSION BUFF(400) 

200270 CHARACTER IFILES*4(7) 

200280 DATA I F I LES( 7 )/ * ; '/ 

200290C 

200300C 

200310 DATA MAXN/35/ 

200320C 

200330 9000 FORMAT ( V ) 

200340 9010 FORMAT ( 6A4 ) 

200350 9100 FORMAT ( 1P6E12.3) 

200360 9200 FORMAT (19A4) 

200370 9300 FORMAT (// 1 X . 70( 1H# )//1X . 18A4// 1 X , 70( 1H# )/ ) 

200380 9400 FORMAT ( A 1 ) 

200390C 

200400 9510 FORMAT ( / ' NDIM = ' , 14/ ' I ANAL Y = ' , I 4/ ' I PLANE= ' , 14/ ' INC0MI=',I4/ 
2004105 ' TREF *',F10.3/' ISOTHE* ' . 14/ ' ITHERM=',I4/ 

2004205 ' NMAT ='.14/' NMT =',I4/' NMTCRP= ' , 14/ ' NORTHO= ' , I 4/ 

2004305 ' NSKEW =',I4/' NCONST= ' , 14// 

2004405 ' KF I XT Y = ',14/' KPDIS ='.I4/ / KSKWBC= ' , 14/ ' KNTHK ='.14/ 

2004505 ' KLOCAL= ' . 14/' KPRCNT= ' , 14// ' MIDXY2=',I4/ 

2004605 ' KXYZCD= ' , 14/ ' KGAUSP='.I4/ 

2004705 ' W1 = ' , F 10 . 3/ ' W2 =',F10.3/' W3 =',F10.3/ 

2004805 ' G 1 = ' . F 10 . 3/ ' G2 ='.F10.3/' G3 =',F10.3/ 

2004905 ' XBARAX= ' , F 10. 3/ ' XBARAY= ' , F 10 . 3/ ' XBARAZ= ' , F 10. 3/ ) 

200500 9520 FORMAT ( / ' NELT YP= ' , 14/ ' MAXEL =',I4/' NUMEL=',I4/ 

2005105 ' MNODE =',I4/' MIDOFN= ' , 14/ ' MSTRAN=',I4/ 

2005205 ' MAXNOD= ' , 14/' MNUMNP= ' , 14/ ' NUMNPS= ' , 14/ ' MAXEQ =',I8/) 

200530C 

200540 9530 FORMAT ( / ' KF I XT Y= ' , 14/ ' NPDIS =',I4/' NSKWBC=',I4/ 

2005505 ' KNTHK =',I4/' NL0CAL= ' , 14/ ' NPRCNT = ' , 14/ ) 

200560C 

200570 9540 FORMAT ( / ' NBLOCK= ' , 14/ ' NEQ =',I8/' NEQBND= ' , 18/ 

2005805 ' NT0TP1 = ' , 18/) 

200590C 

200600 MTOT = LTOT 
2006 IOC 

2006 2 OC SEQUENCIAL FILES IOVSTF AND NSTF 

200630C IOVSTF AND NSTF ARE FILE UNITS FOR OVER-ALL STIFFNESS 

200640C BEFORE AND AFTER BOUNDARY CONDITION MODIFICATION 

200650C 

200660C RANDOM-ACCESS FILE NRED 
200670C 

200680C NRED IS FILE LOGIC UNIT FOR DECOMPOSED STIFFNESS MATRIX 

200690C 

200700 CALL CREATE ( I FI LE ,5000,0.IER) 

200701 CALL CREATE (23, 100,0, I SER) 

200710 CALL CRE ATE ( I OVSTF , 5000, 0, IER ) 

200720 CALL CREATE ( NSTF ,5000,0. IER) 

200730 CALL CREATE ( NRED ,5000,1, IER) 

200740 CALL CRE ATE ( I FCDNN , 5000, 1 , I ER ) 
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200750 CALL CREATE ( I FLK ,5000,1,IER) 

200760 CALL CREATE ( I FSTI F , 5000, 1 f I ER ) 

200770 CALL CRE ATE ( I FEF ,5000,1,IER) 

200780C 

200790 PRINT , ' I NPUT DATA FILE' 

200800 READ 9010. ( IFILES( I ) , 1= 1 , 6) 

200810 CALL ATT ACH( IFILE, I F I LES , 3,0, I ST AT, BUFF ) 

200820C 

200830 LAXDOF® 10000 
200840C 

200850 METHOD® 1 
200860 EPSTOL® 1 . OE-5 
200870 T0LEPS=0 . 01 
200880 TOLDI S=2 . 5E-5 
200890 BET AA =0.75 
200900 BETAB =0.25 
2009 IOC 

200920 READ( IFILE, 9000 , ERR® 1 ) LINE 

200921 REWIND IFILE 

200922 READ( IFILE, 9200) JUNK, JOBID 

200923 L I NEBS= 1 

200924 GO TO 2 

200925 1 CONTINUE 

200926 REWIND IFILE 

200927 READ( IFILE , 9200) JOBID 

200928 L INEBS=0 

200929 2 CONTINUE 

200930 LINETP=LINEBS+ 1 

200931 LINTP1=LINEBS+2 

200932 LINTP2=LINEBS+3 

200933 LINTP3=LINEBS+4 

200934 LINTP4=LINEBS+5 

200935 LlNTP5=LINEBS+6 

200936 LINTP6=LINEBS+7 

200937 LINTP7=LINEBS+8 

200938 LINTP8=LINEBS+9 

200939 LINTP9=LINEBS+10 

200940 PRINT 9300, JOBID 

200950C RE AD (IFILE, 9000) LINE , NDIM , I ANALY , I PLANE , INCOMI , 

200960C& TREF , I SOTHE , ITHERM , NMAT , NMT , NPS , NMTCRP , NORTHO , NSKEW , NCONST , 
200970C& KFIXTY , KPD I S , KSKWBC , KNTHK , KLOCAL , KPRCNT , MIDXYZ , KXYZCD , KGAUSP , 
200980C8 , W 1 , W2 , W3 , G 1 , G2 , G3 , XBARA 
200990C 

201000 CALL READZR(MAXN, NW0RD5 , A , I ERR ) 

201010C 

201020 NDIM =A ( LI NEBS+ 1 ) 

201030 lANALY=A(LINEBS+2) 

201040 IPLANE=A( L1NEBS+3) 

201050 I NCOMI =A(LINEBS+4 ) 

201060 TREF =A(LINEBS+5) 

201070 IS0THE=A(LINEBS+6) 

201080 ITHERM=A< LINEBS+7) 

201090C 

201100 NMAT =A ( LI NEBS+8 ) 

201110 NMT = A ( LI NEBS+9 ) 

201120 NPS =A( LINEBS+ 10) 

201130 NMTCRP® A ( L I NEBS+ 1 1 ) 

201140 NORTHO® A ( LINEBS+12) 

201150 NSKEW =A( LINEBS+ 1 3 ) 

201160 NCONST = A(LINEBS*M4) 

201 170C 

201180 KFIXTV®A( LINEES-M5) 

201190 KPD I S = A ( L I NEBS+ 16) 

201200 KSKWBC=A(LINEBS+17) 

201210 KNTHK =A( LINEBS+ 18 ) 

201220 KL0CAL=A(LINEBS+19) 

201230 KPRCNT=A(LINEBS+20) 
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201240 MIDXYZ=A ( L I NEBS+2 1 ) 

201250 KXYZCD*A (LINEBS+22) 

201260 KGAUSP=A ( LINEBS+23) 

201270C 

201280 W1 *A(LINEBS+24) 

201290 W2 *A ( L I NEBS+25 ) 

201300 W3 = A ( L I NEBS+26 ) 

201310 G 1 =A ( L I NEBS+27 ) 

201320 G2 =A ( LINEBS+28 ) 

201330 G3 = A ( LINEBS+29 ) 

201340 XBARA ( 1 ) *A ( L I NEBS+30 ) 

201350 XBARA ( 2 ) =A ( LINEBS+3 1 ) 

201360 XBARA ( 3 ) =A ( LI NEBS+32 ) 

201370C 

201380 NGAUSS=2 
201390 NGAUST*0 

201400 CALL GAUSS0( NGAUSS , GAUSS, GAUSSW) 

201410 IF(NDIM. EQ.2) GO TO 4 
201420 NGAUST= NGAUSS 

201430 CALL GAUSSO( NGAUST , GAUST , GAUSTW ) 

201440C 

201450 4 CONTINUE 
201460C 

201470 I F ( KGAUSP . EO . 0 ) GO TO 5 
201480 MGAUSP*NGAUSS 
201490 MGAUTP* NGAUST 

201500 NPL0CA=MGAUSP**2*MAX0( 1 , MGAUTP ) 

201510 GO TO 6 
201520 5 CONTINUE 
201530 MGAUSP= 1 
201540 MGAUTP* 1 
201550 NPLOCA* 1 
201560 6 CONTINUE 
201570C 

201580 PRINT 9510.NDIM, I ANALY , I PLANE , INCOMI , TREF , ISOTHE, I THERM, 

20 1 590& NMAT , NMT , NMTCRP , NORTHO , NSKEW , NCONST , 

201 6006- KFIXTY , KPD I S , KSKWBC , KNTHK , KLOCAL , KPRCNT , MIDXYZ , KXYZCD , 

2016 10& KGAUSP, W1 ,W2,W3,G1 ,G2,G3, XBARA 
201620C 

201630 CALL COUNT 
201640C 

201650 PRINT 9520 , NELTYP , MAXEL , NUMEL , MNODE , MIDOFN , MSTRAN , 

20 1 660& MAXNOD f MNUMNP , NUMNPS , MAXEO 
201670C 

201680 REWIND IFILE 
201690C 

201700 READ( IFILE ,9400) JUNK 
201710 READ( IFILE, 9400 ) JUNK 
201720C 

201730 CALL RANSI Z ( I FCONN , MNODE , 1 ) 

201740 CALL RANSIZ( I FLK ,MIDOFE,l) 

201750 CALL RANSIZ ( I FEE .MID0FE.1) 

201760 CALL RANS I Z ( I F ST I F , MNES ,1) 

201770C 

201780 I F ( NDIM . LE . 2 ) NGAUST=0 
201790C 

201800 I F ( NCONST . EQ . 0 .OR. KFCONS.LE.O) GO TO 10 
201810 CALL CREATE ( I FNCON , 5000 , 1 , I ER ) 

201820 CALL CRE ATE ( I FFCON , 5000 , 1 , I ER ) 

201830 CALL CRE ATE ( I F ECON , 5000 , 1 , I ER ) 

201840C 

201850 CALL RANSI Z ( I FNCON , LCONLK , 1 ) 

201860 CALL RANSI Z ( I FFCON . LCONLK , 1 ) 

201870 CALL RANSI Z ( I FECON , LCONES . 1 ) 

201880C 

201890 GO TO 15 
201900 10 CONTINUE 
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201910 I FNCON*- I FNCON 
201920 I FFCON*- I FFCON 
201930 I FECON*-I FECON 
201940 15 CONTINUE 
201950C 

201960 NDIMP 1 =NDIM+ 1 
201970 NDIMSQ=NDIM*NDIM 
201980C 

201990 W1SQ=W1**2 
202000 W2SQ=W2**2 
202010 W3SQ=W3**2 
202020 WMATRX( 1 , 1 )=W2SQ+W3SQ 
202030 WMATRX (2,2) *W1 SQ+W3SQ 
202040 WMATRX ( 3 , 3 ) =W 1 SQ+W2SQ 
202050 WMATRX ( 1 ,2)=-W1*W2 
202060 WMATRX ( 1 t 3)=-W1*W3 
202070 WMATRX (2,3)=-W2*W3 
202080 WMATRX (2,1) = WMATRX (1,2) 

202090 WMATRX (3,1) = WMATRX (1,3) 

202100 WMATRX (3,2 )= WMATRX (2,3) 

2021 IOC 

202120 GMATRX ( 1 ) = G1 
202130 GMATRX ( 2 ) = G2 
202140 GMATRX ( 3 ) =G3 
202150C 

202160 CALL ADRESS 
202170C 

202180 PRINT 9530 , KF I XTV , NPDI S , NSKWBC , KNTHK , NLOCAL , NPRCNT 
202190C 

202200 PRINT 9540 , NBLOCK , NEQ , NEOBND , NTOTP 1 

2022 IOC 

202220 I F ( I ANALY . EQ . 1 ) GO TO 60 
202230C 

202240 CALL CREATE ( I FSS , 5000, 1 , I ER ) 

202250 CALL RANSIZ(IFSS ,LRSS3D,1) 

202260C 

202270C INITIALIZATION FOR IN-ELASTIC IFSS DATA 

202280C 

202281 TOTCRT = 0 . 

202290 DO 40 I=1,LRSS3D 
202300 A ( I ) =0 . 0 
202310 40 CONTINUE 
202311 A ( 1 12) =TREF 

202320C 

202330 I PRECD=0 
202340 DO 50 I E L= 1 , NUMEL 
202350 DO 50 I=1,NPL0CA 
202360 I PR E CD* IPRECD+ 1 

202370 CALL GENIOF ( 2 , I FSS . IPRECD . LRSS3D , LRSS3D , A , A ( IPSS ) ) 

202380 50 CONTINUE 
202390C 

202400 60 CONTINUE 
2024 IOC 

202420 LCASt = 0 
202430 I F I RST =0 
202440C 

202450C INITIALIZE TOTAL MECHANICAL LOADS AT TIME (T) 

202460C TOTAL DISPLACEMENTS AT TIME (T) 

202470C INCREMENTAL LOAD FOR THE CURRENT INCREMENT TIME 

202480C 

202490 DO 100 1=1 ,NEG 
202500 A( IPFM1B+I )=0.0 
202510 A ( I PDI SB+I ) =0 . 0 
202520 A ( I PDFB +1 )=0.0 
202530 100 CONTINUE 
202540C 

202550C INITIALIZE NODAL TEMPERATURES AND TEMPERATURE GRADIENTS 
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202560C 

202570 I POUTB* I P0UTN- 1 
202580 00 150 I*1,NUMNPS 
202590 I A ( I POUTB+I ) *0 
202600 A ( IPTN1B+I ) *TREF 
202610 A( IPTN2B+I )=TREF 
202620 IF(KNTHK.EQ.O) GO TO 150 
202630 A ( I PDT 1 B+ I ) *0 . O 
202640 A( IPDT2B+I )*0. O 
202650 150 CONTINUE 
202660C 

202670C INITIALIZE ELEMENT MATERIAL # 

202680C 

202690 DO 160 I=1,NUMEL 
202700 IEL*I 

202710 CALL F LAG( 2 , 1 . I DMATN , 1 , I A ( I PFLGB+ 1 EL ) ) 

202720 160 CONTINUE 
20273OC 

202740C TO INITIALIZE ELEMENT STIFFNESS RE-COMPUTATION CODE 

202750C ELEMENT AVERAGE TEMPERATURE 

202760C 

202770 I STIFF* 1 
202780C 

202790 IPTEMB=IPTEME-1 
202800C 

202810 DO 170 JEL= 1 , NUMEL 
202820 I EL* JEL 

202830 CALL FLAG( 2 , 1 , IDSTIF, 1 , IA( IPFLGB+IEL ) ) 

202840 A( IPTEMB+IEL)=TREF 
202850 170 CONTINUE 
202860C 

202870 200 CONTINUE 
202880C 

202890 LCASE=LCASE+ 1 
202900C 

202901 PRINT. "LOAD CASE".LCASE 

202910 DO 300 1=1 .NEQ 

202920 A ( I PFM2B+ I ) =0 . O 

202930 300 CONTINUE 

202940C 

202950C TO INPUT LOAD PARAMETERS AND LOADS 

202960C 

202970 CALL LOADS ( I A ( I PF LAG ) , I A ( I PI CON ) . I A ( I PUCON ) , I A ( I PE I NT ) . I A ( I PNINT ) , 

202980& A ( I PDOFB ) ,A( I PLM ) f A(IPFM1 ),A(IPTN1 ),A(IPDTN1), 

202990& A ( IPFM2 ),A(IPTN2 ) . A ( I PDTN2 ) , TOTCRT , DELT . I ER ) 

203000C 

203001 PRINT. "DONE OF LOADS" 

203010 IF(IER.NE.O) GO TO 990 
203020C 

203030C PRINT . ' FM 1 ' , ( A ( I ) . I=IPFM1 . ( IPFM1+NEQM 1 ) ) 

203040C PRINT. ' FM2 * . ( A(I ) . I=IPFM2, ( I PFM2+NEQM 1 ) ) 

203050C 

203060 DO 900 INC=1.L0DINC 
203070C 

20307 1 TOTCRT=TOTCRT+DELT 

203080C PRINT, ' INC . I ST I FF . NTEMCH ' , INC . I ST I FF , NTEMCH 
203090C 

203100C CURRENT TOTAL MECHANICAL LOADS 

2031 IOC 

203120 DO 360 1*1 ,NEQ 

203130 A ( IPFOB+I ) *A( I PFM1B+I )+A ( I PFM2B+I ) 

203140 360 CONTINUE 

203150C PRINT. ' FO ' . ( A ( I ) , I*IPFO, ( IPF0+NEQM1 ) ) 

203160C PRINT. 'NTEMCH" .NTEMCH 
203170C 

203180 I F ( NTEMCH . EQ . 0 ) GO TO 400 
203190C 
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TITLE- DATE =03/13/85 TIME « 13.03 

203200C CURRENT NODAL TEMPERATURES AND THE GRADIENTS 

2032 IOC 

203220 DO 370 1=1, NUMNPS 

203230 A ( I PTNB+ I )=A( IPTN1 B+I )+A(IPTN2B+I ) 

203240 I F (KNTHK . EO . 0) GO TO 370 

203250 A( IPDTNB+I )=A( IPDT1B+I )+A( IPDT2B+I ) 

203260 370 CONTINUE 
203270C 

203280C PRINT, 'TN' , (A( I ) , I=IPTN, ( IPTN+NUMNPS- 1 ) ) 

203290C 

203300C TO COMPUTE ELEMENT AVERAGE TEMPERATURE 

2033 IOC 

203320 CALL ELTEMP( I A ( IPFLAG) , I A( IPICON ) , I A ( IP JCON ) , I A ( I PNINT ) , 

203330& A(IPTN),A(IPTEME)) 

203340C 

203350C PRINT, 'TEME' . ( A( I ) , I -IPTEME , ( IPTEME+NUMEL- 1 ) ) 

203360C 

203370 400 CONTINUE 
203380C 

203390 I F ( I STIFF . EO . 0 ) GO TO 500 
203400C 

203410 CALL ELSTI F ( A ( IPFLAG) ) 

203420C PRINT, ' DONE OF ELSTI F ' 

203430C 

203440C TO ASSEMBLE OVER-ALL STIFFNESS AND DETERMINE MAX. DIAGONAL VALUE 

203450C 

203460C 

2034 70C CALL ADDSTF ( NEOBND , A( IPFLAG ) , LK , ES , MAXMIN, MAXA , NCOLBV , A ARRAY ) 

203480 CALL ADDSTF ( NEOBND , A ( I PFLAG ) . A ( I PDOFE ) , A ( I PSTI F ) , A ( I PMXMN ) , A ( IPMAXA ) , 
203490& A ( I PCLBV ) , A ( IPMATA ) ) 

203500C PRINT, ' DONE OF ADDSTF ' 

2035 IOC 

203520C TO MODIFY OVER-ALL STIFFNESS DUE TO SPECIFIE DISP. B.C. 

203530C 

203540C CALL MODSTF ( NEOBND , MAXA , NEODI S , NCOLBV , AARRAY ) 

203550 CALL MODSTF ( NEOBND , A ( IPMATA ) ) 

203560C PRINT, ' DONE OF MODSTF * 

203570C 

203580C TO DECOMPOSE THE MODIFIED OVER-ALL STIFFNESS 

203590C 

203600 CALL DECOMP 

2036 IOC 

203620 500 CONTINUE 
203630C 

203640C PRINT , 'FO' , ( A( I ) , I=IPF0, ( I PFO+NEQM 1 ) ) 

203650C PRINT . 'FM1 ' , ( A( I ) , I=IPFM1 , ( I PFM 1 +NEQM 1 ) ) 

203660C TO COMPUTE INCREMENTAL PSEUDO THERMO-MECHANICAL LOADS 

203670C 

203680C A ( I PDF ) IS A WORKING ARRAY 

203690C 

203700 CALL INCLOD( A ( IPFO ) . A( IPFM1 ) , A ( I PDF ) , A ( IPDFO) ) 

2037 IOC 

203720C PRINT, 'INCREMENTAL ELASTIC THERMO-MECHANICAL LOADS * 

203730C PRINT, ( A ( I ) , I = I PDFO , ( I PDFO+NEQM 1 )) 

203740C 

203750 IF ( I ANALY .NE . 1 ) GO TO 600 
203760C 

203770C ELASTIC ANALYSIS 

203780C 

203790 CALL ELASTY 
203800 GO TO 800 

2038 IOC 

203820C IN-ELASTIC ANALYSIS 

203830C 

203840 600 CONTINUE 
203850C 

203860 CALL PLASTY ( MSTRAN , NDI M , NPS , NMT , I A ( I PI CON ) . I A ( I PUCON ) , 
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TITLE” 


DATE ”03/13/85 


TIME = 13.03 


PAGE” 8 


2038708 IA( IPDOFE ) , I A ( I PNINT ) , A ( IPDF ) , A ( IPDFO) . A( IPFP ) , A ( I PDFP ) 
2038808 A ( IPDFP 1 ) . A ( IPDI S 1 ) . A( IPDFC ) , A( IPTN ) , A ( IPOTN ) , A ( IPSKEW) ’ 
2038908 A( IPTSKW) ,A( IPXYZ ) ,A( IPTHKN ) ,A( IPFLAG) ,A( IPELVO ) ,A( IPTEM 
2039008 A(IPTMAT),A(IPORTH) . A ( IPSXSX ) , DELT ) 

2039 IOC 

203920 800 CONTINUE 
203930C 

203940 CALL UPDATE 
203950C 

203960 I F ( KDI SP .LT.O .AND. KSTRES.LT.O) GO TO 850 
203970 CALL OUTPUT( INC . LODINC ) 

203990 850 CONTINUE 
204000C 

204010 IFIRST-1 
204020C 

204030 I F ( NTEMCH . NE . O ) GO TO 900 
204040C 

204050 1ST I FF=0 
204060C 

204070 900 CONTINUE 

204080C * 

204090C INITIALIZATIONS FOR THE NEXT LOAD CASE 

204100C 

204110 I STI FF=0 

204120 DO 910 JEL= 1 , NUMEL 

204130 I EL= JEL 

204140 CALL FLAG ( 2 , 1 , IDSTI F . 0 . I A ( I PF LGB+ 1 EL ) ) 

204150 910 CONTINUE 
204160C 

204170 GO TO 200 
204180C 

204190 990 CONTINUE 
204200 STOP 
204210 END 
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APPENDIX C 


LISTIN G OF INPUT DATA 

Following is a detailed description of the input for the 8 
noded shell programs and 20 noded brick programs, input for the 9 
noded shell program is similar to that for the 8 noded shell. 
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INPUT "OR M$S8 WITH H-A PLA$T 


CCP MS55 w:TH *’ a and simple creep 

:ccc SwSrdltine ma inpg 


4 C 4 CC 
4 C5CC 
1 06 CC 
4 C70C 
4 080C 
109 CC 
1 1C0C 
4 1 IOC 
1 1 ICC 
4 4 3CC 

• 4 4CC 

• 4 £0C 
i 16CC 
1 170C 
1 180C 
1 1 9CC 
1200C 
12 IOC 
1220C 
i 230C 
1240C 
125CC 
126CC 
1270C 

• 2 8 C C 

4 Q nr 

4 300C 
4 3 IOC 
1320C 
133CC 
‘34CC 
4 350C 
1360C 

4 c^cc 

1 38CC 
4 39CC 
‘40C C 
M IOC 
4 -2CC 
1430C 
1 44QC 
1450C 
1460C 

14 7 OC 
1480C 
1 4 90C 
4 500C 

15 IOC 
1 520C 
1530C 
1540C 
1550C 
1560C 
1570C 
1580C 
1590C 
1600C 

16 IOC 
1620C 
1630C 
1640C 
1650C 
1660C 


ITHERM : 


NMAT : 
NMT : 

NPS : 
NMTCRP : 
NORTHO : 


NOTE -he f 0l wOW I NG INPUTS ARE AlL ON ONE LINE FOLLOWING ThE TIT..E 
I ana l ' analysis flag 

-^t.C r eao aaoi t iona1 line containing restart in*o 
= i (in absolute value) elastic 

s 2 ( in aosolute value) clastic 

=3 ( in aosolute value) creep 

_ ni . _ 2:3 Mn absolute value) c^eec ana plasticity comomec 

tre- reference temperature 
I3CTHE; isothermal element * \ ag 
=0 isothermal 

_ _ _ •» . . ne C snao ® functions used to interpolate int.pt. temperatures 
ITHERM: thermal load flag 

*0 no thermal loads 
. ne . 0 thermal loads calculated 
NMAT: numDer of materials 

NMT: numoer of temperatures for material property spec i f i ca t i on 

NPS: numoer of points used to specify stress-strain curves 

NMTCRP: numoer of temperatures where creep coefficients specified 
NORThO: orthotroo i c material flag 
*0 isotropic material 

.ne C ortnotropic material imat’l props, in 3 ai-ecticns - eG a i 
Program stores numoer of different ortno mat Is in NORTEC 
nskew snewea coordinate s, stems 'lag 
=0 no skewed coordinate systems 
ne 0 Skewed coordinate systems 

Program stores numoer of Skewed coord, systems in NSKEw 

mas SCO- mass matrix flag 

*0 no mass matrix calculated 
■i lumoed mass matrix calculated 
*2 consistent mass matrix calculated 
KFivj T f: fixed node flag 
*C no f i xea noaes 
. ne . 0 tnere are some fixed noaes 
System of eans is ^educed Oy numoer o' fixec DOF 
KPDIS prescrioea nodal displacement flag 
KSkwBC . skewed oounaary condition flag 
*0 no skewed BC 
. ne . 0 skewed BC ' s 

KLOCAL : local coordinate system flag 
*C no local coord, systems 

.ne.O there are local coord, systems (5 max) 

KPRCNT : nodal percentage distribution flag 
■0 no percentage di str iput ions 

.ne.O there are percentage distributions of nodes (5 sets max) 
MIDXyZ: flag to generate element mid side nodes 
•0 no mid side node generation 
.ne.O mid side node generation 
kxyzcO: nonstandard coordinate input flag 
*0 standard coordinates used 
.ne.O nonstandard coordinate order is used 
Standard is (X.y.Z) for cartesian 

(R. THETA, Z) for cylindrical 
( P . THETA , PHI ) for spherical 
KGAUSP : gauss integration order flag 
■0 first order integration 
.ne.O second order integration 

Wi.W2.W3: unit vector defining axis of rotation in global system 
G1.G2.G3: unit vector defining gravity direction in global system 
XBAR AX , XBARA Y , XBARAZ : location of origin of rotational coord, 
system in global system 

KEIGTP: eigenvalue/eigenvector computation flag 


MASSCD ■ 


KFlXXTt 


KPDI s 
KSkwBC . 


KPRCNT : 


MIDXYZ : 


KXYZCO : 
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INPUT FOR MSS8 WITH H-A P L AST 


* € TOC 

1 c a cc 

159CC 
■ T COC 


date -10*15. ’85 ^ I ME = 16 

*C ano.MASSCD ne C calculate e i genva 1 ue - e . oenvec t = - 

us’no aete-mmar-t searcr n e r-c 5 

and MmS sw.ne.c calculate e-genvaiue eigenvector 
using suosoace iteration metnoa 


k io: 5 


1 

• T3CC I h I I ANAL r . LT 0 ) ------ - 

* T4CC NOUT . lASCAS . NEXCAS 

1 7 50C 

,’^ 70 - N0UT: ou *Bux < f or later -estart) file creation flag 

’ 111 0 ° outBu ' t *’>e created y 

,-a?r . .a-.r’ 1 , 0utput f ’'* f ° r "* star * created 

•annr ^ ^ oaa =ase on r e=tart file ‘-cm wmcn solut'c- or-ceeas 

■ 800C NMCAS. ioaa case on .nout file .men is next load case to sc’/e 

1820C 

1830C SUBROUTINE ADRESS 

1840C* 

1850C SUBROUTINE EKIN02 * 

1860C 

1 870C I EL 1 . I E L2 . IEL3 . 

1880C 0 

1890C 

1900C IElh; element numoers 

!«2§ In enterea as native, all elements ‘-om orevious 

, 930 ; ® n ' e ea e,emenT *=> me negative element are included 

■ 940C ' • 

*35cc Subroutine nddopo 

1 960C 

1 9 7 0C IDOFN . NODE 1 , N0DE2 . N0DE3 

1 98CC C ■ 

1990C 

2 OOCC IDOFN OOP D e<- noae 

20 ICC NGGEn nodes with this number of DOF 

*C2CC 

2C30C - 

20U0C SUBROUTINE CONNEC 

2060C SUBROUTINE ELGEN2 ’*"*’** 

I E- . t INODE ( I ) , I « 1 .NODE ) . I E END 1 . I E I NC 1 . INI NO 1 . IEEND2 . IEINC2 . ININC2 


2080 C 
209CC 
2 1C0C 
2 110 
2 1 2CC 
2 1 30C 
2 140C 
2 i 50C 
2160C 
2 17QC 
2 180C 
2 190C 
2200C 
22 IOC 1 
2220C 
22 300 
2240C 
2250C 
2260C 
2270C 
2280C 
2290C 
2300C 
2310C 
2320C 
2330C 


IEL : 
INODE : 

IEEND1 

IEINC1 

ININC1 

IEEN02 

IEINC2 

ININC2 


Deg inning element number 

array sized to numoer of nodes oer element and is 
connectivity of element in terms of g local node numbers 
ending element number for generation set 1 
increment in element numoer for generation set i 
1 ncrement in node numbers for generation set i 
ending element number for generation set 2 
increment in element numoer for generation set 2 
increment in node numoers for generation set 2 


SUBROUTINE X Y2C0R 
I P ( KLOCAL . NE . 0 ) 

X Y ZO( 1 , NLOCAL ) , XYZ0( 2 . NLOCAL ) . X Y 20( 3 , NLOCAL ) 

0 

XY20: array 3 by number of local coordinate systems and is 

origin of local cooramate system in global system (5 max) 

I F ( KPRCNT . NE . 0 ) 

NPTS. ( PERCNT ( I , NPRCNT ) . I - i . NPT5 ) 

0 
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t : t w 5 = 
2 34CC 

INPUT FOR MSS 8 k 

2 3 r : 

NP - S '-umpor 

235:c 

PERCNT 

a r r a 7 

- 2 . C 
238CC 


d i s t r 

239CC 

: - 1 k x. zc c 

NE 0 » 

24C0C 
24 iQC 

:co>> z< i 

MCO 

24ZCC 
2 4 30C 

ICOXYZ 

array 

2440 C 

IF( I NCI . EO 

.0 ) 

2450C 

NODENC.X 

i . X2 . i 

2 4 6 CC 
2470C 

0 



" A LAST OATS s iq. «5 35 

:aes aer,e r a:ea *n *^ s eo* 


T - V E - ’6 * 3 


3 ^:icn sets arc 


w<V 2 array specifying coordinate input oraer 


2480C 

2490C 

2500C 

25 IOC 

2520 C 

2530C 

2540C 

255CC 

2560C 

2570C 

2580C 

2590C 

2 cCOC 

25 ICC 

2620C 

2 6 3GC 

2540C 

265CC 

2660C 

2S‘ , CC 

2680C 

269CC 

2 -QCC 

2" IOC 

2“* 20C ' ■ 

2 " 3CC 

2 ~ 40C 

2 T 5CC 

2 7 €0C 

2-70 C 

2780C 

2790C 

2800C 

2810C 

2820C- ‘ 

2820C 

2840C ! 

2850C 

2860C 

2 S 7 CC 

2880C 

2890C 

2900C 

29 IOC 

2920C* * ■ 

2930C 

2940C S 

2950C 

2960C 

2970C 

2980C 

2990C 

3000C 


NOD ENO : ^ Deg i nr, i ng noae numoer 
*1.X2.X3. coordinates 

INC1 ; noae numoer increment for generation set 1 
rtvewn n ?°* ? umoar ' ncrement for generation set 2 

!SS5; ;r Sii’SSlSVKS.-SS' — <— ■ 

s 0 cartesian 
* i Cvli nor i ca 1 
- 2 sonerical 

I PCNf 2 n!^* nta9e a!!tr ’ Bu,,Bn set ‘°n noa e generation set - 
IPCNT2: percentage a, str ,cut - on set for noae generation set 2 

I p < INC i . ME 0 ) 

NODE NO . » 1 . x 2 . * 3 

o D x 5 °xi, ena ’ n9 n ° ae " UmDer for noae veneration set • 
i , ** , a J . enoing coordinates 

I c ( INC2 . NE . O) 

NODENO, x 1 , K2 . \ 3 

x° D x-°^ ena ’ n9 n °°* numDer noae generation set 2 
A • ending coordinates 


; s ( MID* YZ . NE .0 1 

SUBROUTINE NOOMID 

IEL 1 . IEL2. IEL3. . . . 
0 


IELn: elements for mid-side node generation 

subroutine nooein 


TH I CK , NO 1 , N02 . N03 , 
0 


THICK: thickness 

NOn : nodes wi tn this thickness 


SUBROUTINE mtable 

SPECWT . ISOCOD. TMAT( 1 . NMAT ) , TMAT ( 2 , NMAT ) TMAT ( NMT . NMAT ) 

s oec ,f1 c we’gnt * aensity * g (LB/CUBIC INCH) 
ISOCOD: i sot rop 1 c material code 
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TITLE 


3C30C 
3040C 
3C5CC 
3060C 
3C7QC 
3080C 
3090C 
3100C 
3 1 IOC 
3 120C 
3 130C 
314QC 
3 150C 
3 160C 
3170C 
3180C 
3 190C 
3200C 

32 IOC 
3220C 
3230C 
3240C 
3250C 
3260C 
327QC 
3280C 
3290C 
33COC 

33 IOC 
3320C 
3330C 
3340C 
3350C 
3360C 
3370 C 
3380C 
3390C 

34 00C 
34 IOC 
3420C 

34 30C 
3440C 
3450C 
3460C “ ' 
3470C 
3480C 
3490C 
3500C 

35 IOC 
3520C 
3530C 
3540C 
3550C 
3560C 
3570C 
3580C 
3590C 
3600C 

36 IOC 
3620C * • 
3630C 
3640C 
3650C 
3660C 
3670C 


a INPUT FOP M5S6 WITH H-A PLAST DA^E - i C 15- 85 T’ME = ■*€ * 

= 0 ISCt'-ODiC 
ne C o^tncr-oc i c 

• MA , array numoer c f mat 1 temps o, numoer o* mat : s arc *s 
temperatures wnere material properties are SDec'f -ec 

I p l ISOCGD.EO.O) INPUT NMT CF THE FOLLOWING LINES 
E 1 . NU 1 . ALPHA 

El; Young s moau 1 us ( program multiplies o> 1E6) 

NU 1 : Poisson • S rat l O 

Alpha thermal coefficient (program multiplies py IE-6) 

1*1 ISOCOD. NE .0) INPUT NMT OF THE FOLLOWING lINcS 
E 1 . E2 . E3 .NU1 .NU2.NU3. ALPHA 

En; Young 5 modulus in the n direction (multiplies by 1E6> 

NUn : Poisson's ratio in the n direction 

alpha: thermal coefficient (multiplies by IE-6) 

IF(NPS.NE.C) INPUT NMT OF THE FOLLOWING LINES 
STSXSXl 1.1 . NMT , NMAT ) . STSXSX( 2. 1 . NMT . NMAT ) 

STSXSX( 1 .2 .NMT. NMAT ) . STSXSX ( 2 . 2 . NMT . NMAT ) 


S T Sx Sx I 1 .NOS. NMT .NMAT ) . S T S X S X ( 2 . NP S . NMT . NMA T 1 
e ETA ( 1 , NMAT ) , BETA ( 2 . NMAT ) , . . . E E T A ( NMT . NMA T ) 

STSXSx: array 2 by NPS by NMT p. NMA * and is 

Stress , st ra i n pairs defining curve at temp NMT 
BETA ar^ay NMT by NMAT and is 

hardening coefficient at temperature NMT 
*1 isotropic hardening only 
=0 Kinematic hardening only 
(C .LE. BETA .lE. 1) 

I F ( NMTCR P NE O) INPUT NMT OF THE FOLLOWING ' TNES 
T EMP . SNGRM . C . R . STRCUT 

1 EMP : temperature at which these properties apply 

SNORM normalizing stress for tnese properties 
0: creep property 

R ; creep prooerty (creep « 0 - (stress) •• R ) 

STRCUT ; cutoff stress below which no creep occurs 


IF(NSKEW.NE.O) 
subroutine SKEW 

ISKEWN. A1 .A2.A3 
0 

ISKEWN: skew set number 

gt.O A ' s are node numbers defining skewed cooro. system 
.lt.0 A's are successive rotation angles in degrees 
Ai.A2.A3: nodes or angles defining skewed coordinate system 
if nodes, axisi is formed from ai to A2 
axis2 is formed from AI to A3 
ax i s3 is formed as axisi cross axis3 


IFfKSKWBC .NE .0) 

SUBROUTINE SKEWBC 

ISKEWN. NODE 1 . N00E2 , NODE 3 . 
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.£= INPUT FOR MS S3 WITH h-A P L AS~ 


DATE = IQ/ i 5 85 


’ I ME s 1 6 *6 


O ~ J 1 *. 

3 “* ’ OC 
3t::c 
3730C ' 
3 7 40C 
3 7 5CC 
3760C 

3773c 

3?80C 

3790C 

3SC0C 

36 IOC 

3 8 20C 

3330C 

384CC- 

3850C 

3860C 

38 TOC 
3880C 
3890C 
390CC 

39 IOC 
3920C 
3930C 
39-iOC 
3950C * 
3960C 


I5k?wN. st<ew set numpe r 
NODE ^ nooes wnere skew 3C s 


referrea to tnis Skew set 


a'-e acc e: 


I ? ( N0R7HC NE . 0 ) 

subroutine orthop 

I SKEWN . I EL 1 , IEL2 . I El3 . 


ISKSWN. Skew set numoer 

elements wnose material prooerties nave this skew 


IFIKFMTY . NE . C ) 

SUBROUTINE FIXITY 

IDIREC . NODE 1 . NODE 2 . NODE 3 . 

0 


IDIREC; giooal f 1 xea direction 
NODE n : nodes wtn tms fixitv 


39-oc 

I c i KPDIS NE .01 


■ \ 

3980C 

SUBROUTINE PREDIS 



3990C 




40C0C 

IDIREC . VALUE . NOD? 1 , 

NCDE2 . NCDE3 , 

n 

40 ICC 

0 


4C20C 



■ 1 

1 

4030C 

IDIREC: global a * - 

ection of prescrioed displacement 


4C4QC 

VALJE: magnitude 

of orescrioea displacement 


4C50C 

NODE n • nodes with 

tnis displacement 

' — 1 

4C60C 



! 

407 OC - 



i 

4C80C 





SUBROUTINE LOADS 

KD I SP . KSTRES . LODI NC , MATCH . NTEMCH . NR E MOV , NO LOAD . RPM . 
KGLOAO . KTLOAD , KALOAO . KLLOAD . KCLOAD , 

method, epstol.toleps, TOLD IS 


KDISP: 


KSTRES: 


LODINC; 


nodal data print flag 
■1 print noaal data (displacements) 
as specified in this load case 
*0 print nodal data as specified previously 
■-1 no nodal data printing 
element data print flag 
■1 print element data ( stresses , stra 1 ns ) 
as specified in this load case 
“O print element data as specified previously 
*-1 no element data printing 
number of load increments 
• 1 1 . 0 sets ITIME • 1 (reads time data) 
material change flag 
■1 there is material changed 
*0 no material is changed 
termperature cnange flag 
*1 there is change in temperature 
■0 no temperature change 
element removal flag 
*1 some elements are to be removed 
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■LE* INPUT POP MSS8 WITH SIMPLE P u 


J 

I 

30 1 0C 
3 C 2 0 C 

ALPHA thermal coeffi 


30 30C 

: r i r sococ ne o > inpu~ n 


oC*J 

E * .E2.E3.NL* . NL'2 . NU 3 . 

f 

3C50C 



306CC 

you^g s moauus 


3070C 

Nun Pc 1 S SOn s r a t ’ O 


308 OC 

a w p-ia t ne r ma 1 coeffm 


3C9CC 


f 

3 t OCC 

IF1NPS NE Cl INPUT NMT ( 


3 1 1 OC 

$~$ * S < ( 1 . 1 . NMT , NMAT ) . 1 


3 1 2CC 
313CC 

S~ S x S * 1 1.2. NMT . NMAT ) < 

j — - 

3*400 


j 

3 1 50C 


) 

3 160C 

S"SX$X( 1. NPS. NMT. NMAT 


OAT£ = iQ •' *5 8 5 


i£SE _ I \E 5 


" :m£ = 
^E-6 > 


■ £c i 


3 1 TOC 
3 180C 
3 i 90 C 
32C0C 

32 IOC 

J 4. *. w w 

3230C 
3240C 
3250C 
326CC 
327CC 
3 2 8 0 C 
3290C 
3 3 OCC * 

33 IOC 
332CC 
3330C 
334CC 
3350C 
3360C 


STSXSx- array 2 Dv NPS Dy NMT Dy NMAT ana is 

stress , st r a ir pairs defining curve at temp. NM"! 

I c ( NMTCRP NE 0) INPU T NMT OF THESE „INE$ 

T tMP . SNORM. 0 . R . STPCUT 

T EMP- temperature at wm:r tnese properties app 1 , 
SNORM- norma 1 : 2 i ng stress *o r tnese pnooenties 
C creep property 

creep property ( c-eec = z * 'stress) p ) 
STpcu T Cutoff stress oeiow *ricn nc creep occurs 


I F t NSKEW NE C ) 
S'JEPCuTINE SKEW 

ISKEwN. At .42, A3 



j 


337CO 
338CC 
3390C 
3 4C0C 
34 1 0C 
34 20C 
34 3CC 

34 50C 
3460C 
347CC 
348CC 
349CC 
3500C 

35 IOC 
3520C 
3530C 
3540C 
3550C 
3560C 
3570C- 
3580C 
3590C 
3600C 

36 IOC 
3620C 
3630C 
3640C 
3650C 
3660C 
3670C 


ISkewn skew set numoe r 

gt.C As are noae numoe^s aef*n<ng sne^ea :oora s.stem 
w * v - / ^ 5 are successive o t a 1 1 on angles in oeo r ees 

‘•-.A*. A3, nooes or angles aef imng skewec coorairate system 
noaes. axiji is formed Ai tc A2 

axis2 is tf ormea ai to A3 

axis3 is formed as axisi cross ax : «3 


I F ( K$KWBC .NE .0) 

Subroutine skewbc 

I SKEWN . NODE 1 , N0DE2 . N00E3 . 

0 

ISKEWN: skew set numoer 

NODEn: noaes where skew 5C s referred to this skew set are applied 


:f<northo.ne o> 

SUBROUTINE OPTHOP 

ISKEWN. IEL1 . I EL2 . IEL3 

0 

ISKEWN: skew set number 

lEun: elements whose material properties have this skew set omenta 
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' X T L E * INPUT FOR MSS 8 WITH SIMPS-H »u 


DATE = i C ' 1 5 85 


* I WE = 15 2 3 


3680C 


36 9CC 


3 'COO 

: p 1 KF I X 17 

37 iCC 

S ~ERG0 T I NE 

3 7 20C 


3~30C 

I 0 I R E Z . 

3 7c cc 

0 

375CC 


3 7 sec 

IDIREC 

37“>CC 

NODE n 

3780C 


3^900 • 


38COC 


38 IOC 

: c 1 k pq 1 s 

3820C 

subroutine 

3830C 


38J0C 

IDIREC. 

3850C 

0 

3860C 


387QC 

IDIREC 

3880C 

VALUE 

3890C 

NCOEn 

39CCC 


39 IOC- 


392CC 


3930C 

SUBROUTINE 1 

394QC 


3950C 

KOI $P . K! 

3960C 

K G i 

397QC 

ME' 

3980C 


39S0C 

KOI $3 

4Q0CC 


- c 1 cc 


4C20C 


-1C 30C 


4C-0C 

KS'RES 

405 DC 


4 06 OC 


4C7CC 


4 08 OC 


4CS0C 

LOOINC : 

4 iCOC 

MATCH 

4 1 ’OC 


4 12 OC 


4 1 3CC 

NTEMCH : 

4 1 4 CC 


4150C 


4 1 60C 

NREMOV : 

4 17 OC 


4 1 8 OC 


4 1 90C 

NCLOAD: 

4 200C 


42 10 c 


4 2 20C 

PPM : 

4230C 

TEND: 

4240C 


4 250C 

KGLOAD : 

4260C 


4270C 


4280C 

ktload 

4 290C 


4 300C 


4310C 

KALOAD : 

4320C 


4330C 


4340C 

KLLOAD : 


■ *ec ai f, ection 
rn t n > s f 'xit, 


NE C > 
BREDIS 


. N0DE2 , NODE 3 , 


gisDai a^ection cf prescribed displacement 
magnitude cf orescr - oea d’sdacement 
nodes wi tn tr ' s 31 sc ! acement 


■ *00 .E^S 3w,’Cw£PS,’0.::S 

r ooa > data e- r -t * ' ag 
= i or t nt n cca i data i d i so : acemen t s l 
as soec ' f 'ec ' n t**'s cad case 
® r * n 7 r ooa • data as soec i * iea oreviousl . 

= -i no ^ccai data c r •''ting 
e'eme^t data print 4 ' ag 
2 * ° r ip ' element data (stresses , strams ) 
as SDec i * * ed »*■ tn-s oaa case 
= 0 Dr m: element data as soecif ied previous 1 , 

= - i no element data c ■ ' ng 
n umpe p ‘ c* load increments 
material change 4 iaa 
=i there 1 s material changed 
*0 no material is cnanaed 
termperature change * 1 ag 
=1 there is change in temoerature 
=0 no temperature cnange 
element removal flag 

some elements are to oe removed 
*0 no elements removed 
concentrated load * 1 ag 
«i there are concent-area l oaas 
s 0 no concentrated ’oaas 

revolutions per mi%te *or this load case 
time at ena of i oaa step 

used only when time oeoenaent response (creep) desired 

gravity load flag 

*i there are gravit, loads 

*0 no gravity loads 

thermal load flag 

=i thermal loads calculated 

*0 no thermal loacs 

area load flag 

*i there are area loads 

=0 no area loads 

line 1 oad f 1 ag 
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TITLE= INPUT for MSS 8 with SIMPLE PL DATE *10'15‘85 


'IME * 16 23 


4350C 

436DC 

42' , 3C 

438CC 

4 39CC 

440CC * 

44 IOC 

4420C 

4430C 

4440C 

4450C 

4460C 

4 4 7 DC 

4480C 

4490C 

4 5COC 

4 5 iQC 

4520C 

453CC 

4540C- 

4550C 

4560C 

4 5TOC 

4560C 

459DC 

4600C 

46 tOC 

4620C 

4630C 

4640C 

4650C 

4660C 


* 1 tnere are 1 me i oacs 
= 0 no ’ me i cacs 

I5 LE ? 5 convergence tolerance c- st-a-n 
OlCIS convergence tolerance on 3 i SO 1 acement 


: f I KD : SR . 50 . 1 ) 

subrou t :ne outputn 

I0ISP. I RE F , NO 1 , NC2 . NC3 , 


1 2 I SR pr i n t flag 

*C aon • ermt 
* 1 ao print 

I RE c 

NOn noaes for wnicr, output is aesirea 


:p(ks t res . eo n 

SUBROUTINE CUTPuTN 

IElF . ISTRESS. ISTRAIN. LOCA I c ■ i 
C 


" - ' *orce print * ) ag 

= 1 o r t n t element forces 
=C aon t print-fences 
^STRESS' stress print f>ag 


J6-3C 

*0 aon 

t c-mt s tresses 

4 6 8CC 

I 3TR A I N • Strair 

c •" 1 n t c ' a g 

4 69DC 

= t pr 1 1 

nt strains 

4 -DOC 

=0 aon 

t ormt ST-?’ ns 

4 7 ICC 

lD S a 1 oca 1 1 cn 

for whicn aata 1 

4 - 2CC 

* i cente 1 

of gravit. of e 

4 - 3 DC 

=2 e'eme 1 

n t noaes 

4 — 4 DC 

“ 3 Gauss 

points 

4^50C 

I£l~ elements 

for v, men output 

4 ”6 DC 




4 “* 9 DC 
4 79GC 
480CC 

48 tOC 
4 S 2 DC 
4 8 30C 
484CC 
4 850C 
4860C 
4870C 
4880C- 
4890C 
4 900C 

49 1 0C 
4920C 
4930C 
4940C 
4950C 
4 960C 
4S70C 
4980C 
4990C 
5000C * 

50 IOC 


I c I MA^CH NE . 0 ) 

subroutine input2 

MATNO . I EL 1 . I E l2 . I E L 3 . 

0 

MATNO : new material number 

IELn: elements cnangea to tms material 


I F ( NTEMCH . NE .0) 
SUBROUTINE NOOEIN 

TEMP . NO t , N02 . NOD , 


T EMP : new temoerature 

NOn: nodes cnangea to this temperature 
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" L E = INPUT FOR MSS 8 WITH SIMPLE 


C A ~ E = 4 0 ’595 


':me = is :3 


5 22 OC 
5C30C 
5COC 
50500 
506 OC 

c C ~ : 

5CSCC 
5 09 OC 

5 ?:oc 

5 ' 'OC * 
3 4 20C 
5 i 3QC 
5 * -OC 
5 4 50C 
5 1600 
5 4 7 OC 
5 1 8QC 
5 1 90C 
5 2 OOC 

52 4 OC m 
5220C 
323DC 

5 2 -OC 
5 2 5 0 C 
5260C 

5 2 SCO 
5290C 
5 3COC 

53 1 0C 
5320C 
53300 ' ’ 
534 OC 

5 36CC 
oG^OC 
5 3SOC - - 
5 39 OC 

5 4 * OC 

54 0 OC 
3 4 300 


I F I NT EMCH NE 0 I 

su 6 rou t ine nodein 

~~Emd .NCI, NC2 . NC 3 . 


= temperature g^ao’ent 

n w aes ■< 1 tr t'v s te r nDe r a*. j r e g^ac'ent 


I ‘ ( NREMCV NE o i 
SUBROUTINE INPUT 1 

IE- 4 . I EL2 . IEl2 . 

removea e’emerts 


I c < NC ^0 AD NE . 0 > 

Su5kOU t !NE clcads 

ID I RE 0 . VALUE .NC 4 .NC 2 .NC 3 . 


IDIR-C. glopai a ,r ectio n o* : caa 
v A w w £ ■ magnitude c f 1 caa 

NOn noaes wne^e s :once^tratec 1 oaa apoiiec 


: ' k E _ 0 AC NE 0 * 
IROU^INE E LO AOS 


I s ' KA. 0 A 0 .NE 0 i 
SUBROUTINE A LOADS 


i=-j 


5 4 — j\. 

g 

5450C 


5 4 6GC 

I EBEG 

5470C 

I £ E NO . 

548GC 

IEINC: 

5J9DC 

I c ACE 

550CC 

1 0 I RE 

55 4 OC 

P 1 . P2 . 

5520C 

553CC'- 


55JOC 


5550C 

IFfKLLOA 

5560C 

SUBR0UTI 

557QC 


5580C 

I EBEG . 

5 590C 

0 

560CC 


56 IOC 

I EBEG : 

5620C 

I E END . 

5630C 

IEINC : 

5640C 

IFACE : 

5650C 

I D I R : 

5660C 

P 1 . P2 : 

5670C 


5680C * » 



EE EG. I E END . IEINC , IFACE . ID I RE . P 4 .P2.P3.PJ 

Deg inning element 
ending element 
increment in element numoers 


: i on 


Deg inning element 
enamg element 

increment in element nunoers 
element face wnere l me load aopl ied 
direction of load in element local coord, system 
pressures at either ena of line load 
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*: t L£= INPUT FOR MSS8 WITH simple PL 


5- ‘DC 

5’3CC 

5~-OC 

s’soc 
5 - 6 OC 
5^0C 
5 "■ 5 DC 
5*^900 
5B00C 
58 1 CC 
5S20C 
58 30C 
5840C- 


3ATE =1C"5 85 


-IME = 16 23 


p A G £ = 5 


3 9CC' 


KPDIS NE c ANC 
S^ERCU’INE PREDI 5 


:a$e ne i i 


"EC . vii.Lt . MODE 1 . NODE 2 . NODE 3 , 


ICIREC. g>ccal a'^ect -cn c? orescnoea a ’ sc ’• acene^* 
vALwE maan;tuae o* cresc^ioea aisoiacement 
NOOEr. noaes wtn :h, s a i so l acemen t 

\0"'E -HE NOCES with PRESCRIBED DISPLACEMENTS MUS’ SE ’he 

same Throughout and must be specified from ’«e cutset 
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:tle> input POR MSS8 with bcdnep UN 


DATE * i C ' < 5 8 5 


r IME s ■ 15 : : 


i OGO c * 
4 : *oc 

.2CC- 
’ 2 300 
♦c JCC 
1C50C 
1 06 2 C 

*:^cc 

4 08CC 
1C90C 

* iCDC 
1 t ICC 
t * 20C 
i 4 DOC 
i i -CC 
1 1 5 CC 
1 1 60C 
1 170C 

1 ^eoc 

1 19CC 
1200C 

* 2 ' OC 
1220C 
123CC 
1 2 -ACC 
125CC 
*2600 

* 2 7 CC 

* 28CC 

* 290C 
130 CC 
*3 iCC 
1 32CC 
1330C 

* 3 -i C C 
13502 

* 360C 
i 3 T CC 

* 38CC 
1 390C 


i 4 2 CC 
i - 3CC 
1440C 
1 4 5 OC 
146.CC 
1 470C 
1480C 

* 4 90C 
1 5CCC 

15 IOC 
1 520C 

* 530C 
*540C 
1550C 
1560C 
1570C 
1580C 
1590C 
1600C 

16 IOC 
1620C 
1630C 
1640C 
1650C 
1660C 


INPUT FOR M$S8 WITH BCDNER 5 CONST : - jr * yc MODEL 
SUBROUTINE MAINPG 


:*E T-c 'Ol-OWING are A_ on. ONE ^INE A-^EP ’-4£ - : t_E ^ :N(£ 

: a\a w • ana 1 . s i s flag 

t ~ ~eac aaa i t * :na ’ line contain i no r 6 st 3 r t ’^*0 
= < t ' " aesdute value) elastic 
= 2 ( • r aesdute value) inelastic 
reference temoerarjre 
1 so t nerma 1 element # laa 
=C 'setne^mai 

n e C sraoe * unctions usea tc -"te-Doiate ct 
tnerma 1 ioaa f • ag 
~ C no the r ina.l i caas 
ne C thermal l oaas calculated 
numoer of materials 

numoer c* temperatures mate^al property SDeci^cafc^ 
numoer c 4 temperatures wie r e creep coefficients SDec ’ ^ ■ ec 
or t net roD i c material * • ag 
*0 isot-'OD'C materia 1 

ne . v or t nc t '"CD 1 c material ( mat 1 proos . in 3 directions -ea d> 


tre- 
: so t he 


: t -erm 


NMA T 

NM T 

NMTCRP 

NORTHQ 


temDe r atw,'*es 


NSKEw 


MAS SCO 




KPOIS • 
KSKWBC . 


KLOCAL 


kprcnt • 


MIOX V Z 


KX> ZCO : 


KGAUSP : 


W1.W2.W3 
G 1 .G2.G3 
XBARAX . X 

KEIGTP : 


o r tnc mat' is in N0R""-3 


:cord systems m NSKEw 


• xec Our 


3r og'*afr steles numoer cf different 
skewed coordinate systems flag 
: 2 nc skewed coordinate systems 
ne 2 Skewed coordinate systems 
0, "cgram stores numoer c* skewed 
ma ss mat r 'x flag 

= 0 nc mass ma t r 1 v v , ' l ce computed 
r i . umoec mass matrix computed 
=2 ccns'sten* mass matrix computed 
* 1 *ea ncae * i ag 
= 0 r c * * xea nodes 
ne 0 tne r e a^e some f ucd noaes 
System of eons ^s reduced C5v numDGr 0 < f , 
prescr^oed noaa i displacement flag 
skewed oounaary conoit’on * ■ ag 
0 nc skewed EC 
ne C skewed BC s 
local coordinate system flag 
=0 no 'ccal coord s. stems 
ne . C there are local ccora systems (5 max) 
noaal percentage distriputicr *iag 
=0 no percentage distributions 

ne.C tnere are percentage distributions of nodes (5 sets max) 
f * ag to generate element mid side nodes 
=C no mid side node generation 
ne 0 mid side node generation 
nonstandard coordinate' input flag 
*0 standard coordinates used 
ne . 0 nonstandard coordinate order is used 
Standard is (X.v.Z) for cartesian 

(R. THETA. Z) for cylindrical 
( R. f theta . phi ) for spherical 
gauss integration oroer flag 
s C first order integration 
ne 0 second oroer integration 

; unit vector defining axis of rotation in glopal system 
: unit vector defining gravity direction in glooal system 
BARA v , XBARA2 : location of origin of rotational cooro. 

system in global system 
e i genva lue/ei genvector computation f l ag 
«C . and . MASSCD . ne 0 calculate ei genva Vue/ei genvector 
using determinant search method 
* i . and . MASSCD . ne . 0 calculate e i genva l ue/e i genvector 
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*! T LE* INPUT POR MSS8 WITH SOONER UN DATS ='0 '5,85 TIME * '6 D5 

- °'? c using sucscace iteration metres 

* -3 3 OC 

1690C - 

4 "DCC SuBRCU T INE AOQESS 


' 2CC SUBROUTINE E kino 2 
' 30C 

:el? . lEui. iel3. 


element numoers 

I* element enterea as negative, an elements ‘-sip orevous 
erterea element to the negative element are mcluaec 


•'■iOC 
<~5CC 
4 7 S0C 

^70C I E Ln • 

1 “*30C 
179CC 
4 800C 

18 IOC 

1820C SUBROUTINE NDDQF 2 
<830C 

1840C I DC C N , NODE 1 . NODE 2 . NODE 3 . 

1650C 0 

1860C 

1 8 7 OC I0C c N. DOF oer ncoe 

1880C NODEr: noaes witn s numoer of DOF 

1890C 

t 90CC 

1 5 1 OC SUBROUTINE CONNEC 

320C — 

193CC SUBROUTINE £ LGEN2 
1940C 

’SSOC I EL . ( INODE i I > . I * i . NODE ) . I EEND 1 . I E I NC 1 . ININC 1 . I SEND2 IEINC2 T NIN~2 

1 S60C C 

1970C 

I EL Beginning element n U noe r 

-NODE a ^ ^ a v sized to numoe'- cf neces oer element ana is 

connect i v i t > of e*eme-t *n terms of giccai n C ae n.umcers 
hoe , ‘ *or generation set 4 


2D 1 0C 

I EEND 1 . 

end i ng e ' 

ement 

222CC 

I E INC 1 • 

i ncrement 

t r e 

2 D 3 0 C 

ININC 1 . 

i ncrement 

1 r n( 

2040C 

IEEND2 ■ 

ending element 

2 0= DC 

IEIND2 

i ncrement 

i n e 

2 OcOC 
2 C~CC 
2C80C * 

ININC2 . 

i ncrement 

i n nc 

20SGC 

subroutine 

x > ZCOR 



2 iQOC 
2 1 IOC 
2 120C 
2 130C 
2 14 OC 
2 4 50C 
2160C 
2170C 
2180C 
2 190C 
2200C 

22 1 0C 
2220C 
2230C 
224CC 
2250C 
2260C 
2270C 
2280C 
2290C 
2300C 

23 IOC 
2320C 
2330C 


I c ( KLOCAL . NE . 0) 

XYZO( 1 , NLOCAL ) . XYZ0< 2 . NLCCAL ) . > v ZO( 3 . NLOCAL ) 


XYZO: array 3 by numoer of local coordinate systems ana i s 

origin of local cooramate system in global system (5 max) 

I F ( KPRCNT . NE .0) 

NPTS. ( PERCNT ( I , NPRCNT ) , 1= i .\ers ) 

0 

NPTS: number of noaes gen#r a tea m this set 

PERCNT: array NPTS x numoer of oercent distribution sets and is 
distribution of noaes oetween coordinates input pe 1 ow 

I F ( KX YZCD NE .0) 

ICDXYZI 1 ) . ICDXYZ< 2 ) . IC.DX YZ< 3 ) 

ICDXYZ: array specifying coordinate input order 
I F f I NCI . E0.0) 

NODENO . X1-.X2. X3. INC1 . INC2 . ISYSNO. ISYSTP, IPCNT 1 . IPCNT2 
0 
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TIT L£- INPUT FOR M$S8 WITH 60DNER UN 


0 A T c = i 


* 5 85 


->s 


:34cc 

235CC 

2360C 

23 TCC 
238CC 
: 39CC 
240CC 

24 ICC 
24CCC 
2 4 30C 
244QC 
2 4 SCC 
246CC 
2470C 
248CC 
2490C 
2500C 

25 IOC 
252CC 
2 5 30C 
2540C 
2 5 SCC 
256CC 
2570C 
2580C 
2590C* 
25C0C 

26 IOC 
262CC 
2 6 30C 
2640C 
2 c 5C C 
2 66 OC 
25 TOC 
268CC 
259CC* 
27CCC 

27 iQC 
2^202 
2 7 30C 
2 " 4QC 
2 7 5 OC 
27SOC 
2 ”7CC 

27 8CC 
2790C- 
28CCC 

28 1 OC 
2820C 
2B30C 
2840C 
2850C 
206OC 
2870C 
2880C 
2890C 
2900C 

29 ICC 
2920C 
2930C 
2940C 
2950C 
2960C 
2970C 
2980C 
2990C 
3000C 


’ : me 


NIOOENO 

A 1 . .< 2 . A 

INC * 

INC 2 
: S • 5 NO 

: s ■ stc 


aoeve coofci-.a: 


IPCNT 1 
I P CNT 2 


oegmn:ng node numoer 
3 coora snares 

noae numoer increment * c generation set 
R °ae numoer increment f'or generator, set 
local coora i nate system romper usee f c r 
9 1 °oa 1 cocra-.nate system tvoe 
S C cartesian 
= 1 Cv linar :;a' 

s 2 sonencai 

percentage a i s t r i ou t ’ on set ^cr noae generation set 
oercentage cistr 'cut ten set for noae generation set 


: ( INCi . NE 0) 
NCDENO . a 1 . *2 . x 2 


NODENO ending noae numDe*" *or noae generation set i 
X1 .*2. *3; enamg coora*nates 

1*1 I NC2 NE 0) 

NOCENO , * * , X2 . X3 

NODiNO enoing noae numoer fo" noae generation set 2 
y .i.*2.x3 enamg coorainates 


I F ( MIOx V 2 NE . C > 
SUBROUTINE NODMID 

. I E 12 . IEL3. 


elements mid-siae noae generation 


SUBROUTINE nodein 

T^I CK , NO 1 . NO 2 . NO 3 

n 

Tt “!CK tnickness 

NCn. noaes witn this thickness 


Subroutine mtable 

SPECWT . I SOCOD . TMAT M . NMAT ) . TMAT'l 2 , NMAT ) , . ,TMAT(NMT NMAT) 

0 

SPECWT : specific weight * density * g 

ISOCOD: isotropic material coae 

*0 isotropic 
.ne.O or t ho t rop i C 

TMAT : array numoer of mat 1 temDS Oy number of mat' is and is 

temperatures wnere material properties are specified 

I F ( ISOCOD. EO 0) INPUT NMT OF THESE LINES 
E 1 . NU 1 . ALPHA 

El: Voung's modulus (program multiplies Dy 1E6) 

NU 1 : Poisson's ratio 

ALPHA: tnermal coefficient (program multiplies Dy IE-6) 

IF( ISOCOD. NE .0) INPUT NMT OF THESE LINES 
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~:tle* input for MSS8 WITH sooner UN 


DA~E ■ 1C. 1585 


"ime * * 6 : = 




20 ICC 

2C20C 

3C3CC 

3C40C 

2050C 

3060C 

OC^OC 

3GBOC 

3C90C 

3 i OCC 

3 1 IOC 

3120C 

3 i 30C 

3 14QC 

O'SCC 

3 1 60C 

3 170C 

3 1 80C * 

3 i 90C 

320CC 

32 IOC 

3220C 

3230C 

324CC 

32 5 OC 

32SOC 

32 TOC 
3280C 
3290C 
33COC 

33 ICC 

3 3 n ^ 

22 3CC 

2 34QC ’ 

3 350C 
3 3 6 0 C 
3 3 "*CC 
33SCC 
3 39CC 

34 OCC 
34 iQC 
342CC 
34 30C 
3440C 
345CC • - 
346GC 
34 7 OC 
3480C 

34 90C 
3500C 

35 IOC 
3520C 
3530C 
3540C 
3550C 
3560C * • 
3570C 
3580C 
3590C 
3SOOC 

36 IOC 
3620C 
3630C 
3640C 
3650C 
3660C 
3670C • * 


£ 1 . E2 . E 3 . NU i . NL‘2 . NU3 . AL°ha 


t r 
♦ Nu r 
A w - 


■Cung s modulus 
~ c ' s s on s -at ’o 


* n fro n 


-'e-na 1 c oe* f i c • er 


a i ’ 
r>e n a t * 

i mu 1 t • 


5Ct ’ on I mu 1 t ipl i es 
ec T 1 on 

Cl 1 es Oy t E -6 1 


* E6 > 


I z ( NMTCRP NE Ci INPUT NMT CF Th£ 5 E ^INES 
TEMP 0 AN 2 C 21 22 AM A R 

t£mp * t emoe r 3 fur e m degrees ~ 

D. AN, ZO. 21 . 22 , AM. a, s = material Darameters reoui'etj 
4 - r Bearer's tioae 1 at temperature temp 

Zw.a. ,2-. in KSI. am t r * KSI. remaining are aime^S'Cn’ess 


:f(n$kew.ne o) 
SUBROUTINE SKEW 


: SKEWN . A 1 . A2 . A3 


I SKEWN 
A i , A2 . A3 


Skew set mumper 

gt.C Ac are rcae r umDer . s ae^mng Skewed cooro. svstem 
t ^ As a-e Successive rotation angles m degrees 
nodes or angles ae^'^ng skewed cco-'dnate SvStem 
if ncaes. axis 4 *s formed from A * to A2 
axis! ’S formed from A1 :o A3 
a>is3 ’s * cornea as axisi c-css axis3 


r p f KSKWEC NE oi 
SUBROUTINE SKEwBC 

: SKEWN. NODE 1 . NODE 2 . NODE 3 . 


IjK^wn. skew set -'umoen 

NOCm nodes wne-e Skew BC s referred tc xms Skew set are acc’-ec 


: r ( NORT ho nE 0) 

Subroutine crthop 

ISKEWN, I EL 1 , I EL2 , I E L3 
0 

ISKEWN: skew set number 

IELn: elements whose maternal orooer 1 1 es have this skew set omenta 


I f < kf i x it y . ne . O ) 

SUBROUTINE FIXITY 

IDIREC . NODE 1 . NOOE2 . NODE 3 , 

0 

IDIREC: global fixed direction 

NODEn: nodes with this fixity 
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INPUT FOP MSS 8 WlTh SOONER UN 


} A 7 £ = 10 / i 0 5 


' IME 


3680C 

3S9CC 

3TCCC 

3 n • oc 

3^2 CC 
3"3DC 
3 7 4 C C 
3 7 5CC 
3760C 
37-OC 

37 80C 
u^SDC- 
3800C 

38 4 DC 
3820C 
3830C 
3840C 
3850C 
3860C 

38 ~ DC 
3880C 
389CC 
3900C 

39 IOC 
3920C 
3930C 
3940C 
3950C 
3960C 

3 S 7 OC 
39S0C 
399CC 

4 00CC 
4C 1 OC 


: c (KDDIS NE 0) 

SUBROUTINE precis 


;:re 


-UE . NODE 4 . NODE 2 . NODE 3 . 


I 0 I REC 
viwUE 
NODE n 


Qiocal a ,r ec: ’on of o r escr - oea a ’ so > acemer.i 
magmtuae of orescr’ced aiso'acemert 
nodes w ’ ?h tn’s oisd! acement 


SUBROUTINE LOADS 

KDI SP . KSTRES . uODINC , MATCH . NTEMCH . NR £ MOV , NO LOAD . RPM 
KGLOAD . K T LOAD . KALOAO . KL-DAO . KCLCAD . 

METHOD . EPSTOL . TOLEPS , TQlDI S 


4280C 
4290C 
4 300C 
43 IOC 
4320C 
4330C 
4340C 


kd: sp 

KSTRES 

-CD INC 
M£ "■ C 

N 7 £ M 0 H 


prev ious! 


) 


4 D 3 DC 


- 0 4 c c 


4 C50C 
4 0 o w C 

NR E MO v 

*+ w 9 

408DC 
4C90C 
- 4 ODC 

NC LD AC : 

^ i 

RDM • 

4 1 2 DC 

KGLOAD 

- ’30C 
4 140C 
4 1 5CC 

K^lOAD . 

4 160C 
4 17 OC 
4 18 OC 
4 190C 
4 200C 

KALOAD : 

4 2 4 0C 
4220C 
4 2 30C 

KLLOAO : 

4240C 

TOLEPS: 

4 250C 
4260C 
4270C * * 1 

TCLOIS: 


nodal data onnt *’ag 
= 1 crmt nodal data ( a i sp \ acement s i 
as specified in tn^s loaa case 
=0 print nodal data as specified 
=-i no nodal data printing 
element data c-int flag 
=i print e’ement data ' stresses, 
as specified m tr ’s l oaa case 
=0 print element data as specified previously 
* “ ■ n 0 element data printing 
numoer o* lead increments 
it 0 sets I^IM^si (reads t’me aata ) 
material mange *iag 
= 1 tnere * s material Changed 
~ 0 nc material is cnangea 
t ermpe r a tu r e cnange flag 
=i tne^e s mange in temperature 
*C no temperature cnange 
element removal flag 
= r some elements are tc ce remevea 
*C nc elements removed 
concentrated load flag 
s1 tnere are concentrated loaas 
*0 no concentrated 1 oacs 

revolutions per minute fc* tms lead case 
grav 1 1 y i oad . * l ag 
=< there are gravity loads 
=0 no gravity loads 
tnermal load flag 
*i tnermal loads calculated 
s 0 no thermal loads 
area load flag 
»i there are area loads 
*0 no area loads 
line 1 oad f 1 ag 
*1 there are line loads 
*0 no line loads 
convergence tolerance on 
convergence tolerance on 


Strain 

displacement 


I F ( KD ISP . EO . 1 ) 

SUBROUTINE QUTPUTN 

■IDISP. I REF . NO 1 . N02 , N03 , 

C 
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'-E = INPUT PQR MSS8 WITH BODNER JN 


OA~E = * C ' 15/85 


4350C 
4 36 :c 
4 3TCC 
4 38 DC 
439CC 

IDISP 
: r e c 

\0r, n, 

4420C 

4430C 

: r ( *S~RES 

4440C 

SUBROUTINE 

4450C 

4460C 

I El~ . IS’ 

4 4 ?0C 


448CC 
4 4 90C 

I E L ? • 

4 500 C 
45 «0C 
4 520C 

ISTRESS 

453DC 
4540C 
4 5 50C 

ISTRAIN : 

4 560C 
4 5 TQC 



aon t c-irt 
ao print 


caes f; 


wnicn output is aesirea 


458CC 
459 CC 
4 6CCC 
46 iOC 
4 6 2 CC 
4 6 30C 
4640C * 
46 5 DC 


-c3 DC 
46 r DC 


* 3CC 
*4CC 


484CC 

4850C 

4860C 

4870C 

4880C 

4890C 

4 9GOC * 

4 9 1 0C 

4920C 

4930C 

4940C 

4950C 

4960C 

4970C 

4980C 

4990C 

5000C 

50 1 CC • 


:el-.:stress. is^ra in. loca .ieli.ieli.iels. 


force o r ’ nt f l ag 
= i print e'ement ^O'-ces 
=C aon ' t print sprees 
stress print flag 
=• print stresses 
=C aon t print stresses 
strain p»*int flag 
=i print strains 
=0 aon t print strains 
’ ocat ion f o'* wmc^ aata 'S pointed 
s1 center c* gravity of e ; ement 
s 2 element noaes 
=3 Gauss points 

elements for wnic r output aesireo 


I r ( : T IME NE 0 ) ' 

SU5PCU T INE DvNTJM 

n:m . T CRP.T lN:T E r WAXt S IGMAX .ERMA' . CELTMIN . DE-TMUl t 

N2M numDer o * equal t'ne steps m tn»s loaa case 
-C oynanme t 'me mcremert'ng is usea 
wRP total t 1 me in tnis 1 oaa case 
* 0 1 lowing inputs apply fc 


4-50C 

t :n:t’ ■ 

initial 1 1 me s teo 


4 T 8 C C 


=D anq initial load case. 

4 "“’DC 


=C and not initial 

i oaa c 

4T8CC 


previous load case 

479CC 

ECMAX : 

maximum inelastic 

strain 

48C0C 


aef au 1 t * . 000 i 


48 iOC 

sigmax 

maximum cnange in 

stress 

4620C 


aef aul t* 1000 os i 


4 8 30C 

ERMA X ; 

maximum estimated 

Integra 


i me step of 


smed i n any time step 
n er*-or allowed in any time step 

defaui t* .01 ( i%) 

DELTMIN: minimum allowaple t 1 me step 
defaui t* .001 * TCRP 

OELTMULT : maximum .multiple of current time step allowed for next time 
aef au 1 t * i . 5 


I*<match.ne .0) 

SUBROUTINE INPUT2 

MATNC. I EL 1 , IEL2. IEL3. 

0 


MATNO: new material numoer 

IELn: elements changed to this material 
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title 

« INPUT FOR MSS8 wl r ~ SOONER jn 24 T E *10 IS 85 

^ I me = « E 25 

5C20C 



5C3CC 

I f ( NTEMCH NE o ) 


5 0403 

subroutine NOOEIN 


5 05 or 



506GC 

"EMP . NC 1 . NC2 . NC 3 . 


53^03 



5C80C 



5090C 

~Emo ^ew temoe r ature 


5 •■zee 

5 • * OC 

NOr ncaes cnangea to tn, s te-noerat^re 


si:oc 

5130C- 



5 ' -IOC 



5 ■ 5CC 

I'intcmch NE 3) 


5 -50C 

SUBROUTINE NOOEIN 


5*7CC 



5 1 80C 

OTEMP . NO 1 . N02 . NC3 


5 1 90C 

C 


520CC 



52 'OC 

CTEmp temperature gradient 


5220C 
5230C 
5240C • 

NCn nodes w i t r this temoe^a tu-'e graaient 


5 250C 



52SOC 

I r 1 NREMOy NEC) 


52-tOC 

subroutine inpu~‘ 


5280C 



S290C 

I E L 1 . I E L 2 . I E„3 . 


5 300C 

n 



5320C 

5 33CC 

S3J0C 

5 35CC 

53SOC 

53"0C 

5 38 0C 

53SCC 

5400C 

54 # OC 

5 J 20C 

543CC 

5-4QC 

545 CC 

54603' 

5470C 

548GC 

5490C 

5 5 OOC 

5510C- 

5520C 

5530C 

5540C 

5550C 

5560C 

5570C 

5580C 

5590C 

5600C 

56 IOC 

5620C 

5630C 

5640C 

5650C 

5660C* 

5670C 

5680C 


I E L r -'emovea elements 


I c I NCLCAO NE 0 » 

subroutine c^oaos 

1 0 I P EC . VAlUE . N 3 1 . NCI NC 3 


I0IRE3. 9 1 CDa I a^ect’cr c w ’;ac 
vAuJE magnitude of ; oaa 

NOn noaes wn € re tn-s cense" t -a tea loa 3 aedtea 


: F ( K E lO AO NE . 0 ) 
SUBROUTINE eloads 


IF(KALOAD.NE.O) 

SUBROUTINE ALOADS 

I EBEG . I EEND . I E INC . I F ACE .1 2 I RE . 5 1 P2 p 3 P4 
0 


IEBEG: beg i nn i ng element 

IE END: ending element 

IEINC: increment in element numbers 

^ face : element face where pressure load applied 

IDIRE: direction of load in element local cooro. system 

P1.P2.P3.P4: corner pressure values in clockwise? direction 


IF(KLLOAD NE .0) 
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t: t le* 

INPUT FOf 

5S90C 

SUBROUT] 

5'OCC 


5 " • OC 

I E B E G . 

f ** 2 h r 


5 ~ 3CC 


S"40C 

IEEE3 

5 **30C 

: E E ND ■ 

5 ■’6CC 

:e:nc . 

5~ 7 0C 

:-ACE : 

s-soc 

I D I R : 

5 *90C 

Pi , p 2 

5800C 


58 'OC- • 

^ n o.nr . > 


583CC 

5840C 

: * i kpdis 

5850C 

SUBROUTIN 

5860C 


5870C 

I D I REC 

588GC 

C 

569CC 


590CC 

:cirec 

59 iCC 

VAuliE . 

592CC 

NODE n . 

5930C 


5340C 

NOTE 

595CC 


596CC 


597QC - - 



date * ic *5 55 


beginning element 


T — WVW OWL- ’ 

a,rec? '°" = * ,oaa ,n element '.seal cooro s.stem 
Dressures at ettner eno of i me icaa 








D I REC , VALUE . NODE 1 . NODE 2 . NODE 3 . 


g coal airecticn o* oresc^cea aisdacement 
magnitude of oresencea a-soiacement 
nodes with this 3 *. so ’ acement 


***** * = „• » 3 ^LMwC“tNb VUE BE TI 

SAME .«R„u u hout anc MS’S* BE SPECIFIED from T-iE OU 


■iC 

’SET 
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n 


TITLE* INPUT FOR T3CYAN WITH H-A PLA DATE *10/16/85 TIME * 10.27 PAGE* 

1000C ORGANIZATION OF INPUT 

1010C FOR T3CYAN WITH H-A PLASTICITY AND SIMPLE CREEP 

1020C 

1030C I HEADING AND CONTROL INFORMATION 

1040C 1.1 TITLE CARD 

1050C 1.2 PROBLEM SIZING 

1060C 1.3 ANALYSIS AND RESTART OPTIONS 

1070C 1.4 EQUATION NUMBERING ANO BANOING OPTIONS 

1080C 

1090C II NODE COORDINATES ANO TRANSFORMATIONS 
1100C II. 1 NODE COORDINATES 

11 IOC 1 1. 2 LOCAL NODE COORDINATE SYSTEM TRANSFORMATIONS 
1 120C 

1130C III ELEMENT DEFINITION 
1140C III.1 HEADER LINE FOR ELEMENT 

1150C 1 1 1. 2 20 NODED SOLID DEFINITION 

1 160C 

1170C IV LOAD CASE INFORMATION. INITIAL CONDITIONS 
1180C IV. 1 LOAD CASE CONTROL CARD 

1190C IV. 2 ACCELERATION SPECIFICATIONS 

1200C 

12 IOC V MATERIAL PHYSICAL PROPERTIES 

1220C V.1 MATERIAL PHYSICAL PROPERTIES 

1230C V.1.1 ISOTROPIC ELASTIC PROPERTIES 

1240C V.1.2 ORTHOTROPIC ELASTIC PROPERTIES 

1250C V . 2 ORTHOTROPIC AXES ORIENTATION TABLE 
1260C V . 3 INELASTIC MATERIAL CHARACTERIZATION 

1270C 

12B0C VI TIME AND TIME INCREMENTING CONTROL 
1290C 

1300C VII CONVERGENCE CRITERIA 

1 3 IOC 

1320C VIII INITIAL NODAL CONSTRAINED DISPLACEMENTS 
1330C 

1340C IX INITIAL NODAL APPLIED FORCES 
1350C 

1360C X INITIAL NODAL TEMPERATURES 
1370C 

1380C XI INITIAL ELEMENT PRESSURE LOADS 
1390C 

1400C XII LOAD CASE INFORMATION, FINAL CONDITIONS 

14 IOC XII 1 LOAD CASE CONTROL CARD 
1420C XII. 2 ACCELERATION SPECIFICATION 
<430C 

1440C XIII FINAL NODAL CONSTRAINED DISPLACEMENTS 
1450C 

1460C XIV FINAL NODAL APPLIED FORCES 
1470C 

1480C XV FINAL NODAL TEMPERATURES 
1490C 

1500C XVI FINAL ELEMENT PRESSURE LOADS 

1 5 IOC 
1520C 

1530C 

1540C I HEADING AND CONTROL INFORMATION 
1550C 

1560C 1.1 TITLE CARD 

1570C 

1S80C LINE ITITLE 
1590C 

1600C ITITLE • ANY 1 TO 72 CHARACTER TITLE FOR THE ANALYSIS 

16 IOC 

1620C •«*• 

1630C 1.2 PROBLEM SIZING 

1640C 

1650C LINE NUMNP NM IT RTEM NLC 
1660C 



I 


n 

! 


n 

■ i 


1 


n 
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TITLE* INPUT FOR T3CYAN WITH H-A PLA 


DATE -10/16/85 


TIME » 10.27 


0 AGE - 2 


1 6700 
16800 

16 900 
1 7000 

17 IOC 
1720C 
1 7300 
1740C 
17500 

1 7600' 
1 7700 
17800 
1790C 
1800C 
18100 
1 8200 
1830C 
1 8400 
18500 
1860C 
18700 
18800 
18900 
1 9000 

19 IOC 
19200 
19300 
1 9400 
1950C 
1 9600 
19700 
19800 
1 9900 
20000 

20 IOC 
20200 
2030C 
20400* 
20500 
2060C 
2070C 
2080C 
2090C 
21000 

21 IOC 
21200 
21300 
2 1400 
2 1500 
2 1600 
2 1700 
2 1800 
21900 
22000 
2210C- 
22200 
22300 
2240C* 
22500 
22600 

2 2700 
22800 
2 2900 
23000 
23100 
2 3200 
2 3300 


NUMNP 

NM 

IT 


RTEM 

NLC 


NUMBER OF STRUCTURAL NODES ( ENTER AS A NEGATIVE NUMBER FOR 
TIMING SUMMARY ) 

NUMBER OF DIFFERENT MATERIALS ( MAXIMUM-3 ) 

THERMAL STRESS OPTION 
* 0 INCLUDE THERMAL LOADS 
» 1 IGNORE THERMAL LOADS 
REFERENCE TEMPERATURE ( DEGREES F ) 

NUMBER OF LOAD CASES 


1.3 ANALYSIS AND RESTART OPTIONS 
LINE LAWCRP NOUT NRESTA INREST MASSCD 

LAWCRP : TYPE OF INELASTIC ANALYSIS 

* 0 ELASTIC ANALYSIS 

* 1 HAISLER-ALLEN PLASTICITY 

* 2 SECONDARY CREEP MODEL 

*12 PLASTICITY AND CREEP COMBINED 
NOUT : OUTPUT FILE CREATION OPTION 

- 0 DO NOT CREATE OUTPUT FILE 

* 1 CREATE OUTPUT FILE 
NRESTA : RESTART OPTION 

* 0 THIS IS NOT A RESTART RUN 

> 0 INPUT THE LOAD CASE FROM WHICH THE RESTART IS TO PROCEED 
NOTE .(OUTPUT FROM THIS CASE MUST HAVE BEEN PREVIOUSLY PUT ON AN OUTPUT 
FILE. THE FIRST NEW LOAD CASE WILL BE LABELED AS NRESTA ♦ i . ) 

INREST : LOAD CASE NUMBER IN THE CURRENT INPUT FILE WHICH BECOMES 
THE FIRST NEW LOAD CASE TO BE SOLVED WHEN RESTARTING, WHERE 
( 1 .LE. INREST. LE. NLC ). IF 0 IS INPUT INREST - 1 IS ASSUMED 

MASSCO : MASS MATRIX FLAG 

* 0 DO NOT CREATE MASS MATRIX 

* 1 CREATE LUMPED MASS MATRIX 

* 2 CREATE CONSISTENT MASS MATRIX 


1.4 EQUATION NUMBERING AND BANDING OPTIONS 

LINE N IBAND IPBAND 

N : KEY CODE 

* 0 NO NUMBERING OR BANDING 

■ -1' ACTIVATE NUMBERING AND BANDING OPTION 
IBAND : BANDING OPTION 

■ 0 USE DEFAULT OPTION 

- 1 ASSUME NODE NUMBER IS THE SAME AS MATRIX POSITION 
• 2 ASSUME INPUT NODE ORDER DEFINES MATRIX POSITION 
IPBAND : PRINTOUT OPTION 

■ 0 NO PRINTOUT OF EQUATION NUMBERS 

■ 1 PRINT OUT THE EQUATION NUMBER FOR EACH DEGREE OF FREEDOM 


II NODE COORDINATES AND TRANSFORMATIONS 


1 1.1 NODE COORDINATES 
(ENTER THE FOLLOWING LINE FOR EACH NODE) 
LINE N X Y Z 


N 

X 

Y 


NOOE NUMBER 
NODE X COORDINATE 
NODE Y COORDINATE 
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TITLE* INPUT FOR T3CYAN WITH H-A PLA DATE *10/16/85 TIME ■ 10.27 PAGE* 3 

Z : NODE Z COORDINATE 

(AFTER THE LAST NODE ENTER THE FOLLOWING LINE) 

2380C LINE 0 NSKEW 

2390C NSKEW : INDICATOR FOR LOCAL NODE COORDINATE SYSTEM TRANSFORMATION 

2400C * 0 NONE WILL 8£ INPUT 

24 IOC * 1 ONE OR MORE WILL BE INPUT 

2420C 

2430C ............... 

2440C 1 1. 2 LOCAL NODE COORDINATE SYSTEM TRANSFORMATIONS 

2450C 

2460C I F ( NSKEW . GT . 0 ) 

2470C ( ENTER ALL NODES WITH LOCAL COORDINATE SYSTEM TRANSFORMATIONS) 

2480C 

2490C LINE N NO NI NO 
2500C LINE 0 

N : NODE NUMBER HAVING A ROTATED LOCAL COORDINATE SYSTEM 
NO : NODE ON LOCAL X AXIS 

NI : NODE ON LOCAL X AXIS IN THE +X DIRECTION FROM NO 

NJ : NODE IN LOCAL XY PLANE SUCH THAT THE +Z AXIS IS IN THE 

DIRECTION OF THE VECTOR P3 * (NO-NI) X (NO-NJ). THE *Y AXIS IS 

IN THE DIRECTION OF THE VECTOR P2 • P3 X (NO-NI) 


2620C III ELEMENT DEFINITION 
2630C 

2640C**--***--*******«»«*»** 

2650C III.1 HEAOER LINE FOR ELEMENT DATA 
2660C 

2670C LINE IPRINT 

2680C IPRINT ; ELEMENT DATA EXTENDED PRINTOUT OPTION 
2690C * 0 NONE 

2700C - 1 DUMP VOLUMES AND DISTORTION PARAMETERS 

27 IOC * 2 ALSO OUMP ELEMENT STIFFNESSES 

2720C ■ 3 ALSO OUMP ELEMENT EQUIVALENT NODAL FORCES AND MASSES 

2730C 

2740C**"-***--**«*-**-**«-»*****---. 

2750C III. 2 20 NODED SOLID DEFINITION 
2760C 

2770C (ENTER THE FOLLOWING LINES FOR EACH ELEMENT) 

2780C 

2790C LINE NEL NI N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 

2800C LINE N 1 3 N14 N15 N16 N17 N18 N19 N20 IMAT IOR 

2810C LINE NO NP NO ( OPTIONAL LINE. ENTER ONLY IF IOR ■ 100) 

2820C 

2830C LINE 0 
2840C 
2850C 
2B60C 
2870C 
2880C 
2890C 
2900C 
29 IOC 
29200 
2930C 
2940C 
2950C 
2960C 

2970C IV LOAD CASE INFORMATION. INITIAL CONDITIONS 
2980C 

2990C mmmmmmmmmmmmm * mmmmmmmmmmmmmmm + mmmmmmmm ******** 

3000C IV. 1 LOAD CASE CONTROL CARD 


NEL : ELEMENT NUMBER 

NI ... N20 : NODES DEFINING THE ELEMENT ( SEE FIG III. 2 ) 

IMAT : MATERIAL NUMBER 

IOR : ORTHOTROPIC MATERIAL AXIS SYSTEM INDICATOR 
- 0 ISOTROPIC MATERIAL 

(1.LE.I0R.LE.99) IOR IS ORIENTATION IDENTIFIER OF SYSTEM 
GIVEN UNDER MATERIAL DATA SECTION V.2 
* 100 ORTHOTROPIC AXES DEFINED BY OPTIONAL LINE OF INPUT 


25 IOC 
2520C 
2530C 
2540C 
2550C 
2560C 
2570C 
2580C 
2590C 
2600C 


2340C 

2350C 

2360C 

2370C 


\ 1 


ri 




Fi. 


V I 


n 


n 
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TITLE- INPUT FOR T3CYAN WITH H-A PLA 


DATE -10/16/85 


TIME - 10.27 


30 IOC 
3020C 
3030C 
3040C 
3050C 
3060C 
3070C 
3080C 
3090C 
3100C 

31 IOC 
3120C 
3130C 
3140C 
3 150C 
3 160C 
3170C 
3180C 
3190C 
3200C 
321 OC 
3220C 
3230C 
3240C 
3250C 
3260C 


LINE RPM IAXIS IACC NPCI NT I ITHST 


RPM 

IAXIS 


IACC 


NPCI 


NOTE : 
NT I 


ITHST 


: ROTATIONAL SPEED IN RPM 

: GLOBAL AXIS ABOUT WHICH STRUCTURE IS ROTATING 

- 1 X AXIS 

* 2 Y AXIS 
-32 AXIS 

: INDICATOR FOR* INPUT OF ACCELERATON LOADS IN GROUP IV. 2 

* 0 NO INPUT 

* 1 ACCELERATION LOADS ARE INPUT 
INDICATOR FOR A CHANGE IN MATERIAL PROPERTIES 

- 0 NO CHANGE TO MATERIAL PROPERTIES 

- 1 CHANGE ELASTIC PROPERTIES ( GROUPS V.1 AND V 2) 

» 2 CHANGE INELASTIC PROPERTIES ( GROUP V.3 ) 

* 12 CHANGE ELASTIC AND INELASTIC PROPERTIES 

IF THIS IS FIRST LOAD CASE, PROGRAM SETS NPCI-12 

NUMBER OF NODAL TEMPERATURES WHICH ARE RESPECIFIED. 

* 0 NO RESPECIFIED TEMPERATURES . STIFFNESS IS RECOMPUTED 

* -1 NO RESPECIFIED TEMPERATURES. STIFFNESS IS NOT RECOMPUTED 
> 0 NT I RESPECIFIED TEMPERATURES, STIFFNESS IS RECOMPUTED 

: TOTAL STRAIN PRINTOUT OPTION 

» 0 INCLUDE THERMAL STRAINS IN TOTAL STRAIN PRINTOUT 
° 1 00 N0T INCLUDE THERMAL STRAINS IN TOTAL STRAIN PRINTOUT 


3270C- 


IV. 2 ACCELERATION SPECIFICATION FOR INERTIAL OR GRAVITY LOADS 


ACCEL2 


3260C 
3290C 

3300C I F ( IACC.GT.O) 

33 IOC LINE ACCELX ACCELY 
3320C 

3330C ACCELX : ACCELERATION OF STRUCTURE 

3340C ACCELY : ACCELERATION OF STRUCTURE 

3350C ACCELZ : ACCELERATION OF STRUCTURE 

3360C 
3370C 


IN/SEC* *2 
IN/SEC--2 
IN/SEC* -2 


IN GLOBAL 
IN GLOBAL 
IN GLOBAL 


V MATERIAL PHYSICAL PROPERTIES 


V.1 ELASTIC CONSTANTS 

I F ( NPC I . EO. 1 .OR. NPCI. EO. 12) 

(ENTER NM OF THE FOLLOWING LINES) 


MTN NMT DEN 


) 


3380C- 
3390C 
3400C 

34 IOC* 

3420C 
3430C 
3440C 
3450C 
3460C 

3470C LINE 
3480C 

3490C MTN 
3500C 

35 IOC 

3520C NMT 
3530C 

3540C DEN 
3550C 

3 b6OC*#00#00 000*000*0##** 0000 *0*0 00000000 0 00*0*00*00000000000 

3570C 
3580C 

35900 
3600C 

36 IOC 

3620C LINE 
3630C 
3640C 
3650C 
3660C 
3670C 


MATERIAL NUMBER ( MTN . GE . 1 . AND . MTN . LE . NM 
(IF MATERIAL IS ISOTROPIC, INPUT MTN AS A NEGATIVE NUMBER 
TO SIMPLIFY INPUT) 

NUMBER OF TEMPERATURES AT WHICH ELASTIC PROPERTIES WILL BE 
GIVEN FOR THIS MATERIAL 

WEIGHT DENSITY OF THE MATERIAL ( P0UNDS/IN--3 ) 


V.1.1 ISOTROPIC MATERIAL 


I F ( MTN . LT .0) 

( ENTER NMT OF THESE LINES 


TEMP E PR AL 


TEMP 

E 

PR 

AL 


DEGREES F. ) 

10**6 P.S.I 


TEMPERATURE ( 

ELASTIC MODULUS ( 

POISSON'S RATIO 

MEAN COEFFICIENT OF THERMAL EXPANSION ( 10**-6 IN/IN-DEG . F. 


) 
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H 

; 1 


DATE * 10/ 16/85 


TIME * 10.27 


TITLE* INPUT FOR T3CYAN WITH H-A PLA 
3S80C 




PAGE- 5 


3700C 

37 IOC 
3720C 
3730C 
3740C 
3 750C 
3760C 
3770C 
3780C 
3790C 
3800C 

38 IOC 
3820C 
3830C 
3840C 
3850C 
3860C 
3870C 
3880C 
3890C 
3900C 

39 IOC 


V.1.2 ORTHOTROPIC MATERIAL 
IF(MTN.GT .0) 

(ENTER NMT OF THE FOLLOWING LINES) 

(DIRECTIONS 1.2,3 CORRESPOND TO THE MATERIAL ORTHOTROPIC AXES 
LINE TEMP Ell E22 E33 NU12 NU13 NU23 G12 G23 G3 1 AL1 AL2 AL3 


. Y' ,2' ) 


TEMP 

Ell 

E22 

E33 

NU12 

NU 13 

NU23 

G 12 

G23 

G31 

AL 1 

AL2 

AL3 


TEMPERATURE AT WHICH PROPERTIES ARE GIVEN 
ELASTIC MODULUS IN THE 1 DIRECTION 
ELASTIC MODULUS IN THE 2 DIRECTION 
ELASTIC MODULUS IN THE 3 DIRECTION 
POISSON'S RATIO RELATING DIRECTIONS 
POISSON'S RATIO RELATING DIRECTIONS 
POISSON'S RATIO RELATING DIRECTIONS 
SHEAR MODULUS IN THE 1-2 PLANE 
SHEAR MODULUS IN THE 2-3 PLANE 
SHEAR MODULUS IN THE 3-1 PLANE 
MEAN COEFFICIENT OF THERMAL EXPANSION IN THE 

MEAN COEFFICIENT OF THERMAL EXPANSION IN THE 

MEAN COEFFICIENT OF THERMAL EXPANSION IN THE 


AND 

AND 

AND 


( DEG. F. ) 


DIRECTION 

DIRECTION 


3 DIRECTION 




3920C- 

3930C 

3940C 

3950C 

3960C 

3970C 

3980C 

3990C 

4000C 

40 IOC 

4020C 

4030C 

4040C 

4050C 

4060C 

4070C 

4080C 

4090C 

4 100C 

4 11 OC * 

4 120C 

4130C 

4140C 

4 150C 

4160C 

4 170C 

4 180C 

4190C 

420OC 

42 IOC 
4220C 
4230C 
4240C 
4250C 
4260C 
4270C 
4280C 
429DC 
4300C 

43 IOC 
4320C 
4330C 
4340C 


V.2 ORTHOTROPIC AXIS ORIENTATION TABLE 
LINE NOR 

NOR : NUMBER OF ORIENTATION SPECIFICATIONS 

(O.LE.NOR.LE. 10) i.e. MAX OF 10 ORIENTATION SYSTEMS ALLOWED 
(NOTE: IF THE MATERIALS ARE ISOTROPIC OR THE ORTHOTROPIC AXES 
COINCIDE WITH THE GLOBAL AXES. ENTER 0) 

ENTER O. ( MAXIMUM OF 10 SPECIFICATIONS ALLOWED ) 

I F ( NOR . GT . 0 ) 

( ENTER NOR LINES OF THE FOLLOWING ) 

LINE I NO NP NO 

I : ORIENTATION IDENTIFIER ( IOR IN ELEMENT INPUT ) 

NO, NP, NO : NODE NUMBERS IDENTIFYING ORTHOTROPIC AXES ( SEC. II. 2 ) 

V.3 INELASTIC MATERIAL PROPERTIES 

(IF NPCI . EO . 2 . OR . NPCI . EO - 12 ) 

(IF LA WCRP . EO . 1 . OR . LAWCRP . EO . 12) 

(ENTER MTN OF THE NEXT FOUR LINE SETS - MTN , SSTEMP , PPV , AND BET) 
LINE MTN NPTS NTM 


MTN 

NPTS 

NTM 


MATERIAL NUMBER 

NUMBER OF STRESS-STRAIN POINTS PER CURVE 

NUMBER OF TEMPERATURES FOR WHICH STRESS-STRAIN CURVES ARE GIVEN 


LINE SSTEMP(I) SSTEMP(2) 


SSTEMP(NTM) 


SSTEMP : ARRAY OF TEMPERATURES WHERE STRESS-STRAIN CURVES ARE GIVEN 
(IN INCREASING ORDER) 

LINE PPV ( MTN ,1,1,1) PPV(MTN. I . 1,2) 

LINE PPV (MTN, 1,2,1) PPV< MTN , I , 2 , 2 ) 
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TITLE* INPUT FOR T3CYAN WITH H-A PLA DATE *10/16/85 TIME * 10.27 PAGE* 

4350C LINE PPV(MTN. I .NPTS. 1 ) PPV(MTN, I .NPTS.2) 

4360C 

4370C PPVtMTN, I . J, 1 ) : STRESS VALUES IN INCREASING ORDER FROM J*1 TO NPTS 

4380C FOR TEMPERATURE I ANO MATERIAL MTN 

4 390C PPV( MTN . I . <J , 2 ) : STRAIN VALUES IN INCREASING ORDER FROM 0*1 TO NPTS 

4400C FOR TEMPERATURE I AND MATERIAL MTN 

44 IOC 

4420C LINE BET ( 1 ) BET(2) BET(NTM) 

4430C 

4440C BET : ARRAY OF HARDENING COEFFICIENTS FOR SAME TEMPERATURES WHERE 

4450C STRESS-STRAIN CURVES GIVEN 

4460C 

4470C (IF L AWCRP . EO . 2 . OR . LAWCRP . EO . 1 2 ) 

4480C (ENTER MTN OF THE FOLLOWING TWO LINE SETS) 

4490C 

4500C LINE NCTP TOUT SNORM 
451 OC 

4520C NCTP : NUMBER OF TEMPERATURES WHERE CREEP PROPERTIES SPECIFIED 
4530C TCUT : CUTOFF TEMPERATURE BELOW WHICH NO CREEP OCCURS 

4540C SNORM : NORMALIZING STRESS (IN PSI) FOR THE FOLLOWING CREEP PROPERTIES 

4S50C 

4560C (ENTER NCTP OF THE FOLLOWING LINE) 

4570C 

4580C LINE TEMP 0 R STRCUT 
4590C 

4600C TEMP : TEMPERATURE WHERE THESE CREEP PROPERTIES APPLY 

46 IOC 0 : CREEP PROPERTY 

4620C R : CREEP PROPERTY (CREEP STRAIN « 0 • STRESS •• R) 

4630C STRCUT : CUTOFF STRESS (IN PSD BELOW WHICH NO CREEP OCCURS 

4640C 

4650C 

4660C 

4670C VI TIME AND TIME INCREMENTING CONTROL INPUT 
4680C 

4690C LINE N2M TCRP TINIT ECMAX SIGMAX ERMAX DELMIN DELMUL 
4700C 

4 7 IOC N2M ; NUMBER OF EOUAL TIME STEPS IN THIS LOAD CASE, IF N2M*0 

4720C DYNAMIC TIME INCREMENTING WILL BE USED 

4730C TCRP : TOTAL TIME IN THIS LOAD CASE ( SECONDS ) . 

4740C TINIT : INITIAL TIME STEP. IF THIS IS NOT THE FIRST LOAD CASE AND 

4750C 0.0 IS INPUT, .5 TIMES THE LAST CALCULATED TIME STEP OF THE 

4760C PREVIOUS LOAD CASE IS USED. IF THIS IS THE FIRST LOAD CASE OR 

4770C IF THIS IS THE FIRST LOAD CASE OF A RESTART AND 0.0 IS INPUT, 

4780C DELMIN IS USED. 

4790C ECMAX : MAXIMUM INELASTIC STRAIN INCREMENT DESIRED IN ANY TIME STEP. 

4800C DEFAULT VALUE IS .000100. 

48 IOC SIGMAX : MAXIMUM CHANGE IN STRESS OESIRED IN ANY TIME STEP. 

4820C DEFAULT VALUE IS 1000 PSI. 

4830C ERMAX : MAXIMUM ESTIMATED INTEGRATION ERROR ALLOWED IN ANY TIME STEP. 

4840C DEFAULT IS .01 ( 1% ). 

4850C DELMIN : MINIMUM ALLOWABLE TIME STEP. DEFAULT VALUE IS .001 TIMES 
4860C TCRP. 

4870C DELMUL : MAXIMUM MULTIPLIER ALLOWED ON CURRENT TIME STEP IN 

4880C CALCULATING THE NEXT TIME STEP. DEFAULT • 1.5 

4890C 

4900C 

49 IOC 

4920C VII CONVERGENCE CRITERIA 


4930C 
4940C 
4950C 
4960C 
4970C 
4980C 
4990C 
50OOC 
50 IOC 


LINE I DM ICON DELSIG DELEPS 

I DM : MATERIAL NUMBER 

ICON : CONVERGENCE CRITERIA CODE 

■ 1 CHANGE IN EFFECTIVE STRESS FOR SUBSEOUENT ITERATIONS 
MUST BE LESS THAN DELSIG 

- 2 CHANGE IN EFFECTIVE INELASTIC STRAIN INCREMENT FOR 

SUBSEOUENT ITERATIONS MUST BE LESS THAN DELEPS DELEPS 
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TITLE* INPUT FOR T3CYAN WITH H-A PLA DATE *10/16/85 TIME « 10.27 PAGE* 

5020C * 3 BOTH CONDITIONS 1 AND 2 MUST BE SATISFIED. DEFAULT • 3 

5030C DELSIG : CONVERGENCE TOLERANCE ON EFFECTIVE STRESS 
5040C DEFAULT ■ ,01*SIGMAX 

5050C OELEPS : CONVERGENCE TOLERANCE ON EFFECTIVE INELASTIC STRAIN INCREMENT. 

S060C DEFAULT * . 01 *ECMAX 

5070C 

5080C 

S090C 

5100C VIII INITIAL CONSTRAINED DISPLACEMENTS 

51 IOC 

5 1 20C LINE NOB 
5130C 

5140C NDB * NUMBER OF CONSTRAINED DISPLACEMENT SPECIFICATIONS 
5 1 50C 

51S0C ( ENTER NDB OF THE FOLLOWING LINES ) 

5170C 

5180C LINE N IDIR VALUE NEND NINC 
5190C 

5200C N : NODE NUMBER 

52 IOC IDIR : DIRECTION CONSTRAINED 

5220C * 1 X CONSTRAINED ( ANY COMBINATION OF CODES 

5230C * 2 Y CONSTRAINED i.«. 12. 13. OR 123 

5240C *32 CONSTRAINED MAY ALSO BE USED) 

5250C VALUE : NUMERICAL VALUE OF CONSTRAINED DISPLACEMENT. DEFAULT IS 0.0 

5260C NEND : LAST NODE NUMBER HAVING THIS CONSTRAINT. IF OMITTED N IS ASSUMED 

5270C NINC : INCREMENT TO BE USED FOR CONSTRAINT GENERATION FROM N TO NEND 

5280C DEFAULT IS 1 . ( OPTIONAL ) 

5290C 

5300C 

53 IOC 

5320C IX INITIAL NODAL APPLIED FORCES 

5330C 

5340C LINE NFB 
5350C 

5360C NFB : NUMBER OF APPLIED NODAL FORCE SPECIFICATIONS 
5370C 

5380C ( ENTER NFB OF THE FOLLOWING LINES ) 

5390C 

S40OC LINE N IDIR VALUE NENO NINC 

54 IOC 

5420C N : NODE NUMBER 

5430C IDIR : DIRECTION OF APPLIED FORCE 

5440C * 1 X DIRECTION 

5450C - 2 Y DIRECTION 

5460C • 3 Z DIRECTION 

5470C VALUE : NUMERICAL VALUE OF FORCE TO BE APPLIED ( LBS ) 

5480C NEND : LAST NODE HAVING THIS APPLIED FORCE. IF OMITTED N IS ASSUMED 

5490C NINC : INCREMENT TO BE USED FOR FORCE GENERATION FROM N TO NEND. 

5500C DEFAULT IS 1 . ( OPTIONAL ) 

55 IOC 
5520C 

5530C ... 

5540C X INITIAL NODAL TEMPERATURES 
5S50C 

5560C LINE NTEMPS 
5570C 

5580C NTEMPS : NUMBER OF TEMPERATURE INPUT LINES 
5S90C 

5600C ( ENTER NTEMPS OF THE FOLLOWING LINES ) 

56 IOC 

5620C LINE N TEMP I TYPE NENO NINC 
5630C 

5640C N : NODE NUMBER 

5650C TEMP : TEMPERATURE ( DEGREES F. ) 

5660C ITYPE : INDICATOR FOR TYPE OF TEMPERATURE REVISION 

S670C - 0 CHANGE NODAL TEMPERATURE TO TEMP 

5680C ■ 1 INCREMENT NODAL TEMPERATURE BY TEMP 
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TITLE- INPUT FOR T3CYAN WITH H-A PLA DATE -10/16/85 TIME » 10 27 PAGE- 8 

r- 

5690C NEND : LAST NODE HAVING THE GIVEN TEMPERATURE. IF OMITTED N IS ASSUMED 

5700C NINC : INCREMENT TO BE USED ON TEMPERATURE GENERATION FROM N TO NEND. 

57 IOC IF OMITTED 1 IS ASSUMED 

5720C 

f 5730C 

| 

S750C XI INITIAL ELEMENT PRESSURE LOAOS 
5760C 

5770C LINE NPL 
5780C 

! ; 5790C NPL : NUMBER OF PRESSURE INPUT LINES 

5800C 

58 IOC ( ENTER NPL OF THE FOLLOWING LINES ) 

r~ 5820C 

5830C LINE NEL I FACE PI P2 P3 P4 P5 P6 P7 P8 NEND NINC 
5840C 

5850C NEL : ELEMENT NUMBER 

5860C IFACE : FACE NUMBER ( I FACE . GE . 1 . AND . I FACE . LE . 6 ) 
f ' 5870C PI : PRESSURE AT NODE 1 ( PSI ) 

i 5880C P2 . P3 , P4, PS, P6 , P7 . P8 : PRESSURES AT NODES 2 8 ON THE FACE. 

5890C IF OMITTED THEY ARE SET EOUAL TO Pi. 

5900 (POSITIVE PRESSURES INDUCE COMPRESSION IN THE ELEMENT ) 

r- 59 IOC NEND : LAST ELEMENT HAVING THIS PRESSURE LOADING. IF OMITTED NEL IS 

E 5920C ASSUMED 

5930C NINC : INCREMENT TO BE USED FOR ELEMENT PRESSURE GENERATION FROM 

5940C NEL TO NEND. IF OMITTED NEL IS ASSUMED 

S950C 

r 5960C 

( 5970C- - 

5980C XII LOAD CASE INFORMATION. FINAL CONDITIONS 
5990C 

6000C-- 

f 60 IOC XII. 1 LOAD CASE CONTROL CARD 

i ) 6020C 

6030C LINE RPM IAXIS IACC NT I 
6040C 

6050C - 

'| 6060C XII. 2 ACCELERATION SPECIFICATIONS FOR INERTIA OR GRAVITY LOADS 

j 6070C 

6080C LINE ACCELX ACCELY ACCEL2 
6090C 

r" 6iooc 

J 61 IOC 

! 6120C XIII NODAL CONSTRAINED DISPLACEMENTS 

6130C 

^ 6140C LINE NOB 

f 6 1 50C 

! 6160C ( ENTER NDB OF THE FOLLOWING LINES ) 

1 6 1 70C 

6180C LINE N IDIR VALUE NEND NINC 
r 6 1 90C 

i 620OC 

I 62 IOC - - 

6220C XIV NODAL APPLIED FORCES 
6230C 

r 6240C LINE NFB 

j 6250C 

6260C ( ENTER NFB OF THE FOLLOWING LINES ) 

6270C 

6280C LINE N IDIR VALUE NEND NINC 

6290C 

6300C 

63 IOC - 

6320C XV NODAL TEMPERATURES 
* — eionr 

6340C LINE NTEMPS 
63S0C 
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TITLE* INPUT FOR T3CYAN WITH H-A PLA DATE *10/16/85 TIME • 10.27 PAGE* 9 

6360C ( ENTER NTEMPS OF THE FOLLOWING LINES ) 

6370C 1 

6380C LINE N TEMP ITYPE NEND NINC 
6390C 
6400C 

64 IOC 

6420C XVI ELEMENT PRESSURE LOADS 

6430C 

6440C LINE NPL 
6450C 

6460C ( ENTER NPL OF THE FOLLOWING LINES ) 

6470C 

6480C LINE NEL IFACE Pi P2 P3 P4 P5 P6 p7 pg NEND NINC 


n 



i 






' ! 
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TITLE* INPUT FOP T3CVAN WITH SIMPLE DATE *10/16/85 TIME * 10.28 

lOOOC ORGANIZATION OF INPUT 

1 0 IOC FOR T3CYAN WITH SIMPLE PLASTICITY AND SIMPLE CREEP 

1020C 

1030C I HEADING AND CONTROL INFORMATION 

1040C I . 1 TITLE CARD 

1050C 1.2 PROBLEM SIZING 

1060C 1.3 ANALYSIS AND RESTART OPTIONS 

1070C 1.4 EQUATION NUMBERING AND BANDING OPTIONS 

1080C 

1090C II NODE COORDINATES AND TRANSFORMATIONS 
1100C II. 1 NODE COORDINATES 

11 IOC 1 1. 2 LOCAL NODE COORDINATE SYSTEM TRANSFORMATIONS 
1 1 20C 

1130C III ELEMENT DEFINITION 
1140C III.1 HEADER LINE FOR ELEMENT 
1150C III. 2 20 NODED SOLID DEFINITION 

1 160C 

1170C IV LOAD CASE INFORMATION, INITIAL CONDITIONS 
1180C IV. 1 LOAD CASE CONTROL CARD 

1190C IV. 2 ACCELERATION SPECIFICATIONS 

1200C 

1 2 IOC V MATERIAL PHYSICAL PROPERTIES 

1220C V.1 MATERIAL PHYSICAL PROPERTIES 

1230C V.1.1 ISOTROPIC ELASTIC PROPERTIES 

1240C V.1.2 ORTHOTROPIC ELASTIC PROPERTIES 

1250C V , 2 ORTHOTROPIC AXES ORIENTATION TABLE 
1260C V . 3 INELASTIC MATERIAL CHARACTERIZATION 

1270C 

1280C VI TIME AND TIME INCREMENTING CONTROL 
1290C 

1300C VII CONVERGENCE CRITERIA 

1 3 IOC 

1320C VIII INITIAL NODAL CONSTRAINED DISPLACEMENTS 
1330C 

1340C IX INITIAL NODAL APPLIED FORCES 
1350C 

1360C X INITIAL NODAL TEMPERATURES 
1370C 

1380C XI INITIAL ELEMENT PRESSURE LOADS 
1390C 

1400C XII LOAD CASE INFORMATION. FINAL CONDITIONS 

14 IOC XII. 1 LOAD CASE CONTROL CARD 
1420C XII. 2 ACCELERATION SPECIFICATION 
1430C 

1440C XIII FINAL NODAL CONSTRAINED DISPLACEMENTS 
1450C 

1460C XIV FINAL NODAL APPLIED FORCES 
1470C 

1480C XV FINAL NODAL TEMPERATURES 
1490C 

150OC XVI FINAL ELEMENT PRESSURE LOADS 

1 5 IOC 

1520C 

1530C - 

1540C I HEADING AND CONTROL INFORMATION 
1550C 

1560C I . 1 TITLE CARD 

1570C 

1580C LINE ITITLE 
1590C 

1600C ITITLE - ANY 1 TO 72 CHARACTER TITLE FOR THE ANALYSIS 
16 IOC 

1620C****** 

1630C 1.2 PROBLEM SIZING 

1640C 

1650C LINE NUMNP NM IT RTEM NLC 
1660C 


— ge/mhinep 


PAGE* 1 


PAGE- i 


145 



TITLE* INPUT FOR T3CYAN WITH SIMPLE 


DATE *10/16/86 


TIME « 10.28 


PAGE- 


16700 

1680C 

16900 

1 7000 

1710C 

1 7200 

1 7300 

17400 

1 7500 

17600* 

1 7700 

17800 

1790C 

1800C 

18 IOC 
1 8200 
1 8300 
18400 
1850C 
1860C 
18700 
1 8800 
1 8900 
19000 

1 9 IOC 
19200 
1 9300 
19400 
1 9500 
19600 
1 9700 
1980C 
19900 
20000 

20 IOC 
2020C 
20300 
20400* 
2050C 
2060C 
2070C 
20800 
2090C 
21000 

2 1 100 
2 1200 
2130C 
2140C 
2 150C 
2 1600 
2170C 
2 1800 
21900 
22000 
2210C- 
2220C 
22300 


NUMNP 

NM 

IT 


RTEM 

NLC 


NUMBER OF STRUCTURAL NODES 
TIMING SUMMARY ) 

NUMBER OF DIFFERENT MATERIALS 
THERMAL STRESS OPTION 

• 0 INCLUDE THERMAL LOADS 

* 1 IGNORE THERMAL LOADS 
REFERENCE TEMPERATURE ( DEGREES 
NUMBER OF LOAD CASES 


( ENTER AS A NEGATIVE NUMBER FOR 


( MAX IMUM*3 ) 


F ) 


1.3 ANALYSIS AND RESTART OPTIONS 
LINE LAWCRP NOUT NRESTA INREST MASSCD 

LAWCRP : TYPE OF INELASTIC ANALYSIS 

* 0 ELASTIC ANALYSIS 

» 1 HAISLER-ALLEN PLASTICITY 

- 2 SECONDARY CREEP MODEL 

*12 PLASTICITY AND CREEP COMBINED 
NOUT ; OUTPUT FILE CREATION OPTION 

* 0 DO NOT CREATE OUTPUT FILE 

■ 1 CREATE OUTPUT FILE 
NRESTA ; RESTART OPTION 

* 0 THIS IS NOT A RESTART RUN 

* 2 INPUT THE L0AD CASE FR0M * HICH THE RESTART IS TO PROCEED 
NOTE: (OUTPUT FROM THIS CASE MUST HAVE BEEN PREVIOUSLY PUT ON AN OUTPUT 

FILE. THE FIRST NEW LOAD CASE WILL BE LABELED AS NRESTA +1.) 

INREST : LOAD CASE NUMBER IN THE CURRENT INPUT FILE WHICH BECOMES 
THE FIRST NEW LOAD CASE TO BE SOLVED WHEN RESTARTING , WHERE 
( 1 . LE . INREST. LE .NLC ). IF 0 IS INPUT INREST ■ 1 IS ASSUMED 

MASSCD : MASS MATRIX FLAG 

* 0 DO NOV CREATE MASS MATRIX 

- 1 CREATE LUMPED MASS MATRIX 

■ 2 CREATE CONSISTENT MASS MATRIX 


1.4 EQUATION NUMBERING AND BANDING OPTIONS 
LINE N IBAND IPBAND 


N 


IBAND 


IPBAND 


KEY CODE 

* 0 NO NUMBERING OR BANDING 

■ -1 ACTIVATE NUMBERING AND BANDING OPTION 
BANDING OPTION 

• 0 USE OEFAULT OPTION 

1 ASSUME NOOE NUMBER IS THE SAME AS MATRIX POSITION 

2 ASSUME INPUT NODE OROER DEFINES MATRIX POSITION 
PRINTOUT OPTION 

0 NO PRINTOUT OF EOUATION NUMBERS 

1 PRINT OUT THE EOUATION NUMBER FOR EACH DEGREE OF FREEDOM 


II NOOE COORDINATES AND TRANSFORMATIONS 


22400" 




2250C 

22600 

II 

1 NODE COORDINATES 

2270C 

22800 

(ENTER 

THE FOLLOWING 

22900 
2 3000 

LINE 

N X 

Y 2 

23100 

N : 

NOOE 

NUMBER 

23200 

X : 

NOOE 

X COORDINATE 

2330C 

Y ; 

NODE 

Y COORDINATE 


GE/MHINEP 
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2340C Z : NODE 2 COORDINATE 
2350C 

2360C (AFTER THE LAST NOOE ENTER THE FOLLOWING LINE) 

2370C 

2380C LINE 0 NSKEW 

2390C NSKEW : INDICATOR FOR LOCAL NODE COORDINATE SYSTEM TRANSFORMATION 

2400C * 0 NONE WILL BE INPUT 

24 IOC * 1 ONE OR MORE WILL BE INPUT 


2420C 

2430C' 

2440C 

24 50C 

2460C 

2470C 

2480C 

2490C 

2500C 


2510C 


2520C 

N 

2530C 

NO 

2540C 

NI 

2550C 

NJ 

2560C 


2570C 


2580C 


2590C 


2600C 


2610C- 



2620C 

III 

2630C 
2640C* 1 


2650C 

I] 

2660C 



2670C 
2680C 
2690C 
2700C 
27 IOC 
2720C 
27 30C 
2740C" 


II -2 LOCAL NODE COORDINATE SYSTEM TRANSFORMATIONS 
I F ( NSKEW . GT . 0 ) 

( ENTER ALL NODES WITH LOCAL COORDINATE SYSTEM TRANSFORMATIONS) 

LINE N NO NI NJ 
LINE 0 

NODE NUMBER HAVING A ROTATED LOCAL COORDINATE SYSTEM 
NODE ON LOCAL X AXIS 

NODE ON LOCAL X AXIS IN THE +X DIRECTION FROM NO 
NODE IN LOCAL XY PLANE SUCH THAT THE ♦ 1 AXIS IS IN THE 
DIRECTION OF THE VECTOR P3 - (NO-NI) X (NO-NJ). THE +Y AXIS IS 
IN THE DIRECTION OF THE VECTOR P2 » P3 X (NO-NI) 


III ELEMENT DEFINITION 


IIX.1 HEADER LINE FOR ELEMENT DATA 
LINE IPRINT 

IPRINT : ELEMENT DATA EXTENDED PRINTOUT OPTION 

* 0 NONE 

* 1 DUMP VOLUMES AND DISTORTION PARAMETERS 
■ 2 ALSO DUMP ELEMENT STIFFNESSES 

* 3 ALSO DUMP ELEMENT EQUIVALENT NODAL FORCES AND MASSES 


2750C III. 2 20 NODED SOLID DEFINITION 
2760C 

2770C (ENTER THE FOLLOWING LINES FOR EACH ELEMENT) 

2780C 

2790C LINE NEL NI N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 

2800C LINE N 1 3 N14 N15 N16 N17 N18 N19 N20 IMAT IOR 

2810C LINE NO NP NO ( OPTIONAL LINE, ENTER ONLY IF IOR - 100) 
2820C 

2B30C LINE 0 
2840C 

2850C NEL : ELEMENT NUMBER 

2860C NI ... N20 : NODES DEFINING THE ELEMENT ( SEE FIG III. 2 ) 
2870C IMAT : MATERIAL NUMBER 

2880C IOR : ORTHOTROPIC MATERIAL AXIS SYSTEM INDICATOR 
2890C - 0 ISOTROPIC MATERIAL 

2900C ( 1 .LE.I0R.LE.99) IOR IS ORIENTATION IDENTIFIER OF SYSTEM 

2910C GIVEN UNDER MATERIAL DATA SECTION V.2 

2920C - 100 ORTHOTROPIC AXES DEFINED BY OPTIONAL LINE OF INPUT 

2930C 

2940C 

2950C 

2960C - - 

2970C IV LOAD CASE INFORMATION, INITIAL CONDITIONS 
2980C 



3000C IV. 1 LOAD CASE CONTROL CARD 
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TITLE- INPUT FOR T3CYAN WITH SIMPLE 


DATE *10/16/85 


30 IOC 
3020C 
3030C 
3040C 
3050C 
3060C 
3070C 
3080C 
3090C 
3100C 

31 IOC 
3120C 
3130C 
3 140C 
3 150C 
3160C 
3 170C 
3 180C 
3 190C 
3200C 

32 IOC 
3220C 
3230C 
3240C 
3250C 
3260C 


LINE RPM IAXIS IACC NPCI NT I ITHST 


RPM 

IAXIS 


ITHST 


; ROTATIONAL SPEED IN RPM 

; GLOBAL AXIS ABOUT WHICH STRUCTURE IS ROTATING 

* 1 X AXIS 

* 2 Y AXIS 

- 3 Z AXIS 

: INDICATOR FOR INPUT OF ACCEtERATON LOADS IN GROUP IV. 2 

- 0 NO INPUT 

■ 1 ACCELERATION LOAOS ARE INPUT 
INDICATOR FOR A CHANGE IN MATERIAL PROPERTIES 
1 0 NO CHANGE TO MATERIAL PROPERTIES 

* 1 CHANGE ELASTIC PROPERTIES ( GROUPS V.1 AND V 2) 

* 2 CHANGE INELASTIC PROPERTIES ( GROUP V.3 ) 

* 12 CHANGE ELASTIC AND INELASTIC PROPERTIES 

IF THIS IS FIRST LOAD CASE, PROGRAM SETS NPCI-12 

NUMBER OF NODAL TEMPERATURES WHICH ARE RESPECIFIED, 

- 0 NO RESPECIFIED TEMPERATURES, STIFFNESS IS RECOMPUTED 

- - 1 NO RESPECIFIED TEMPERATURES, STIFFNESS IS NOT RECOMPUTED 
> 0 NT I RESPECIFIED TEMPERATURES, STIFFNESS IS RECOMPUTED 

: TOTAL STRAIN PRINTOUT OPTION 

■ 0 INCLUDE THERMAL STRAINS IN TOTAL STRAIN PRINTOUT 

- 1 DO NOT INCLUDE THERMAL STRAINS IN TOTAL STRAIN PRINTOUT 


3260C IV. 2 ACCELERATION SPECIFICATION FOR INERTIAL OR GRAVITY LOADS 
3290C 

3300C IF(IACC.GT.O) 

33 IOC LINE ACCELX ACCELY ACCELZ 
3320C 

3330C ACCELX : ACCELERATION OF STRUCTURE ( IN/SEC*»2 ) IN GLOBAL X 

3340C ACCELY : ACCELERATION OF STRUCTURE ( IN/SEC**2 ) IN GL08AL Y 

3350C ACCELZ : ACCELERATION OF STRUCTURE ( IN/SEC--2 ) IN GLOBAL Z 


3340C ACCELY 
3350C ACCELZ 
3360C 
3370C 


3380C- 


3390C V MATERIAL PHYSICAL PROPERTIES 
3400C 


3420C V.1 ELASTIC CONSTANTS 

3430C 

3440C IF (NPCI . EO. 1 .OR. NPCI . EO . 12) 

3450C (ENTER NM OF THE FOLLOWING LINES) 

3460C 

3470C LINE MTN NMT DEN 
3480C 

3490C MTN : MATERIAL NUMBER < MTN . GE . 1 . AND . MTN . LE . NM ) 

3500C (IF MATERIAL IS ISOTROPIC. INPUT MTN AS A NEGATIVE NUMBER 

3510C TO SIMPLIFY INPUT) 

3520C NMT : NUMBER OF TEMPERATURES AT WHICH ELASTIC PROPERTIES WILL BE 
3530C GIVEN FOR THIS MATERIAL 

3540C DEN : WEIGHT DENSITY OF THE MATERIAL ( PQUNDS/IN-*3 ) 

3550C 

3570C V.1.1 ISOTROPIC MATERIAL 


IF(MTN.LT.O) 

( ENTER NMT OF THESE LINES ) 
(NE TEMP E PR AL 


3640C 

3650C 

3660C 

3670C 


TEMPERATURE ( DEGREES F. ) 

ELASTIC MOOULUS ( 10--6 P.S.I. ) 

POISSON'S RATIO 

MEAN COEFFICIENT OF THERMAL EXPANSION ( 10**-6 IN/IN-OEG. 
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TITLE- INPUT FOR T3CYAN WITH SIMPLE 


DATE -10/16/85 


TIME * 10.28 


PAGE « 


3680C 

36 90C* * 0 # * * #0 * 0 * M # 0 0 H 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 * 0 * 0 0 0 * 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 


3700C 

37 IOC 
3720C 
3730C 
3740C 
3750C 
3760C 
3770C 
3780C 
3790C 
3800C 

38 IOC 
3820C 
3830C 
3840C 
3850C 
3860C 
3870C 
3880C 
3890C 
3900C 

39 IOC 


V.1.2 ORTHOTROPIC MATERIAL 
IF(MTN.GT.O) 

(ENTER NMT OP THE FOLLOWING LINES) 

(DIRECTIONS 1.2,3 CORRESPOND TO THE MATERIAL ORTHOTROPIC AXES X'.Y 
LINE TEMP Ell E22 E33 NU12 NU13 NU23 G12 G23 G3 1 AL1 AL2 AL3 


.Z' ) 


TEMP 

Ell 

E22 

£33 

NU 1 2 

NU 1 3 

NU23 

G 1 2 

G23 

G31 

AL1 

A L2 

AL3 


TEMPERATURE AT WHICH PROPERTIES ARE GIVEN ( DEG. F. ) 

ELASTIC MODULUS IN THE 1 DIRECTION 

ELASTIC MODULUS IN THE 2 DIRECTION 

ELASTIC MODULUS IN THE 3 DIRECTION 

POISSON ' S RATIO RELATING DIRECTIONS 1 AND 2 

POISSON'S RATIO RELATING DIRECTIONS 1 AND 3 

POISSON'S RATIO RELATING DIRECTIONS 2 AND 3 

SHEAR MODULUS IN THE 1-2 PLANE 

SHEAR MODULUS IN THE 2-3 PLANE 

SHEAR MODULUS IN THE 3-1 PLANE 

MEAN COEFFICIENT OF THERMAL EXPANSION IN THE 1 DIRECTION 
MEAN COEFFICIENT OF THERMAL EXPANSION IN THE 2 DIRECTION 
MEAN COEFFICIENT OF THERMAL EXPANSION IN THE 3 DIRECTION 


3920C- 

3930C 

3940C 

3950C 

3960C 

3970C 

3980C 

3990C 

4000C 

40 IOC 
4020C 
4030C 
404 OC 
4050C 
4060C 
4070C 
4080C 
4090C 
4 100C 

4 1 IOC 1 
4 120C 
4 130C 
4 140C 
4 150C 
4 160C 
4 170C 
4 180C 
4 190C 
4200C 

42 IOC 
4220C 
4230C 
4240C 
4 250C 
4260C 
4270C 
42800 
4290C 
4300C 

43 IOC 
4320C 
4330C 
4340C 


V . 2 ORTHOTROPIC AXIS ORIENTATION TABLE 
LINE NOR 

NOR : NUMBER OF ORIENTATION SPECIFICATIONS 

(O.LE.NOR.LE. 10) i.e. MAX OF 10 ORIENTATION SYSTEMS ALLOWED 
(NOTE; IF THE MATERIALS ARE ISOTROPIC OR THE ORTHOTROPIC AXES 
COINCIDE WITH THE GLOBAL AXES. ENTER 0) 

ENTER 0. ( MAXIMUM OF 10 SPECIFICATIONS ALLOWED ) 

IF(NOR.GT.O) 

( ENTER NOR LINES OF THE FOLLOWING ) 

LINE I NO NP NO 

I : ORIENTATION IDENTIFIER ( IOR IN ELEMENT INPUT ) 

NO, NP, NO : NODE NUMBERS IDENTIFYING ORTHOTROPIC AXES ( SEC. II. 2 ) 

V . 3 INELASTIC MATERIAL PROPERTIES 

(IF NPCI . EO . 2 . OR . NPCI . EO . 1 2 ) 

(IF LAWCRP . EO . 1 . OR . LAWCRP . EO . 12) 

(ENTER MTN OF THE NEXT FOUR LINE SETS - MTN , SSTEMP . PPV ) 

LINE MTN NPTS NTM 


MTN 

NPTS 

NTM 


MATERIAL NUMBER 

NUMBER OF STRESS-STRAIN POINTS PER CURVE 

NUMBER OF TEMPERATURES FOR WHICH STRESS-STRAIN CURVES ARE GIVEN 


LINE SSTEMP(I) SSTEMP(2) SSTEMP ( NTM ) 


SSTEMP 


ARRAY OF TEMPERATURES WHERE STRESS-STRAIN CURVES ARE GIVEN 
(IN INCREASING ORDER) 


LINE PPV( MTN ,1,1,1) PPV( MTN ,1,1,2) 
LINE PPV ( MTN ,1,2,1) PPV( MTN , I , 2 . 2 ) 
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TITLE* INPUT FOR T3CYAN WITH SIMPLE 


DATE *10/16/85 


TIME 


10.28 


PAGE* 


4350C 
4 360C 
4 370C 
4380C 
4390C 
44000 

44 IOC 
4420C 
4430C 
444QC 
4450C 
4460C 
4470C 
4480C 
4490C 
4500C 

45 IOC 
4520C 
4530C 
4540C 
4550C 
4560C 
4570C 
4580C 
4590C 
4600C- 
4€ IOC 
4620C 
4630C 
4640C 
4650C 
4660C 
4670C 
4680C 
4690C 
4700C 

47 IOC 
4720C 
4730C 
4740C 
47SOC 
4760C 
4770C 
4780C 
4790C 
4800C 

48 IOC 
4820C 
4830C 
4840C 
4850C 
4860C- 


LINE PPV ( MTN , I , NPTS , 1 ) PPV( MTN , I , NPTS , 2 ) 
PPV ( MTN , I f J , 1 ) 


STRESS VALUES IN INCREASING OROER FROM j-i TO NPTS 
DDw ( utm t , £2* TEMPERATURE I ANO MATERIAL MTN 

PPV ( MTN , I , J f 2 ) : STRAIN VALUES IN INCREASING ORDER FROM J-i TO NPTS 
FOR TEMPERATURE I AND MATERIAL MTN 

(IF LAWCRP.E0.2.0R.LAWCRP.E0. 12) 

(ENTER MTN OF THE FOLLOWING TWO LINE SETS) 

LINE NCTP TCUT SNORM 


NCTP 

TCUT 

SNORM 


NUMBER OF TEMPERATURES WHERE CREEP PROPERTIES SPECIFIED 

CUTOFF TEMPERATURE BELOW WHICH NO CREEP OCCURS 

NORMALIZING STRESS (IN PSI ) FOR THE FOLLOWING CREEP PROPERTIES 




(ENTER NCTP OF THE FOLLOWING LINE) 
LINE TEMP 0 R STRCUT 


TEMP 

0 

R 

STRCUT 


TEMPERATURE WHERE THESE CREEP PROPERTIES APPLY 
CREEP PROPERTY 

CREEP PROPERTY (CREEP STRAIN • 0 * STRESS •• R) 

: CUTOFF STRESS (IN PSI) BELOW WHICH NO CREEP OCCURS 


VI TIME AND TIME INCREMENTING CONTROL INPUT 

LINE N2M TCRP TINIT ECMAX SIGMAX ERMAX DELMIN DELMUL 

N2M : NUMBER OF EQUAL TIME STEPS IN THIS LOAO CASE IF N2M-0 
DYNAMIC TIME INCREMENTING WILL BE USED 
TCRP : TOTAL TIME IN THIS LOAD CASE ( SECONDS ) 

Tl " 1 Z = ^I 1 ! 1 TIME STEP * IF THIS IS NOT THE FIRST LOAD CASE AND 
0.0 IS INPUT, .5 TIMES THE LAST CALCULATED TIME STEP OF THE 
PREVIOUS LOAD CASE IS USED. IF THIS IS THE FIRST LOAD CASE OR 

DELMIN S IS S USED FIRST CASE ° F * RESTART ANID °-° « INPUT. 

EC default A SI!!ue Js EL ^i$o STRAIN INCREMENT DESIRED IN ANV TIME step. 
SIGMAX : MAXIMUM CHANGE IN STRESS DESIRED IN ANY TIME STEP 
DEFAULT VALUE IS 1000 PSI . 

^DEF AULT^I S^O EE J^ A ^* E ^ INTEGRATION ERROR ALLOWED IN ANY TIME STEP 

DE !f£pp : MJ NINUM ALLOWABLE TIME STEP. DEFAULT VALUE IS .001 TIMES 

DELMUL : MAXIMUM MULTIPLIER ALLOWED ON CURRENT TIME STEP IN 
CALCULATING THE NEXT TIME STEP. DEFAULT ■ 1.5 


4870C 

4880C 

4890C 

490OC 

49 IOC 
4920C 
4930C 
4940C 
4950C 
4960C 
4970C 
4980C 
4990C 
5000C 

50 IOC 


VII CONVERGENCE CRITERIA 
LINE I DM ICON DELSIG DELEPS 


I DM 
ICON 


DELSIG 

DELEPS 


MATERIAL NUMBER 
CONVERGENCE CRITERIA CODE 

* 1 CHANGE IN EFFECTIVE STRESS FOR SUBSEQUENT ITERATIONS 

MUST BE LESS THAN DELSIG 

■ 2 CHANGE IN EFFECTIVE INELASTIC STRAIN INCREMENT FOR 

SUBSEQUENT ITERATIONS MUST BE LESS THAN DELEPS DELEPS 

* 3 BOTH CONDITIONS 1 AND 2 MUST BE SATISFIED, DEFAULT - 3 
CONVERGENCE TOLERANCE ON EFFECTIVE STRESS. 

DEFAULT * . 01 "SIGMAX 

CONVERGENCE TOLERANCE ON EFFECTIVE INELASTIC STRAIN INCREMENT 
DEFAULT * . 01 *£CMAX 


n 
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TITLE* INPUT FOR T3CYAN WITH SIMPLE 


DATE *10/16/85 


TIME * 10.28 


PAGE* 7 



r 
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; 5020C 

I 5030C 

! 5040C - 

5050C VIII INITIAL CONSTRAINED DISPLACEMENTS 
5060C 

i , 5070C LINE NDB 

| ' 5080C 

5090C NDB * NUMBER OF CONSTRAINED DISPLACEMENT SPECIFICATIONS 
5100C 

r~ ■ 5 1 IOC ( ENTER NOB OF THE FOLLOWING LINES ) 

S120C 

! 5130C LINE N IDIR VALUE NEND NINC 

5140C 

5150C N : NODE NUMBER 

r 5160C IDIR : DIRECTION CONSTRAINED 

5 1 TOC » 1 X CONSTRAINED ( ANY COMBINATION OF COOES, 

5180C * 2 Y CONSTRAINED 1.«. 12, 13, OR 123 

5190C * 3 2 CONSTRAINED MAY ALSO BE USED) 

5200C VALUE : NUMERICAL VALUE OF CONSTRAINED DISPLACEMENT. DEFAULT IS 0.0 

j 52 IOC NEND : LAST NODE NUMBER HAVING THIS CONSTRAINT. IF OMITTED N IS ASSUMED 

5220C NINC : INCREMENT TO BE USED FOR CONSTRAINT GENERATION FROM N TO NEND 
'• 5230C DEFAULT IS 1 . ( OPTIONAL ) 

5240C 

r 5250C 

• 5260C 

■ S270C IX INITIAL NODAL APPLIED FORCES 

5280C 

5290C LINE NFB 
f~ 5300C 

I 53 IOC NFB : NUMBER OF APPLIED NODAL FORCE SPECIFICATIONS 

I 5320C 

5330C ( ENTER NFB OF THE FOLLOWING LINES ) 

5340C 

j 5350C LINE N IDIR VALUE NEND NINC 

5360C 

1 5370C N : NODE NUMBER 

5380C IDIR : DIRECTION OF APPLIED FORCE 
r ' 5390C • 1 X DIRECTION 

j 5400C * 2 Y DIRECTION 

! 54 IOC *32 DIRECTION 

5420C VALUE : NUMERICAL VALUE OF FORCE TO BE APPLIED ( LBS ) 

5430C NEND : LAST NODE HAVING THIS APPLIED FORCE. IF OMITTED N IS ASSUMED 

r 5440C NINC : INCREMENT TO BE USED FOR FORCE GENERATION FROM N TO NEND. 

! 5450C DEFAULT IS 1 . ( OPTIONAL ) 

1 ' 5460C 

5470C 

5480C 

J 5490C X INITIAL NODAL TEMPERATURES 

I 5500C 

55 IOC LINE NTEMPS 
5520C 

i— 5530C NTEMPS : NUMBER OF TEMPERATURE INPUT LINES 

| 5540C 

! 5550C ( ENTER NTEMPS OF THE FOLLOWING LINES ) 

5560C 

5570C LINE N TEMP ITYPE NENO NINC 
f 5580C 

f 55900 N : NODE NUMBER 

1 5600C TEMP : TEMPERATURE ( DEGREES F. ) 

56 IOC ITYPE : INDICATOR FOR TYPE OF TEMPERATURE REVISION 

r 5620C * 0 CHANGE NODAL TEMPERATURE TO TEMP 

j S630C ■ 1 INCREMENT NOOAL TEMPERATURE BY TEMP 

I S640C NEND : LAST NODE HAVING THE GIVEN TEMPERATURE. IF OMITTED N IS ASSUMED 

5650C NINC : INCREMENT TO BE USED ON TEMPERATURE GENERATION FROM N TO NENO. 
S660C IF OMITTED 1 IS ASSUMED 

r~ 5670C 

i 5680C 
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TITLE* INPUT FOR T3CYAN WITH SIMPLE 


DATE *10/16/85 


TIME • 10.28 


PAGE * 8 


5690C 

5700C XI INITIAL ELEMENT PRESSURE LOAOS 
57 IOC 

5720C LINE NPL 
5730C 

5740C NPL : NUMBER OF PRESSURE INPUT LINES 
5750C 

5760C ( ENTER NPL OF THE FOLLOWING LINES ) 

5770C 

5780C LINE NEL IFACE PI P2 P3 P4 PS P6 


' ! 


5790C 


5800C 

NEL 

58 IOC 

IFACE 

5820C 

PI 

5830C 

P2. P3 

5840C 


5850 


S860C 

NEND 

5870C 


5880C 

NINC 


P7 P8 NEND NINC 


ELEMENT NUMBER 

FACE NUMBER ( IF ACE . GE . 1 . AND . I F ACE . LE . 6 ) 

PRESSURE AT NODE 1 ( PSI ) 

P4, P5 , P6. P7 , P8 : PRESSURES AT NODES 2 8 ON THE FACE. 

IF OMITTED THEY ARE SET EQUAL TO PI. 
(POSITIVE PRESSURES INDUCE COMPRESSION IN THE ELEMENT ) 

LAST ELEMENT HAVING THIS PRESSURE LOADING. IF OMITTED NEL IS 
ASSUMED 

INCREMENT TO BE USED FOR ELEMENT PRESSURE GENERATION FROM 
NEL TO NEND. IF OMITTED NEL IS ASSUMED 


'"I 


XII. 1 LOAD CASE CONTROL CARD 
RPM IAXIS I ACC NT I 

■aawaaaaaaa 

‘ION SPECIF 
ACCELX ACCELV ACCEL2 


5890C 
5900C 
59 IOC 

5920C 

5930C XII LOAD CASE INFORMATION. FINAL CONDITIONS 
5940C 

5 950C* - •••**" , **»*«**** ■ ■«•••••••*»•*•• mmmmmmmmmmmm 

5960C 
5970C 

5980C LINE 
5990C 

6000C — *-*- 

6010C XII. 2 ACCELERATION SPECIFICATIONS FOR INERTIA OR GRAVITY LOADS 
6020C 

6030C LINE 
6040C- 
S050C 

6060C 

6070C XIII NODAL CONSTRAINED DISPLACEMENTS 
6080C 

6090C LINE NDB 
€ 1 0OC 

61 IOC ( ENTER NOB OF THE FOLLOWING LINES ) 

6120C 

6130C LINE N IDIR VALUE NEND NINC 
6140C 
6 150C 

6160C — 

6170C XIV NODAL APPLIED FORCES 
6180C 

6190C LINE NFB 
6200C 

62 IOC ( ENTER NFB OF THE FOLLOWING LINES ) 

6220C 

6230C LINE N IDIR VALUE NEND NINC 

6240C 
6250C 

6260C - 

6270C XV NODAL TEMPERATURES 
6280C 

6290C LINE NTEMPS 
6300C 

63 IOC ( ENTER NTEMPS OF THE FOLLOWING LINES ) 

6320C 

6330C LINE N TEMP ITYPE NEND NINC 

6340C 
6350C 




n 


n 

n 
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TITLE* INPUT FOR T3CYAN WITH SIMPLE 


DATE « 10/ 16/85 


TIME * 10.28 


PAGE* 


6360C 


6370C XVI ELEMENT PRESSURE LOADS 
6380C 

6390C LINE NPL 
6 4000 

64 IOC ( ENTER NPL OF THE FOLLOWING LINES ) 
6420C 

6430C LINE NEL I FACE PI P2 P3 P4 P5 


P6 


P7 


P8 NEND NINC 
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TITLE- INPUT FOR T3CYAN WITH SOONER DATE -10/16/85 


ORGANIZATION OF INPUT 

I HEAOING AND CONTROL INFORMATION 
I • 1 title CARO 

1 . 2 PROBLEM SIZING 

1.3 ANALYSIS AND RESTART OPTIONS 
14 EQUATION NUMBERING AND BANOING OPTIONS 

II NODE COORDINATES AND TRANSFORMATIONS 

11. 1 NODE COORDINATES 

1 1. 2 LOCAL NODE COORDINATE SYSTEM TRANSFORMATIONS 

III ELEMENT DEFINITION 

111.1 HEADER LINE FOR ELEMENT 

111. 2 20 NODED SOLID DEFINITION 

IV LOAD CASE INFORMATION, INITIAL CONDITIONS 

IV. 1 LOAD CASE CONTROL CARO 

IV. 2 ACCELERATION SPECIFICATIONS 

V MATERIAL PHYSICAL PROPERTIES 

V. 1 MATERIAL PHYSICAL PROPERTIES 

V.1.1 ISOTROPIC ELASTIC PROPERTIES 
V.1.2 ORTHOTROPIC ELASTIC PROPERTIES 

V • 2 ORTHOTROPIC AXES ORIENTATION TABLE 
V • 3 INELASTIC MATERIAL CHARACTERIZATION 

VI TIME AND TIME INCREMENTING CONTROL 

1290C VII CONVERGENCE CRITERIA 
1300C 

1310C VIII INITIAL NODAL CONSTRAINED DISPLACEMENTS 
1320C 

1330C IX INITIAL NODAL APPLIED FORCES 
1340C 

1350C X INITIAL NODAL TEMPERATURES 
1 360C 

1370C XI INITIAL ELEMENT PRESSURE LOADS 
1380C 

1390C XII LOAD CASE INFORMATION. FINAL CONDITIONS 
1400C XII. 1 LOAD CASE CONTROL CARD 

14100 XII. 2 ACCELERATION SPECIFICATION 

1420C 

1430C XIII FINAL NODAL CONSTRAINED DISPLACEMENTS 
1 44QC 

1450C XIV FINAL NODAL APPLIED FORCES 
1460C 

1470C XV FINAL NODAL TEMPERATURES 
1480C 

1490C XVI FINAL ELEMENT PRESSURE LOADS 

1500C 

15 IOC 

1520 END 


1000C 

1010C 

1020C 

1030C 

1040C 

1050C 

1060C 

1070C 

1080C 

1090C 

I 100C 

II IOC 
1 120C 
1 130C 
1 140C 
1 1 50C 
1 160C 
1 170C 
1 180C 
1 1 90C 
120OC 
1 2 IOC 
1220C 
1230C 
1240C 
1250C 
1 260C 


TIME 


10.28 


PAGE- 


n 




n 

( 
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TITLE- INPUT FOP T3CYAN WITH BODNER DATE -10/16/85 TIME • 10.28 

1 5300 I HE AO I NG AND CONTROL INFORMATION 
1540C 

1550C I . 1 TITLE CARD 
1560C 

1570C LINE ITITLE 

1580C ITITLE - ANY 1 TO 72 CHARACTER TITLE FOR THE ANALYSIS 
1590C 

1600C 1.2 PROBLEM SIZING 

16 IOC 

1620C LINE NUMNP NM IT RTEM NLC 

1630C NUMNP • NUMBER OP STRUCTURAL NODES ( ENTER AS A NEGATIVE NUMBER FOR 
1640C TIMING SUMMARY ) 

1650C NM . NUMBER OF DIFFERENT MATERIALS ( MAXIMUM-3 ) 

1660C IT • THERMAL STRESS OPTION 

1670C « 0 INCLUDE THERMAL LOAOS 

1680C - 1 IGNORE THERMAL LOADS 

1690C RTEM • REFERENCE TEMPERATURE ( DEGREES F ) 

1700C 

17 IOC 1.3 ANALYSIS AND RESTART OPTIONS 

1720C 

1730C LINE LAWCRP NOUT NR E ST A INREST 
1740C LAWCRP * TYPE OF INELASTIC ANALYSIS 
1750C » 0 ELASTIC ANALYSIS 

1760C ■ 1 ISOTHERMAL SOONER MODEL 

1770C NOUT • OUTPUT FILE CREATION OPTION 
1780C * 0 DO NOT CREATE OUTPUT FILE 

1790C - 1 CREATE OUTPUT FILE 

180QC NRESTA » RESTART OPTION 

18 IOC - 0 THIS IS NOT A RESTART RUN 

1820C > 0 INPUT THE LOAD CASE FROM WHICH THE RESTART IS TO PROCEED 

1830C ( OUTPUT FROM THIS CASE MUST HAVE BEEN PREVIOUSLY PUT ON AN OUTPUT 

1840C FILE ). THE FIRST NEW LOAD CASE WILL BE LABELED AS NRESTA + 1. 

18S0C INREST ■ LOAD CASE NUMBER IN THE CURRENT INPUT FILE WHICH BECOMES 

1860C THE FIRST NEW LOAO CASE TO BE SOLVED WHEN RESTARTING. WHERE 

1870C ( 1 . LE . INREST . LE .NLC ). IF 0 IS INPUT INREST « 1 IS ASSUMED 

1880C 

1890C 1.4 EOUATION NUMBERING AND BANOING OPTIONS 

1900C 

1 9 IOC LINE N IBAND IPBAND 

1920C N • KEY CODE, ENTER -1 TO ACTIVATE THIS OPTION 
1930C IBAND ■ BANDING OPTION 

1940C 0 - USE DEFAULT OPTION 

1950C 1 « ASSUME NODE NUMBER IS THE SAME AS MATRIX POSITION 

1960C 2 ■ ASSUME INPUT NODE ORDER DEFINES MATRIX POSITION 

1970C IPBAND » PRINTOUT OPTION 
1980C 0 ■ NO PRINTOUT OF EOUATION NUMBERS 

1990C 1 « PRINT OUT THE EOUATION NUMBER FOR EACH DEGREE OF FREEDOM 

2000C 

20 IOC 

2020 END 
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DATE -10/16/85 


TIME • 10.28 


PAGE- 3 
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n 
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2030C 
20400 
2050C 
2060C 
2070C 
2080C 
2090C 
2 100C 

21 IOC 
2120C 
2130C 
2140C 
2150C 
2160C 
2170C 
2180C 
2190C 
2200C 

22 IOC 
2220C 
2230C 
2240C 
2250C 
2260C 
2270C 
2280C 
2290C 
2300C 

23 IOC 
2320C 
2330C 
2340C 
2350C 
2360C 
2370C 
2380C 
2390 


II NODE COORDINATES AND TRANSFORMATIONS 

11. 1 NODE COORDINATES 

LINE N X Y 2 
N - NODE NUMBER 

X « NODE X COORDINATE 

v - NODE Y COORDINATE 

2 • NODE 2 COORDINATE 

AFTER THE LAST NOOE ENTER THE FOLLOWING LINE 
LINE 0 NSKEW 

NSKEW - INDICATOR FOR LOCAL NOOE COORDINATE SYSTEM TRANSFORMATION 

0 • NONE WILL BE INPUT 

1 * ONE OR MORE WILL BE INPUT 

1 1.2 LOCAL NODE COORDINATE SYSTEM TRANSFORMATIONS 

( IF NSKEW > 0. ENTER ALL NODES WITH LOCAL COORDINATE SYSTEM 
TRANSFORMATIONS ) 

LINE N NO NI NJ 

N - NODE NUMBER HAVING A ROTATED LOCAL COORDINATE SYSTEM 
NO » NODE ON LOCAL X AXIS 

NI • NODE ON LOCAL X AXIS IN THE +X DIRECTION FROM NO 

NU » NODE IN LOCAL XY PLANE SUCH THAT THE +2 AXIS IS IN THE 

DIREC i ION OF THE VECTOR P3 « (NO-NI) X (NO-NU). THE +Y AXIS IS 
IN THE DIRECTION OF THE VECTOR P2 - P3 X (NO-NI) 

( AFTER THE LAST NOOE HAVING A LOCAL COORDINATE SYSTEM ENTER A LINE 
WITH N ■ 0 ) I . E , 

LINE 0 


END 


n 


i 

i 


n 


i 


I i 


1 



) 



I 


! 


n 
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TITLE* INPUT FOR T3CYAN WITH BODNER DATE *10/16/85 TIME « 10.28 

2400C III ELEMENT DEFINITION 
24 IOC 

2420C 1 1 1 - 1 HEADER LINE FOR ELEMENT OATA 

24 30C 

2440C LINE IPRINT ( 2 ) 

2450C IPRINT ( 2 ) « ELEMENT DATA EXTENOED PRINTOUT OPTION 
2460C 0 « NONE 

2470C 1 * DUMP VOLUMES AND DISTORTION PARAMETERS 

2480C 2 « ALSO DUMP ELEMENT STIFFNESSES 

2490C 3 * ALSO DUMP ELEMENT EQUIVALENT NODAL FORCES AND MASSES 

2S00C 

25 IOC 1 1 1. 2 20 NODED SOLID DEFINITION 
2520C 

2530C LINE NEL N1 N2 N3 N4 N5 N6 N7 N8 N9 N10 Nil N12 

2540C LINE N 1 3 N14 N15 N16 N17 N18 N19 N20 IMAT IOR 

2550C LINE NO NP NO ( OPTIONAL LINE. ENTER ONLY IF IOR • 10O) 

2560C NEL * ELEMENT NUMBER 

2570C N 1 ... N20 * NODES DEFINING THE ELEMENT ( SEE FIG III. 2 ) 

2580C IMAT « MATERIAL NUMBER 

2S90C IOR » ORTHOTROPIC MATERIAL AXIS SYSTEM INDICATOR 
2600C IOR * 0 ISOTROPIC MATERIAL 

26 IOC IOR .GE . 1 . AND . IOR . LE .99 ORTHOTROPIC AXES SPECIFIED BY THE TABLE 

2620C GIVEN UNDER MATERIAL DATA SECTION V.2 

2630C IOR- 100 ORTHOTROPIC AXES DEFINED BY OPTIONAL LINE OF INPUT 

2640C 

2650C ( AFTER THE LAST ELEMENT. ENTER THE FOLLOWING LINE ) 

2660C 

2670C LINE 0 

2680C 

2690C 

2700 END 
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DATE *10/16/85 


TIME * 10.28 


IV LOAD CASE INFORMATION. INITIAL CONDITIONS 


IV. 1 LOAD CASE CONTROL CARD 

LINE RPM IAXIS IACC NPCI NT I ITHST 
RPM « ROTATIONAL SPEED IN RPM 

IAXIS » ?v?e AL AXIS AB0UT WHICH STRUCTURE IS ROTATING 

1 * X AXIS 

2 » y AXIS 

3 « Z AXIS 

IA 0 C **NO N INPUT° R F0R 1NPUT 0F ACCELERAT0N tOADS IN GROUP IV. 2 
1 * INPUT 

NPCI * INDICATOR FOR A CHANGE IN MATERIAL PROPERTIES 

0 • NO CHANGE TO MATERIAL PROPERTIES OP6RTI *5 

1 « CHNAGE ELASTIC PROPERTIES ( GROUPS V.1 AND V 2) 

2 * CHANGE INELASTIC PROPERTIES ( GROUP V.3 ) 

12 * CHANGE ELASTIC AND INELASTIC PROPERTIES 

NTI - NUMBER OF NODAL TEMPERATURES WHICH ARE RESPECIFIED 

0 * NO RESPECIFIED TEMPERATURES BUT STIFFNESS IS RECOMPUTED 

■’* N ? n ??p P !S^^ D 0 I EMPERATUR£S AND STIFft NESS IS NOT RECOMPUTED 
ITHST * TOTAL STRAIN PRINTOUT OPTION 

0 ■ INCLUDE THERMAL STRAINS IN TOTAL STRAIN PRINTOUT 

1 * DO NOT INCLUDE THERMAL STRAINS IN TOTAL STRAIN PRINTOUT 

IV. 2 ACCELERATION SPECIFICATION FOR INERTIAL OR GRAVITY LOADS 
( ENTER ONLY IF IACC > 0 ) 

LINE ACCELX ACCELY ACCEL2 

*^ELX - ACCELERATION OF STRUCTURE ( IN/SEC-2 ) IN GLOBAL X 

ACCELY ■ ACCELERATION OF STRUCTURE ( IN/SEC — 2 ) IN GLOBAL Y 

ACCELZ - ACCELERATION OF STRUCTURE ( IN/SEC-2 ) IN GLOBAL Z 


-- T30RG 


158 


date *10/16/85 


r -■ 
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TITLE* INPUT FOR T3CYAN WITH BODNER 

V MATERIAL PHYSICAL PROPERTIES 

V.1 ELASTIC CONSTANTS ( ENTER NM GROUPS ) 


TIME » 10.28 


3070C 
3080C 
3090C 
3 100C 
31 IOC 
3 120C 
3130C 
3 140C 
3 150C 
3160C 
3170C 
3 180C 


3 190C 

V. 1 

3200C 

32 IOC 

( IF 

3220C 

3230C 

LINE 

3240C 

TEMP 

3250C 

E « 

3260C 

PR - 

3270C 

AL ■ 

3280C 

3290C 

V. 1 

3300C 

33 IOC 

( IF 1 

3320C 

CORREJ 

3330C 

3340C 

LINE Tl 

3350C 

TEMP 

3360C 

Ell - 

3370C 

E22 > 

3380C 

£33 * 

3390C 

NU 1 2 

3400C 

NU 1 3 

3410C 

NU23 

3420C 

G 1 2 ■ 

3430C 

G23 ■ 

3440C 

G31 * 

3450C 

AL 1 * 

3460C 

AL2 - 

3470C 

A L3 • 

3480C 

3490C 

V.2 

3500C 


LINE MTN NMT DEN 

MTN • MATERIAL NUMBER ( MTN . GE . 1 . ANO . MTN . LE . NM ) 

SIMPLIFY R INPUT) IS0TR0PIC * INPUT MTN AS A NEGATIVE NUMBER TO 

NM Gi:E^OR R T SlS T MrTERllL RES ^ ^ ELASTIC PR0PERTI “ WILL BE 
DEN « WEIGHT DENSITY OF THE MATERIAL ( P0UNDS/IN**3 ) 

1 ISOTROPIC MATERIAL 

( IF MTN < 0/ ENTER NMT OF THESE LINES ) 


TEMP E PR AL 

TEMPERATURE ( DEGREES F. 


) 


10**6 P . S. I . ) 


ELASTIC MODULUS ( 

POISSON'S RATIO 
MEAN COEFFICIENT OF THERMAL EXPANSION ( 

.2 ORTHOTROPIC MATERIAL 


10* *~6 IN/IN-DEG. F. ) 


0, ENTER NMT OF THESE LINES. THE DIRECTIONS 1,2.3 


X' ,Y' ,Z' 

MP Ell E22 E33 NU12 NU13 NU23 G12 G23 G3 1 AL 1 
■ TEMPERATURE AT WHICH PROPERTIES ARE GIVEN ( 
ELASTIC MODULUS IN THE 1 DIRECTION 
ELASTIC MODULUS IN THE 2 DIRECTION 
ELASTIC MODULUS IN THE 3 DIRECTION 

* POISSON'S RATIO RELATING DIRECTIONS 

* POISSON'S RATIO RELATING DIRECTIONS 

* POISSON'S RATIO RELATING DIRECTIONS 
SHEAR MODULUS IN THE 1-2 PLANE 
SHEAR MODULUS IN THE 2-3 PLANE 
SHEAR MODULUS IN THE 3-1 PLANE 

MEAN COEFFICIENT OF THERMAL EXPANSION IN THE 

MEAN COEFFICIENT OF THERMAL EXPANSION IN THE 

MEAN COEFFICIENT OF THERMAL EXPANSION IN THE 


AL2 AL3 
OEG . F. ) 


AND 

AND 

AND 


DIRECTION 

DIRECTION 

DIRECTION 


V.2 ORTHOTROPIC AXIS ORIENTATION TABLE 


35 IOC 
3520C 
3530C 
3540C 
3550C 
3560C 
3570C 
3580C 
3590C 
3600C 
3610C 
3620C 
3630C 
3640C 
3650C 
3660C 
3670C 
3680C 
3690C 
3700C 
37 IOC 
3720C 
3730C 


LINE NOR 

NOR * NUMBER OF ORIENTATION SPECIFICATIONS, IF THE MATERIALS ARE 
ISOTROPIC OR THE ORTHOTROPIC AXES COINCIDE WITH THE GLOBAL AXES 
ENTER 0. ( MAXIMUM OF 10 SPECIFICATIONS ALLOWED ) 


( ENTER NOR LINES OF THE FOLLOWING ) 


LINE I NO NP NO 

I • ORIENTATION IDENTIFIER ( I OR IN ELEMENT INPUT ) 

NO. NP, NO - NODE NUMBERS IDENTIFYING ORTHOTROPIC AXES ( SEC. II. 2 ) 

V.3 INELASTIC MATERIAL PROPERTIES 

LINE MTN NCTEM 

MTN - MATERIAL NUMBER 

NCTEM * NUMBER OF TEMPERATURES AT WHICH INELASTIC MATERIAL PROPERTIES 
WILL BE GIVEN FOR THIS MATERIAL 

( ENTER NCTEM OF THE FOLLOWING LINES ) 

LINE TEMP 0 AN 20 Z1 22 AM A R 
TEMP • TEMPERATURE IN DEGREES F 

0. AN, 20, 21. Z2 , AM, A, R • MATERIAL PARAMETERS REQUIRED FOR BODNER' 
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n 


TITLE. INPUT FOR T3CYAN WITH BOONER DATE -10/16/85 TIME ■ 10.28 PAGE ■ 7 

37400 INELASTIC CONSTITUTIVE MODEL AT THIS TEMPERATURE 

3760C 

3770 END 


n 






n 


n 


n 



r*\ 


i 

i 
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TITLE* INPUT for T3CYAN WITH SOONER 


DATE *10/16/85 


TIME * 10.28 


VI TIME AND TIME INCREMENTING CONTROL INPUT 


LINE N2M TCRP UNIT ECMAX SIGMAX ERMAX DELMIN DELMUL 

N2M • NUMBER OF EQUAL TIME STEPS IN THIS LOAD CASE, IF N2M-0 
DYNAMIC TIME INCREMENTING WILL BE USED 

TCRP » TOTAL TIME IN THIS LOAD CASE ( SECONDS ). 

TINIT • INITIAL TIME STEP. IF THIS IS NOT THE FIRST LOAD CASE AND 
0.0 IS INPUT, .5 TIMES THE LAST CALCULATED TIME STEP OF THE 
PREVIOUS LOAD CASE IS USED. IF THIS IS THE FIRST LOAD CASE OR 
IF THIS IS THE FIRST LOAD CASE OF A RESTART AND 0.0 IS INPUT 
DELMIN IS USED. 

ECMAX * MAXIMUM INELASTIC STRAIN INCREMENT DESIRED IN ANY TIME STEP 
DEFAULT VALUE IS .000100. 

SIGMAX - MAXIMUM CHANGE IN STRESS DESIRED IN ANY TIME STEP 
DEFAULT VALUE IS 1000 PSI . 

ERMAX « MAXIMUM ESTIMATED INTEGRATION ERROR ALLOWED IN ANY TIME STEP 
DEFAULT IS .01 ( 1% ). 

DELMIN - MINIMUM ALLOWABLE TIME STEP . DEFAULT VALUE IS .001 TIMES 
TCRP . 

DELMUL - MAXIMUM MULTIPLIER ALLOWED ON CURRENT TIME STEP IN 
CALCULATING THE NEXT TIME STEP. DEFAULT - 1.5 
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TITLE- INPUT FOR T3CYAN WITH SOONER 


DATE -10/16/05 


TIME • 10.28 


PAGE- 


40200 
403OC 
4040C 
4050C 
4060C 
4070C 
408OC 
4090C 
4 1000 
4 1 IOC 
4 1200 
4 1300 
4 1400 
4 1500 
4 160C 
4 170 


VII CONVERGENCE CRITERIA 


DELEP$ 


LINE I DM ICON OELSIG 
IDM « MATERIAL NUMBER 
ICON - CONVERGENCE CRITERIA CODE 

* 1 CHANGE IN EFFECTIVE SRESS FOR SUBSEQUENT ITERATIONS < DELSIG 
ITERATIONS ^DE^EPS 7 * ^ INELASTIC STRAIN INCREMENT FOR SUBSEQUENT 

* 3 BOTH CONDITIONS 1 AND 2 MUST BE SATISFIED DEFAULT « 3 
DELSIG » CONVERGENCE TOLERANCE ON EFFECTIVE STRESS 

DEFAULT - . 01 *SIGMAX 

0E DEFlulT C ° NV "Sx T0LeRANCE ° N EFFECTIVE ^ELASTIC STRAIN INCREMENT 


END 
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DATE *10/16/85 


TIME 


10.28 


page • 10 


4 180C 
4 190C 
4 200C 
42 IOC 
4220C 
4230C 
4240C 
4 250C 
4 260C 
4 2 70C 
4280C 
4290C 
4 300C 
43 IOC 
4320C 
4 330C 
4340C 
4350C 
4360C 
4370 


VIII INITIAL CONSTRAINED DISPLACEMENTS 
LINE NDB 

NDB • NUMBER OF CONSTRAINED DISPLACEMENT SPECIFICATIONS 

( ENTER NDB OF TH FOLLOWING LINES ) 

LINE N ID I R VALUE NEND NINC 
N * NODE NUMBER 
ID IR * DIRECTION CONSTRAINED 

1 * x (ANY COMBINATION OF CODES, l.e. 12. 13 OR 123 MAY 

2 * Y ALSO BE USED) 

3 * Z 

VALUE - NUMERICAL VALUE OF CONSTRAINED DISPLACEMENT. DEFAULT IS 0 0 
NEND * LAST NODE NUMBER HAVING THIS CONSTRAINT. IF OMITTED N IS ASSUMED 
NINC * INCREMENT TO BE USED FOR CONSTRAINT GENERATION FROM N TO NEND 
DEFAULT IS 1. ( OPTIONAL ) 


END 
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TITLE* INPUT FOR T3CYAN WITH SOONER 


DATE *10/16/85 


TIME « 10.28 


4380C IX INITIAL NODAL APPLIED FORCES 
4 390C 

4400C LINE NFB 

44 IOC NFB * NUMBER OF APPLIED NODAL FORCE SPECIFICATIONS 
4420C 

4430C ( ENTER NFB OF THE FOLLOWING LINES ) 

4440C 

4450C LINE N ID I R VALUE NEND NINC 
4460C N - NODE NUMBER 

4470C ID I R ■ DIRECTION OF APPLIED FORCE 
4480C 1 - X 

4490C 2 ■ Y 

4500C 3 « Z 

4 5 1 OC 

4520C VALUE ■ NUMERICAL VALUE OF FORCE TO BE APPLIED ( LBS ) 

4530C NEND » LAST NODE HAVING THIS APPLIED FORCE. IF OMITTED N I 

4540C NINC - INCREMENT TO BE USED FOR FORCE GENERATION FROM N TO 

4550C DEFAULT IS 1. ( OPTIONAL ) 

4560C 

4570C 

4580 END 


i ASSUMED 
NEND. 
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DATE ■ 10/ 16/85 


TIME 


10.28 


PAGE* 12 


4590C 

4600C 

46 IOC 

4620C 

4630C 

4640C 

4650C 

4660C 

4670C 

4680C 

4690C 

4700C 

47 IOC 

4720C 

4730C 

4740C 

4750C 

4760C 

4770 


X INITIAL NODAL TEMPERATURES 
LINE NTEMPS 

NTEMPS » NUMBER OF TEMPERATURE INPUT LINES 

( ENTER NTEMPS OF THE FOLLOWING LINES ) 

LINE N TEMP I TYPE NEND NINC 
N * NODE NUMBER 

TEMP • TEMPERATURE ( DEGREES F. ) 

ITYPE • INDICATOR FOR TYPE OF TEMPERATURE REVISION 

0 « CHANGE NODAL TEMPERATURE TO TEMP 

1 • INCREMENT NODAL TEMPERATURE BY TEMP 

™ * ^ AST N0DE MAVING the given temperature. IF OMITTED n is assumed 

if om!t^d M i N Is T as“med E0 0N TEMPERATURE generation prom N to NEND 


END 
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TITLE* INPUT FOP T3CYAN WITH BODNER 


DATE -10/16/85 


TIME - 10.28 


4780C 

4 7900 

4800C 

48 IOC 

48200 

4830C 

4840C 

4 8 500 

48600 

48700 

48800 

48900 

4900C 

49100 

49200 

4930C 

4940C 

49500 

49600 

4970C 

4980 


XI INITIAL ELEMENT PRESSURE LOADS 
LINE NPL 

NPL • NUMBER OF PRESSURE INPUT LINES 
( ENTER NPL OF THE FOLLOWING LINES ) 

LINE NEL IFACE Pi P 2 P3 P4 P5 P6 p7 pg NEND NINO 
NEL * ELEMENT NUMBER ° NINC 

IFACE - FACE NUMBER ( I FACE . GE . 1 . AND . I FACE LE 6 ) 

PI » PRESSURE AT NODE 1 ( PSI ) 

P2, P3, P4 , P5, P6 , P7 , PS ■ PRESSURES AT NODES 2 8 ON THF fa cf 

IF OMITTED THEY ARE SET EQUAL TO Pi . ( POSITIVE PRESSURES INDUCE 
COMPRESSION IN THE ELEMENT ) iVE PRESSURES INDUCE 

NE assumeS ST element having this pressure loading, if OMITTED NEL IS 

NI NFL*Tn N wc5? EN Ie T SufT T ^ £D F ° R ELEM ENT PRESSURE GENERATION FROM 
NEL TO NEND. IF OMITTED NEL IS ASSUMED 


END 
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DATE -10/16/85 
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PAGE- 


49900 XII LOAD CASE INFORMATION, FINAL CONDITIONS 
5000C 

5010C XII. 1 LOAO CASE CONTROL CARD 
5020C 

5030C LINE RPM IAXIS IACC NT I 
5040C 

5060C XI1 - 2 ACCELERATION SPECIFICATIONS FOR INERTIA OR GRAVITY LOADS 

5070C LINE ACCELX ACCELY ACCELZ 

5080C 

5090C 

5100 END 
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TITLE. INPUT FOR T3CYAN WITH BOONER DATE .10/16/85 TIME - 10.28 PAGE-1S 

SI IOC XIII NODAL CONSTRAINED DISPLACEMENTS 
5 1 20C 

5130C LINE NOB 
5 140C 

5150C ( ENTER NDB OF THE FOLLOWING LINES ) -1 

5160C ; 

5170C LINE N IDIR VALUE NEND NINC ! 

5 1 80C 
5 190C 

5200 END 

■ t 


i 

1 


n 


T30RG 


PAGE- 15 


168 


TITLE* INPUT FOR T3CYAN WITH BODNER 


DATE -10/16/05 


TIME - 10.28 


PAGE- 16 


52100 XIV NODAL APPLIED FORCES 
f 52200 

j 5230C LINE NFB 

i 5240C 

5250C ( ENTER NFB OF THE FOLLOWING LINES ) 
5260C 

5270C LINE N IDIR VALUE NEND NINC 
! 5280C 

1 5290C 

5300 END 

n 


i 


— T30RG 


PAGE* 16 


169 



TITLE- INPUT FOR T3CY4N WITH BOONER DATE -10/16/85 TIME - 10.26 

53 IOC XV NODAL TEMPERATURES 
5320C 

5330C LINE NTEMPS 
5340C 

5350C ( ENTER NTEMPS OF THE FOLLOWING LINES ) 

5 360C 

5370C LINE N TEMP ITYPE NEND NINC 

5380C 

5390C 

5400 ENO 
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